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DOCUMENT- IDENTIFIER: US 6512371 B2 

TITLE: System and method for determining oi 1 , water and gas saturations for 
low-f ipld gradient NMR logging tools 



Abstract Text (1) : 

A well logging system and method are disclosed for detecting the presence and 
estimating the quantity of gaseous and liquid hydrocarbons in the near wellbore 
zone. The system uses a gradient-based, multiple -frequency NMR logging tool to 
extract signal components characteristic for each type of hydrocarbons. To this end, 
new data acquisition methods are proposed in which measurements at different 
frequencies are interleaved to obtain, in a single logging pass, multiple data 
streams corresponding to different recovery times and/or diffusivity for the same 
spot in the formation. The resultant data streams are processed to determine 
mineralogy- independent wahpr and hydrocarbon saturations and porosity estimates. Gas 
and oi 1 saturations are used to obtain accurate estimates of the water content, 
permeability and other parameters of interest. In another aspect, a novel 
diffusion- enhanced data acquisition sequence is disclosed for use with low field 
gradient tools. 

Brief Summary Text (2) : 

The present invention relates to nuclear magnetic resonance (NMR) logging and is 

directed more specifically to a system and method for detecting the presence and 
estimating the quantity of gaseous and liquid hydrocarbons in the near wellbore 
zone. 

ttr-jpf Summary Tpxt (3) : 

Petrophysical parameters "of a geologic formation which are typically used to 
determine whether the formation will produce viable amounts of hydrocarbons include 
the formation porosity PHI,_£lnid saturation S, the volume of the formation, and its 
permeability K. Formation porosity is the pore volume per unit volume of formation; 
it is the fraction of the total volume of a sample that is occupied by pores or 
voids. The saturation S of a formation is the fraction of a its pore volume occupied 
by the f Inid of interest. Thus, wat-pr saturation S.sub.W is the fraction of the pore 
volume which contains wat-pr . The water saturation of a formation can vary from 100% 
to a small val hp which cannot be displaced by <">i T r and is referred to as irreducible 
watPT- saturation S.sub.Wirr. For practical purposes it can be assumed that the oil 
or hydrocarbon saturation of the formation S.sub.O is equal to S.sub.O =1 -S.sub.W. 
Obviously, if the formation's pore space is completely filled with water r that is if 
S.sub.W =1, such a formation is of no interest for the purposes of an oil search. On 
the other hand, if the formation is at S.sub.Wirr it will produce all hydrocarbons 
and no wat.pr . Finally, the permeability K of a formation is a measure of the ease 
with which_fliiida can flow through the formation, i.e., its producibility. 

Brief Summary Text (4) : 

Nuclear magnpt-ic rpgonanrp (NMR \ logging is among the most important methods which 
have been developed to determine these and other parameters of interest for a 
geologic formation and clearly has the potential to become the measurement of choice 
for determining formation porosity. At least in part this is due to the fact that 
unlike nuclear porosity logs, the nmr measurement is environmentally safe and is 
unaffected by variations in matrix mineralogy. The nmr logging method is based on 
the observation that when an assembly of magnetic moments, such as those of hydrogen 
nuclei, are exposed to a static magnetic field they tend to align along the 
direction of the magnetic field, resulting in bulk magnetization. The rate at which 
equilibrium is established in such bulk magnetization upon provision of a static 
magnetic field is characterized by the parameter T.sub.l, known as the spin-lattice 
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rpiavatinn time. Anothe^^related and frequently used_MME^R>gging parameter is the so 
called fipin-Bpin rplayarion time ^atanh T.suh.2 (also known as tTflnSYRrRFi 
-rPlpyai-ion time) which is an expression of the relaxation due to non-homogeneities 
in the local magnetic field over the sensing volume of the logging tool. 

R-ripf Summary Text (5) : 

Another measurement parameter used in NMR well logging is the formation diffusion D. 
Generally, diffusion refers to the motion of atoms in a gaseous or liquid state due 
to their thermal energy. The diffusion parameter D is dependent on the pore sizes of 
the formation and offers much promise as a separate permeability indicator. In an 
uniform magnetic field, diffusion has little effect on the decay rate of the 
measured_JIMR echoes. In a gradient magnetic field, however, diffusion causes atoms 
to move from their original positions to new ones, which moves also cause these 
atoms to acquire a different phase shifts compared to atoms that did not move, and 
will thus contribute to a faster rate of relaxation. Therefore, in a gradient 
magnetic field diffusion is a logging parameter which can provide independent 
information about the structure of the geologic formation of interest, the 
properties of the fin-iris in it, and their interaction. 

W-ri^f Summary Tfvh (6) : 

It has been observed that the mechanisms which determine the value a of T.sub.l, 
t . sub . 2 and D depend on the molecular dynamics of the sample being tested. In bulk 
volume liquids, typically found in large pores of the formation, molecular dynamics 
is a function of molecular size and inter-molecular interactions which are different 
for each fluid. Thus, wat-pr f gas and different types of oil each have different 
T.sub.l, T _ mih - 2 and D values . On the other hand, molecular dynamics in a 
heterogeneous media, such as a porous solid which contains liquid in its pores, 
differs significantly from the dynamics of the bulk liquid and generally depends on 
the mechanism of interaction between the liquid and the pores of the solid media. It 
will thus be appreciated that a correct interpretation of the measurement parameters 
T.sub.l, t mih . 2 and D can provide valuable information relating to the types of 
f lii-iris involved, the structure of the formation and other well logging parameters of 
interest . 



EtT-tpf Summary Text (7) : 

On the basis of the t snh.2 Rppri-.ra f two specific methods for gas measurements are 
known in the prior art and will be considered briefly next to provide relevant 
background information. The first method is entitled "differential spectrum method" 
(DSM) . The DSM is based on the observation that often light oil and natural gas 
exhibit distinctly separated t. mih.2 measurements in the presence of a magnetic 
field gradient, even though they may have overlapping T.sub.l measurement values . 
Also, it has been observed that brine and water have distinctly different T.sub.l 
measurements, even though their D.sub.O values, may overlap. The DSM makes use of 
these observations and is illustrated in FIG. 1 in a specific example for a 
sandstone reservoir containing brine, light oil and gas. 

Rripf Summary Text (8) : 

A second method known in the art is called "shifted spectrum method" (SSM) . The SSM 
is illustrated in FIGS. 2A-B. Specifically, FIG. 2A shows synthetic T . suh . 2 decay 
curves in a gas bearing zone. The solid curve is for the short interecho time 
( . apprxeq . 0 . 6 msec) and the dashed curve corresponds to a longer interecho time of 
about 2.4 msec. FIG. 2B illustrates the T.mih.2 spectra obtained from the inversion 
of the synthetic echo trains in FIG. 2A. The solid spectrum corresponds to the 
shorter interecho time, while the dashed spectrum line corresponds to the longer 
interecho time. In FIG. 2B the solid spectrum line corresponds to both brine and 
gas. The signal from gas is shifted out of the detectability range, so that the 
single Rppfit-rum peak is due to brine. 

Ftripf Summary Text {9) : 

While prior art DSM and SSM methods provide a possible working approach to detection 
of gas using solely nmt? data, the methods also have serious disadvantages which 
diminish their utility in practical applications. Specifically, typically two 
separate logging passes are required and therefore the methods show relatively poor 
depth matching properties on repeat runs. Furthermore, subtraction of signals from 
different logging passes is done in the T.siih.a spec t. rum domain which may result in 
losing valuable information in the transformation process, especially when the 
received signals have Inw signal-to-noise ratios (SNRs) . In fact, for a typical 
logging pass, low hydrocarbon i nripy (HI) of the gases in the formation, and 
relatively long T.sub.l times generally lead to low SNR of the received signals. 
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After transformation iMi the T.Buh.2 sppctrum domain e^K more information can be 
lost, thus reducing the accuracy of the desired parameter estimates. 



Brief Summary Teyr (10) : 

In the parent application Ser. No. 08/822,567 filed Mar. 19, 1997, which is 
incorporated herein by reference for all purposes, a well logging system and method 
are disclosed for detecting the presence and estimating the quantity of gaseous and 
liquid hydrocarbons in the near wellbore zone. The approach presented in this 
application effectively addresses some of the concerns associated with prior art DSM 
and SSM methods. In particular, the proposed system uses a gradient -based, 
multiple-f requency_UMR logging tool to extract signal components characteristic for 
each type of hydrocarbons. To this end, a data acquisition method is proposed in 
which measurements at different frequencies are interleaved to obtain, in a single 
logging pass, multiple data streams corresponding to different recovery times and/or 
diffusivity for the same spot in the formation. The resultant data streams are 
processed to determine mineralogy- independent water and hydrocarbon saturations and 
porosity estimates. Gas and n-t 1 saturations are used to obtain accurate estimates of 
the water content, permeability and other parameters of interest. 

Brief S ummary Tpyt (11) : 

In most practical applications the approach used in the parent application is 
completely adequate and has been demonstrated to work. However, this approach relies 
on the presence of a static field gradient of typically 10-20 gauss/cm, and may not 
work well in an emerging class of wireline and LWD nmr logging tools that use 1 nwpr 
field strengths, different field patterns and deeper depths of investigations (DOI) . 



R-rief Summary Tpxt {12) : 

T.nwpr field gradients are important in achieving large sensitive volumes and the 
sensitivity required to perform deeper reading measurements. The combination of a 
Inwpr field strength and deeper DOI's results in dramatically reduced static field 
gradients on the order of 0.3-3 gauss/cm. For the measurement method disclosed in 
the parent application, which employs variable wait times (1 sec to 12 sec) and 
fixed echo spacings (typ. 1 msec) , it will generally be difficult to differentiate 
between the nil s water and gas phases in a weak gradient, because the method relies 
on differences in apparent T.mib.2 decay rates due to diffusion in a gradient field. 
One possible solution is to increase the echo spacing (i.e., from 1 msec to 10 
msec) , but such an approach would result in a significant loss of sensitivity 
because the sampling rate is only one-tenth of its original value. Therefore, it 
will be apparent that low-gradient logging tools require modified methods for 
successful gas/oil detection and quantification. - 

Brief Summary Text (14) : 

In accordance with a preferred embodiment of the present invention a novel system 
and method are disclosed for the interpretation of NMR measurements of T.sub.l, D 
and effective t mih 2 parameters made with a NMR logging tool using low field 
gradient magnet icefields. The present invention is based on a multi- frequency, 
gradient based logging tool providing the capability of conducting substantially 
simultaneous^NMR measurements in adjacent non- overlapping resonant volumes of the 
geologic formation of interest. In particular, by hopping the resonant frequency of 
the device and thus sensing non-overlapping volumes of the formation, in accordance 
with a preferred embodiment of the present invention the time between experiments is 
reduced substantially without compromising the T.snh. i r-pl avai-i nnp; or adopting 
imprecise T.sub.l /t snh 2 empirical relationships. 

Brief Su mmary Text (15) : 

In one aspect, the present invention uses the multi -frequency capabilities of an NMR 
tool along with a data acquisition method suitable for the detection of gas on the 
basis of at least two sets of data points corresponding to a long (T.sub.RL) and 
short (T.sub.RS) recovery times, respectively. To this end, an interleaved pulse 
sequence is proposed in which at least two CPMG pulses associated with a first 
resonant frequency are followed by at least two CPMG pulses associated with a second 
resonant frequency. Due to the fact that each resonant frequency excites protons in 
a separate volume of the formation, pairs of complex data points can be collected at 
substantially the same depth mark, such that the first data point corresponds to a 
short recovery time T.sub.RS while the second data point corresponds to a long 
recovery time T.sub.RL. 

Rripf Summary Tevf (16) : 
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In accordance with a prWerred embodiment of the present^Shvention, the sequence of 
data pairs is used next to form two complex time -domain signal vectors x and y 
corresponding to the long and the short recovery times, respectively. Following 
calibration, a difference and a sum signal vectors (x-y) and (x+y) are formed. The 
difference signal is phase corrected to obtain a real-time signal using phase 
information from the sum signal. Next, matched filters corresponding to a gas phase 
and arL_oil phase are computed, in real time, using information about reservoir 
temperature, pressure and other known probe and/or formation properties. The 
phase -corrected difference signal is then filtered using the matched filters to 
separate oil and gas signal components from the input NMR signals. The output 
signals from the matched filters are next used to obtain gas- and oil-porosity 
estimates and further to reconstruct the nil and gas components in the original 
measurement signals. The reconstructed components are subtracted from the sum signal 
1/2 (x+y) to provide only the brine component of the original signal. This brine 
component is finally subjected to T . sub . 2 inversion to obtain an estimate of the 
irreducible water saturation and water porosity. To further increase the separation 
between different hydrocarbon phases in the formation, dif f usion-WFti ghted 
measurements can also be used according to the SSM approach. 

RrHef Summary Text (17) : 

The system and method in accordance with the present invention have been shown to be 
very sensitive due to the fact that the gas and ni 1 components of the original 
signal are determined from the original signals, prior to t snh ? ape<-i-r-nm 
inversion. The method of the present invention has been demonstrated to give more 
accurate parameter estimates than other presently available NMR logging techniques 
and can be used advantageously in low-porosity formations, where low signal-to-noise 
ratios (SNRs) tend to broaden all t mih o components. Additionally, corrections for 
hydrogen index (HI) and incomplete longitudinal recovery are also provided to 
calculate estimates of gas- filled porosity and to correct both apparent NMR porosity 
(MPHI) and free fluid indpy FFI, which are necessary to obtain accurate estimates of 
formation permeability. 

Brief Summary Text (18) : 

In another aspect, the present invention is based on the use of a modified CPMG 
(Carr-Purcell-Meiboom-Gill) sequence that is characterized by two echo spacings. The 
first echo is subject to a variable echo spacing TD r * ranging from approximately 1 ms 
to hundreds of milliseconds. With such a large possible range for the echo spacing 
TD, the diffusion characteristics of the JLliiid.(s) involved can be made a dominant 
factor for amplitude decay. In a preferred embodiment, the second and all following 
echoes are generated with short spacings (TE, typically 1 msec) . In this period, the 
influence of diffusion is negligible in a low field gradient, and a dense sampling 
regime is maintained, resulting in good sensitivity. The amount of amplitude loss 
incurred in the TD interval is directly related to the dif fusivity of the 
hydrogen-bearing fluids. In a preferred embodiment, TD is a parameter that varies 
from a low of TE (i.e., the case in which the modified sequence reduces to a simple 
CPMG sequence) to higha of hundreds of milliseconds. 

TTra wing Dp script inn Text (3) : 

FIG. 1 is an illustration of the differential spftrtnim method for identifying the 
presence of gaseous components. 

rtrawing Dearr-Lpt i on Text (4) : 

FIG. 2A illustrates synthetic_X^sutiL^2 decay curves used in the Shifted Rpe-rt-rum 

method. FIG. 2B annum t snh o spectra obtained from inversion of the synthetic echo 
trains in FIG. 2A. 

Drawing Dearrript i nn Text (7) : 

FIG. 5A shows the hydrogen_Lndex (HI) of methane as a function of depth at 
temperature gradients of 1, 1.5 and 2. degree. F./100 ft. 

Drawing De arri pti np Text (8) : 

FIG. 5B shows the dependency of the longitudinal r-el ayat-inn time T.sub.l as a 

function of depth at temperature gradients of 1, 1.5 and 2. degree. F./100 ft, and 
pressure gradient of 43.3 psi/100 ft. 

Drawing npsrri p^inn Tpxf. (10) : 

FIG. 5D Shows the apparent tranaverae relaxation time T anb . 2 * based On diffusivity 
D.sub.O as in FIG. 5C, diffusion restriction D/D.sub.O, and magnetic field 
temperature gradient of -0 . 18%/ .degree .C. 
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n -rawing Dpfirri pti on Tpyf. (11) : 

FIG. 5E shows values for the T.sub.l parameter of methane gas and light r»i 1 « as a 
function of logging depth. 

Drawing De script, ion Text (12) : 

FIG. 5F illustrates the self diffusion coefficients D.sub.O for methane, water and 
light oil as a function of logging depth. 

Drawin g npsrri pt-inn Tpyt (13) : 

FIG. 5G illustrates the measured T.sub.2R for gas and oi 1 as a function of logging 
depth . 

Drawing Description Tpxh (14) : 

FIG. 6 is an illustration of an interleaved data acquisition pulse sequence for 
T.sub.l weighted and dif f usion^weight£d saturation recovery CPMG echo trains in a 
specific embodiment of the present invention using two tool frequencies. 

Drawing Description Te.xt. (17) : 

FIG. 9A shows sample response of the gas and the oil matched filters in accordance 
with the present invention; FIG. 9B is the magnitude of the Fourier transform of the 
matched filter responses shown in FIG. 9A. 

Drawing npsrripi-inn Text (18) : 

FIG. 10 shows orthogonal i zed filter functions of the gas and the oil matched filters 
in accordance with the present invention. 

Drawing Description Text (20) : 

FIGS. 12A-C illustrate a modified pulse sequence used in accordance with the present 
invention for low field gradient magnetic tools. FIG. 12A illustrates the case in 
which the modified sequence reduces to a CPMG pulse sequence with a single echo 
spacing. FIGS. 12B and 12C illustrate pulse sequences with increased diffusion time 
TD. 

Drawing Descri ption Text (21) : 

FIGS. 13A, 13B 14A and 14B illustrate simulation results showing computed echo 
ampl i i-ndes for the case of no diffusion and signal attenuation due to diffusion. 

Dft tailed Descri ption Text (2) : 

During the course of the description like numbers will be used to identify like 
elements shown -in the figures. Bold face letters represent vectors, while vector 
elements and scalar coefficients are shown in gt-andai-H print. 

D etai led Descri pt i on Text (4 ) : 

In order to separate signal contributions from different f lui ds, an NMR tool must be 

able to operate in a three-dimensional parameter space: T . suh . 2 f transverse decay 
time), measured by a CPMG pulse-echo sequence: T.sub.l (longitudinal polarization 
time) , measured by variable 'saturation- recovery times; and D (apparent, restricted 
dif fusivity) , measured by varying the CPMG pulse-echo spacing .tau. in the presence 
of a magnetic field gradient. 

Detailed Description Text (5) : 

In a preferred embodiment of the present invention these measurements in a moving 
logging tool are enabled using the system illustrated schematically in FIGS. 3 (A-C) . 
In particular, FIG. 3A illustrates, in relatively general form, apparatus for 
carrying out NMR borehole diffusion coefficient determinations in accordance with a 
preferred embodiment of the present invention. The apparatus includes a first 
portion 106, which is arranged to be 1 owe red into a borehole 107 in order to examine 
the nature of materials in the vicinity of the borehole. 

Detailed Desrr-ri pi- i on Text (10) : 

The complex time -domain signal from the RF receiver preamplifier 126 is supplied to 
an RF receiver 40 which optionally receives input from a phase shifter 44. Phase 
shifter 44 receives an input from variable frequency RF source 36. As discussed in 
more detail next, in a preferred embodiment of the present invention phase 
correction is done using signal processing algorithms instead. Receiver 40 outputs 
via an A/D converter with a buffer 46 to computer 50 for providing desired well 
logging output data for further use and analysis. Pulse programmer 146 controls the 
gradient coil power supply 129 enabling and disabling the flow of current, and hence 
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the generation of stat^^or pulsed field gradients, accoTOing to the commands of the 
computer 50. Some or all of the elements described hereinabove as being disposed in 
an above-ground housing, may instead be disposed below ground. Improved devices and 
measurement methods which can be used for the probe 106 are described generally in 
U.S. Pat. Nos. 4,710,713; 4,717,876; 4,717,877; 4,717,878, 5,212,447; 5,280,243; 
5,309,098 and 5,412,320 all of which are commonly owned by the assignee of the 
present invention. A specific embodiment of the tool which can be used in accordance 
with the present invention is also discussed in detail in Chandler et al., "Improved 
Log Quality with a Dual -Frequency Pulsed NMR Tool," paper SPE 28365, presented at 
the 69-th Annual Technical Conference and Exhibition, Society of Petroleum 
Engineers, New Orleans, Sep. 25-28, 1994. The contents of these patents and the 
Chandler et al . paper are hereby expressly incorporated for all purposes. 

n^tailp ri nfsrri pt-inn T<*vt- (11) : 

FIG. 3B is a block diagram of the system in accordance with a specific embodiment of 
the present invention which shows individual block components for controlling data 
collection, processing of the collected data and displaying the measurement results. 
In FIG. 3B the MRX electronics 30 comprises an MRX probe controller and pulse echo 
detection electronics. The output signal from the detection electronics is processed 
by data processor 52 to analyze the relaxation characteristics of the sample. The 
output of the data processor 52 which in accordance with a preferred embodiment of 
the present invention comprises at least two complex time -domain data sets is 
provided to the parameter estimator 54. Measurement cycle controller 55 provides an 
appropriate control signal to the MKT probe. The processed data from the log 
measurement is stored in data storage 56. Data processor 52 is connected to display 
58 which is capable of providing a graphical display of one or more measurement 
parameters, preferably superimposed on display data from data storage 56. The 
components of the system of the present invention shown in FIG. 3B can be 
implemented in hardware or software, or any combination thereof suitable for 
practical purposes. 

Dgtailpri Deacri phinn Tpyf (13) : 

Specifically, in a preferred embodiment the system of the present is capable of 
"hopping" from one operating frequency to another, the effect of which is to shift 
the radial position of the resonant volume for the tool. The frequency shift is 
selected in such manner that at least two non-overlapping resonant volumes are 
formed; each new resonant volume associated with a different frequency being filled 
with fully_relaxfid protons. Hopping between two or more (i.e., K) frequencies thus 
allows reducing the time between experiments approximately by a factor of K, without 
compromising complete T.sub.l measurements or adopting imprecise empirical T.sub.l 
/T. flnh. o. relationships; the logging speed for the tool can accordingly be increased 
approximately K times. This feature is illustrated in FIG. 3C in which hopping 
between two different frequencies is shown to result in conducting measurements in 
two non-overlapping resonant volumes. In the specific example illustrated in FIG. 3C 
each frequency band is about 6 kHz wide and the two mean band frequencies are offset 
by about 15 kHz. This mode of operation forms two concentric annuli, each 0.04 inch 
(0.1 cm) thick, separated center to center by about 0.09 inches (0.23 cm). 

rtetailR H DPfirript-inn TPxt (15) : 

The MRIL tool used in a preferred embodiment of the present invention generally has 
a vertical excitation/response function that can be represented by a near-perfect 
rectangular aperture. In a specific embodiment, a high vertical resolution, 24" 
(60.96 cm) long aperture, or a Inwpr vertical resolution, 43" (109.22 cm) long, 
aperture are used. In order to perform T.sub.l -weighted signal measurements, as 
discussed in detail below, it is required that the formation volume being sensed 
remains substantially unchanged over the course of a recovery period. Specifically, 
for a moving tool, it has been determined that volume changes of about 10-20% still 
provide adequate measurement accuracy. This condition is illustrated in FIG. 4 which 
shows the relationships between saturation- recovery time interval T.sub.R, logging 
speed v and aperture length L. Using the notations in FIG. 4, it is possible to 
impose either a minimum aperture length or a maximum tool logging speed requirement 
which must satisfy the condition for substantial measurement stationarity . For 
example, in the specific embodiment of a 24" long aperture, assuming recovery time 
T.sub.R =2 s, and imposing a 10% accuracy requirement, it can be seen that the 
maximum allowed tool speed is v=5*2.4"/2 s=6 ft/min (3.05 cm/sec). For the alternate 
configuration using a 43" long antenna, under the same assumptions the maximum tool 
speed is about 11 ft/min (5.588 cm/sec) . 

DRtaileH Tteficr-jpti on Text (16) : 
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In accordance with the^resent invention, for the purpos^e of making T.sub.l 
wpi ghhpd measurements with a moving logging tool at least one long 
saturation- recovery (SR) interval is required, preferably of about 8-10 sec. It 
should be noted that for such an interval logging data is substantially insensitive 
to vertical tool displacement because at the end of the interval the formation 
magnetization is already close to an equilibrium. The t, ran aye r.sfr magnetization left 
after a CPMG sequence is quickly dephased in the strong gradient field. At this 
point, a saturation -recovery (SR) measurement can be started, as known in the art. 
The recovered magnetization is read out by the next CPMG train. 

npt-ailp-d ngRr-THpt-inn Tf>xt (18) : 

Random lateral tool motion is a source of concern for the validity of the downhole 
T.sub.l -wpi ghted measurements in accordance with the present invention. The reason 
is that since the sensitive volume for the tool resembles a cylindrical slice of 
about 1-2 mm thickness, lateral swaying of the tool could cause an influx of fully 
polarized magnetization and thus give incorrect measurements. Studies of actual log 
data repeatedly acquired over the same zones, however, show monotonic recovery 
behavior of sequences with increasing SR intervals, indicating that lateral tool 
motion generally has a negligible effect on the actual measurements. Further 
evidence that side effects due to lateral motion of the tool are insignificant is 
provided by the consistency of MPHI and FFI measurements made with the tool, which 
are both independent of the T.sub.l parameter. 

Dptai 1 ed Des^ri phi on Text (19) : 

Another source of concern -i n wmr logging is the relatively shallow 
depth-of -investigation which, due to the generally cylindrical shape of the 
resonance volume of the tool, also depends on the borehole size. Thus, in some cases 
shallow depth-of -investigation along with the fact that invading fluid in the 
borehole replaces gas can lead to a reduction in the gas effect which can be sensed 
by the tool. It should be noted, however, that the MRIL tool's sensitive volume has 
an approximately 4" (10.16 cm) blind zone extending from the borehole wall. The 
presence of such blind zone effectively limits the influence of f 1 iH d invasion. 
Experimentally, in most cases residual hydrocarbon saturations seen by the tool have 
been shown to be sufficient for hydrocarbon detection purposes and can be close to 
uninvaded saturations. 

Detailed Description Text (20) : 

The CPMG pulse sequences used with the MRIL tool in accordance with the present 
invention have been described generally in U.S. Pat. No. 5,212,447 assigned to the 
assignee of the present application. Also discussed in this patent are specific 
methods of -conduct ing_HMR measurements, including derivations of the diffusion 
coefficient D anH/nr t snh 2 The relevant portions of the disclosure of the U.S. 
Pat. No. 5,212,447 are expressly incorporated herein for all purposes. The MRIL tool 
used in accordance with a preferred embodiment of the present invention stores 
multiple pulse sequences downhole, in a memory (not shown in FIGS. 3A and B) within 
probe 106. These sequences are then activated by commands from the measurement time 
controller 55 of the surface system. At the surface, raw tool data are separated 
into data streams and associated with the correct calibration and correction tables 
in data processor 52. An essentially unlimited number of pulse sequences can be used 
quasi-simultaneously, as described in more detail next. In an alternative preferred 
embodiment of the present invention the operation of the tool can be re -programmed 
on command from surface controller 55. 

Detailed Description Text (21) : 

In accordance with another aspect of the invention, LWD nmr logging tools are 
desirable that use 1 ower field strengths, different field patterns and deeper depths 
of investigations (DOI) . An example of such a tool is disclosed in European Patent 
Application 99200567.8 (Document EP 0 940 688 A2) entitled "NMR logging-while 
drilling tool generating an axisymmetric magnetic field having straight contour 
lines in the resonance region." Another example is provided by the basic magnet 
configuration described in U.S. Pat. No. 4,629,986. The above references are 
incorporated herein by reference for all purposes. 

Detailed npsrHptinn Text (23) : 

In accordance with a preferred embodiment of the present invention several 
parameters which correspond to the gas and the rn 1 phases of the formation are 
computed in real time. Due to the fact that logging conditions dynamically change 
during the course of a pass, correction for various factors which may affect the 
accuracy of the measurements have to be made. In the following paragraphs, a brief 



discussion is presentec^un the specifics of the paramete^computations and the 
required corrections used in accordance with a preferred embodiment of the present 
invention to estimate the relative quantities of gas and oil in a formation of 
interest . 

Detailed Description Text (25) : 

It is known in the art that the static field gradient required for downhole 
diffusion measurements induces stimulated echo effects within a CPMG echo train. 
These stimulated echoes partially undergo T.sub.l relaxation and therefore benefit 
less from the re-focusing effects of repeated .pi. pulses in a standard CPMG 
sequence. In accordance with a preferred embodiment of the present invention this 
problem can be treated by introducing the concept of "effective" rel axation times, 
as described in more detail next. 



Detailed DfiBrription Text (26) : 

In particular, it is known that the classic Carr-Purcell expression for spin echo 
attenuation due to transverse rplayatinn and diffusion in a field gradient which, 
using the_atandard notations above, is given by the expression: ##EQU1## 

Ttetailed Description Text (27) : 

strictly speaking is valid only if: (a) the gradient G is small, or (b) if only the 
on- resonance portion of the spin sper-tnim is utilized. As indicated above, however, 
the MRIL tool operates with a relatively strong gradient field, on the order of 
about 15-25 G/cm. In addition, low signal-to-noise considerations make it necessary 
to utilize the full bandwidth of the tool, so that strong of f -resonance effects are 
necessarily included in the echo signals. Thus, for example, even for the simplest 
t snh ? /d experiment which requires at least two different pulse-echo spacings 
.tau. a correction in the expression in Eq. (1) is required in order to avoid 
systematic errors. Consequently, the observed echo decay signal has to be modeled as 

a complex superposition of longitudinal rftlayaUnn, fransvprsp rplayatinn and 

different diffusion times. 



npta-tlpd DRsm ptinn Text (28) : 

Therefore, in accordance with the present invention, the signal observed at the N-th 
echo is modeled as a superposition of all possible combinations of transitions 
between r ran s verse and longitudinal magnetization and is given by the expression: 
##EQU2## 

Detailed Description Text (29) : 

Using the expression in Eq. (2) , the effect of diffusion dephasing is taken into 
account by introducing "effective" i-ransvprsp -relaxation times T. sub. 1 . sup . .dagger . 
and T. sub. 2 . sup. .dagger . given by the following expressions: 

DPtailpd Description Text (30) : 

It can be shown that direct echoes (i=N) decay with a rate 1/T. sub. 2 . sup. .dagger. ; 
indirect echo decay (i<N) is controlled by l/T. sub. 1. sup. .dagger, and by 
1/T. sub. 2 . sup. .dagger . . Without diffusion, indirect echoes decay either slower or at 
the same rate as direct echoes. With very fast diffusion, however, indirect echoes 
drop out faster than direct ones. The effect on combined p-^hn amplitudes primarily 
depends on the receiver's bandwidth and has been determined to require an about 15% 
correction at high diffusion rates. 

T)pfaiTprf Description Text (31) : 

It should be noted that the expressions for the effective relaxation rates in Eq. 
(3) refer to the echo decay process, and not to the recovery of longitudinal 
magnetization, which is controlled by T.sub.l. For gases, both effective relaxation 
times are dominated by the diffusion term in a gradient field and therefore T.sub.l 
>>T. sub. 1. sup. .dagger. .apprxeq.T. sub. 2. sup. .dagger. /3. In this case, the echo train 
decays slightly faster than expected, and an analysis based on the standard 
Carr-Purcell formula will overestimate the diffusion parameter D. This problem is 
corrected by inserting into the Carr-Purcell formula of an effective pulse-echo 
spacing . tau. . sub. ef f , which incorporates the influence of both pulse width and 
receiver bandwidth: 



Detailed Description Text (34) : 

As evident from Eq. (4), the prediction of T. sub. 2 . sup. .dagger . in the gas phase 
generally requires knowledge of the field gradient G, which is dependent on the 
probe temperature. A specific example of measurements of the depth-of -investigation 
(diameter of the sensitive zone) and the magnetic field gradient values., as 
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functions of probe temperature , are summarized in Table 



np^a^^f> d DeHrripti nn Tpxt (35) : 

Typical valnps used in Eq. (4) are . tau . . sub. ef f =0.65 ins, . tau. =26750 s.sup.-l 
G.sup.-l, and G=18 G/cm. Probe temperature, as reported by a sensor embedded in the 
permanent magnet of the MRIL tool, is always recorded, which allows the calculation 
of the field gradient G at any point on the log. 

Detailed Description Text (36) : 

3) Parameterization of HI, T.sub.l and T..suh.2 

Detailed Description Text (37) : 

The matched filter signal processing method of the present invention, described in 
more detail below, requires the calculation of hydrocarbon {ail and gas) signatures. 
These components are assumed to be the non-wetting phase, i.e., to be generally 
characterized by their bulk relaxation properties. As known in the art, the effects 
of temperature and pressure on T.sub.l and D of the gas phase substantially cancel 
each other, resulting in fairly stable and predictable values for both parameters, 
for which mathematical expressions are available. On the other hand, the 
corresponding values for the nil phase are generally dependent on the formation and 
are determined in accordance with the present invention from sample measurements 
conducted prior to the logging experiment. 

Dft-ai 1 eri Desc^r i phi on Text: (38) : 

In particular, the hydrogen_index (HI) of oi T is assumed to be 1.0. The measured 
drop in NME porosity is typically observed in gas zones, because HI. sub. g <1. Most 
natural gases are predominantly methane. FIG. 5A shows HI. sub. g variations between 
about 0.2 and 06 for a methane gas under typical conditions. In overpressured 
reservoirs HI. sub. g can be about 0.7. Accordingly, the gas HI is sufficient to give 
readily detectable signals from gas. For methane gas, in accordance with a preferred 
embodiment of the present invention the corresponding inrtex HI. sub. g is estimated 
mathematically using the expression 2 . 25 . times . . rho . , where .rho. is the gas density 
in g/cm.sup.3, calculated by solving the equation of state. For gases other than 
methane, or for mixed gases, the multiplying factor is less than 2.25. For example, 
for a typical gas mix, characterized as C. sub. 1.1 H. sub. 4. 2, the factor becomes 
2.17. In an alternative preferred embodiment of the present invention, the hydrogen 
index of different hydrocarbons can also be estimated using the expressions 
presented, for example, in "Schlumberger: Log Interpretation 

Principles/Applications," Schlumberger Educational Services, 1989, pp. 5-20 and 
5-21, the content of which is expressly incorporated herein. 

Detailed Description Text (39) : 

A simple power law has been found sufficient to fit published laboratory data for 
longitudinal rplavahinn time T.sub.l of methane gas, as well as log data. The 
expression used in accordance with the present invention is: 

Detailed Description Text (40) : 

where T.sub.l is measured in seconds, the density .rho. in g/cm.sup.3 and the 
absolute temperature T is in degrees Kelvin. Eq. (5) is valid for gas densities up 
to about 0.3 g/cm.sup.3 ; higher densities generally approaching a liquefied gas 
state . 

Detailed Dpsnr i pt i on Text (41) : 

In accordance with the present invention the non-wetting oil phase relaxes with its 
bulk relaxation T.sub.l,o, which is determined, for example, by using viscosity 
measurements of a sample. It has been determined that in order to successfully 
detect liquid hydrocarbons, for the T.sub.l -weight-ed measurements in accordance 
with the present invention a long T.sub.l component (low viscosity) on the order of 
1-2 s is necessary. The relatively large values for the parameter T.sub.l of light 
hydrocarbons provide a mechanism for distinguishing these fluids frnm wat.ar. since 
T.sub.l of water- in rocks is almost always less than about 500 msec. In partially 
hydrocarbon- saturated water wet rock the hydrocarbon -water contrast is even better 
because T.sub.l (and T.sub.2) of water are shorter, due to the fact that water 
typically resides in the smallest pores. 

net-ailed Description Text (42) : 

FIG. 5E shows, on the other hand, values for the T.sub.l parameter of gas (methane) 
and light oils at depths below 4000 feet. Both curves were computed assuming a 
geothermal gradient of 1.5. degree. F./100 feet. The oil curve was computed assuming 
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a temperature -dependen^^iscosity, using the express ion ^^ta. = . eta. .sub. 0 exp (Q/RT) , 
where .eta.. sub. 0 =0.01 cp, Q=10.5 kJ/mol and R=8.314 J/mol/K. As seen in the 
figure, due to different rplavai-ion mechanisms, T.sub.l for methane can be 
relatively short, i.e., between about 2.5 and 4 seconds, while in the specific 
example T.sub.l for oi 1 can be very long (on the order of 10 sec) . standard logging 
practice requires to set the wait time between successive CPMG pulse trains long 
enough for substantially fall recovery (about 95%) of the longitudinal 
magnetization. Accordingly, waiting times for a particular measurement have to be 
adjusted dependent on the specific oil. 

Detaile d Description Text (43) : 

The apparent diffusivity D of a f 1 "id depends both on the self -diffusion coefficient 
D.sub.o and the restrictions imposed by the pore space. In accordance with a 
preferred embodiment of the present invention, an experimental temperature and 
density relationship for unrestricted gas diffusion D.sub.0,g is used, which can be 
expressed mathematically as: 

Detailed Description Text (44) : 

where D.sub.0,g is measured in cm. sup. 2 /s, the temperature T is measured in degrees 
Kelvin and the density .rho. in expressed in g/cm.sup,3. Below 7,000 ft, the 
opposing effects of temperature and pressure stabilize the diffusion parameter 
D.sub.O at a value of about 10. sup. -3 cm. sup. 2 /s. Diffusion restriction in the pore 
space should also be taken into account since the diffusion length (given by 
sqrt (2 . tau.D. sub. 0) ) is approximately equal to 10 .mu.m. D/D.sub.O ratios in rock 
samples at this length scale have been observed ranging from about 0.55 (Indiana 
limestone) to about 0.9 (oomoldic limestone). Sandstone samples have been found to 
cluster in a tight D/D.sub.O ratio range of 0.7-0.8, which is consistent with 
experimental observations of T. sub. 2 , g. sup. .dagger, from log data. 

Ttetai 1 pd Dpscri pti on Text (45) : 

Because of diffusion, the intrinsic r^T axation rate 1/T.sub.2,g for gas is 
negligible compared to l/T. sub. 2 ,g. sup. .dagger, (see Eq. (3)). Similarly, the 
diffusivity of the„o±l phase is small compared to that of the gas phase. 
Consequently, the parameters T.sub.2,o and T. sub. 2, o. sup. .dagger, which are used in 
the matched filter expression considered next are much larger than both 
T. sub. 2, g. sup. .dagger, and also much larger than the total acquisition time required 
to separate oil from gas signals. As indicated above, numerical_v_alii£s. for these 
parameters can be obtained, for example, from sample measurements. 

Detailed Dfifinri pti on Tsxt (46) : 

FIG. 5F illustrates the self diffusion coefficients D.sub.O for methane, water and 
light oil. All curves are based on geothermal gradient of 1.5. degree. F./100 feet, 
and (for gas) a hydrostatic pressure gressure. As seen, the methane D.sub.O is at 
least about 50 times larger than that of water and light oil. The resulting 
contrasts in the measured T.mih.2 (i.e., T.sub.2r) for gas compared to oil are shown 
in FIG. 5G. The plots include the effect of temperature and pressure on T.sub.l (see 
FIG. 5E) and D (see FIG. 5F) for both fluids and the effect of temperature on the 
parameter G for the tool. Moderate restriction diffusion effect on gas and no 
restriction effects on ni T diffusion was assumed, i.e., (D/D.sub. 0) .sub.g =0.7; 
(D/D. sub. 0) . sub.o =1, Comparison of FIGS. 5E and 5G reveals that gas has a high 
ratio T.sub.l /t sub 2 (larger than about 200) at all depths, which is a 
characteristic signature of gas. For light oil, however, this ratio is approximately 
equal to one and rises slowly with increased depth. 

Derailed Description Text (47) : 

Table 2 summarizes expressions for the parameter estimates of different fluids used 
in accordance with a preferred embodiment of the present invention. 

Deta i ] p H Dear ri pt i on Tsxt (48) : 

For illustrative purposes, examples of pre-calculated values for HI. sub.g 
diffusivity and the relaxation time parameters T.sub.l,g and T . sub . 2 , g . sup .. dagger . 
as functions of depth are shown in FIGS. 5A to 5D. In particular, FIG. 5A shows the 
hydrogen index (HI) of methane as a function of depth at different temperature 
gradients; FIG. 5B shows the dependency of the longitudinal relaxation time T.sub.l 
as a function of depth at temperature gradients of 1, 1.5 and 2. degree. F./100 ft, 
and pressure gradient of 43.3 psi/100 ft; FIG. 5C shows the dependency of the 
self -diffusion coefficient D.sub.O of methane as a function of depth; and FIG. 5D 
shows the apparent transverse relaxation time T . sub . 2 . sup . . dagger . based on 
diffusivity Do as in FIG. 5C, diffusion restriction D/D.sub.O, and magnetic field 
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temperature gradient or^O . 18%/ .degree .C. 



TV>t-a-i1gri Description Text (51) : 

As indicated above, the MRIL tool of the present invention is capable of performing 
separate, quasi -simultaneous measurements in different sensitive volumes by simply- 
switching the operating frequency of the tool by a small amount. In accordance with 
a preferred embodiment of the present invention, this multi- frequency capability of 
the MRIL tool is used to provide a new data acquisition method which is particularly 
suitable for the detection of gas on the basis of NMR measurements with different 
recovery times T.sub.Ri. To this end, with reference to FIG. 6, a novel interleaved 
pulse sequence is proposed in which at least two CPMG pulses 210 and 22 0 associated 
with resonant frequency F.sub.2 are followed by at least two CPMG pulses 230, 24 0 
associated with a different resonant frequency F.sub.l. As shown in FIG. 6, the NMR 
measurement is continued next using at least two new pulses 250, 260 at the initial 
resonance frequency F.sub.2, followed by at least two separate pulses 2 70, 280 at 
the F.sub.l frequency. Due to the fact that resonant frequency F.sub.l excites 
protons only in volume 1 of the formation and resonant frequency F.sub.2 excites 
protons only in volume 2 of the formation, pairs 225, 245, 265, etc., of independent 
complex data points can be collected at each depth mark. As shown for illustrative 
purposes in FIG. 6, the first data point in each pair, generally designated as 
B.sub.i corresponds to a short recovery time T.sub.Ri, while the second data point, 
generally designated as A.sub.i, corresponds to a long recovery time T.sub.R2. 

Detaile d rteFtrrr-ipt-.irm T^xt (52) : 

Thus, using the data acquisition sequence illustrated in FIG. 6, by "hopping" the 
resonance frequency F.sub.i of the tool, and alternating between adjacent resonant 
volumes of the formation one can obtain a sequence of signal pairs, each pair 
corresponding to substantially the same depth mark in the formation, but measured at 
different recovery times. It should further be noted that data for two different 
recovery times need not necessarily be obtained from only two different frequencies. 
For example, two or more measurements associated with different frequencies can be 
combined (i.e., avpragp.d l to result in a single data stream corresponding to either 
a short, or a long recovery time. Furthermore, it should be clear that by using more 
than two resonance frequencies, and applying a correspondingly larger number of 
pulses in each resonant volume, the data acquisition method of the present invention 
can easily be extended to the more general case of M-tuple measurement data sets, 
each measurement point corresponding to a different recovery time T.sub.Ri. 

nRi-a-ilgd rescript-. ion Text (54) : 

The data acquisition method was described above with reference to identical CPMG 
sequences, which mode is referred to in -the present application as T.sub.l -wfti ghtfid- 
acquisition. Data from this acquisition mode is suitable for the Differential 
Rpprtmm Method (DSM) described in the Akkurt et al. paper. Notably, however, the 
method is also suitable for direct signal subtraction in the time domain, as 
described in more detail next. 

DRt-a-i 1f>d Ttefirrr-ipt-ion Text (55) : 

In an alternative preferred embodiment of the present invention, a novel data 
acquisition mode referred to as T.sub.l - and dif f usion^stfeighted acquisition can 
also be used. As indicated above with reference to the SSM method, the contrast 
between liquid and gas signals can be enhanced by using a slightly larger pulse-echo 
spacing for the CPMG train associated with the shorter recovery interval. This 
embodiment is illustrated in FIG. 6 using two different intervals .tau..sub.i for 
each successive pulse in the same resonance volume. It has been found, however, that 
it is not necessary to eliminate the gas signal altogether. For example, an increase 
by only 40% in the pulse echo .tau. has been found to cause a 50% decrease in the 
diffusion- induced part of T . sub . 2 . As indicated above with reference to the SSM 
method, because of diffusion dominance, the effect is much more pronounced for gases 
than for liquids, and can accordingly be used to enhance the separation of the two 
phases . 

Dpta-jlgd Description Text (56) : 

In another aspect of the invention, a modified pulse sequence is disclosed 
preferably for use with low field gradient tool measurements. The key idea is to use 
a modified CPMG (Carr-Purcell-Meiboom-Gill) sequence that is characterized by two 
echo spacings. The first echo is subject to a variable echo spacing TD; ranging from 
approximately 1 ms to hundreds of milliseconds. With such a large range possible for 
TD, the diffusion characteristics of the fluid (s) involved can be made a dominant 
factor for ampl -i tnrig decay. The second and all following echoes are generated with 
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short spacings (TE, typically 1 msec) . In this period, the influence of diffusion is 
negligible in a low field gradient, and a dense sampling regime is maintained, 
resulting in good sensitivity. The amount of amplitude loss incurred in the TD 
interval is directly related to the diffusivity of the hydrogen -bearing fluids. TD 
is a parameter that varies from a 1 nw of TE (i.e. a simple CPMG sequence) to highs 
of hundreds of milliseconds. 



Detailed Description Text (57) : 

An example of the sequence used in a preferred embodiment is shown in FIGS. 12A-C. 
In particular, FIG. 12A illustrates the trivial case of TD=TE, which is a standard 
CPMG sequence used as a reference. The TE interval is chosen short enough, such that 
the influence of diffusion on the echo amplitude decay is negligible. As shown in 
FIGS. 12B and 12C, in accordance with this aspect of the invention the TD interval 
is increased in successive measurements. In a preferred embodiment, the granularity 
of the TD interval is chosen such that the later echoes line up in time with those 
of the CPMG scan and a one-to-one correspondence can be established. In a preferred 
embodiment, the second TD valnp can be chosen to attenuate or eliminate signals from 
the gas phase only. The difference between corresponding echoes (as shown in FIGS. 
12A and 12B) would therefore be signal corresponding to the presence of gas. In 
accordance with the preferred embodiment the third TD valnp can be selected to 
attenuate or eliminate both gas and free water signals. Accordingly, the remaining 
signal (illustrated in FIG. 12C) could come only from the oil phase, which typically 
has the slowest self -diffusion rate. 

Detailed De script inn Teyt (60) : 

where A is the echo ampl i tude, g is the gyromagnetic ratio, G is the field gradient, 
Tx is the interval TD or TE, respectively, and D is the self -diffusion constant. For 
the purpose of this discussion, spin- lattice and spin- spin relaxation are ignored 
for the moment. To illustrate the method, the following values can be used: 

Detailed Description Text (61) : 

It will be appreciated that other values can be selected in different applications. 
With the above assumptions, in the TE interval and in the case TD=TE, even for the 
fastest diffusion case, the diffusion- induced signal decay rate is only: 

npfailprf Dpflpriptinn Text (62) : 

which is small compared to typical intrinsic t. snh- ? rates (.about. 0.25 sec. sup. -1) 
and typical pulse-echo train lengths (. about. 1 sec). Therefore, in accordance with 
this aspect of the invention during the TE interval the influence of diffusion can 
be ignored. Further, it will be appreciated that the case TD=TE can serve as a 
reference without echo attenuation due to diffusion. 

Detailed Description Text (63) : 

In accordance with the present invention, the next TD valnp should attenuate the gas 
signal by approximately 86% (i.e., two time constants). It follows that TD should 
obey the following equation: 

Detail eri Description Text (65) : 

The ainplitiide of the free^Hafcer signal is attenuated by a factor of 
Detailed Description Text, (66): 

The echo train with TD=32 ms contains 14% of the gas signal, 82% of the free -water 
signal and 98% of the oi l signal (from an analogous calculation) . An echo-by-echo 
subtraction from the TD=1 ms train isolates the gas signal as follows: 86% gas 
signal, 18% free -water signal and 2% oil signal. The bound -water signal is 
eliminated by considering only echoes that occur after 32 ms. 

np<-ai1pH Description Text (67) : 

For the current example, a third TD valnp is chosen equal to 70 ms to isolate the 
oi 1 signal. In this case, the gas signal is completely eliminated, the free -water 
signal is attenuated by 86% and the oil signal remains at 82% of its value. 

Detailed Description Text (69) : 

It will be appreciated that the above valnpR are used as guidelines only and can be 
replaced in practical application using the described computation techniques. 

Detailed Description Text (71) : 

In a first aspect, the method for fluid detection in accordance with a preferred 
embodiment of the present invention assumes data acquisition in the presence of a 
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static magnetic field gradient in the range 10-30 G/cm. The method further requires 
at least two separate measurements with different saturation- recovery times T.sub.Ri 
and/or different echo-to-echo spacings, and is implemented using the data 
acquisition sequence illustrated in FIG. 6. In addition, the very low 
signal-to-noise (SNR) levels which are due to HI losses and incomplete magnetization 
recovery in a typical nmr measurement necessitate signal detection using two-channel 
complex data stream. Therefore, in a preferred embodiment of the present invention 
data is acquired in two orthogonal channels, and avpragpri over a vertical logging 
interval to increase the SNR for the acquired measurement data. 

Detailed Descri ption Tsxt. (72) : 

Turning next to FIG. 7, it shows in semi -block diagram form the signal processing 
method in accordance with a preferred embodiment of the present invention. 
Specifically, the determination of water, nil and gas saturations in the sensitive 
volume begins by performing at least two interleaved T.sub.l -weighted measurements 
to separate the wetting phase {brine, surface -dominated relaxation) from the 
non-wetting phase (light hydrocarbons, bulk-dominated relaxation) . Optionally, these 
measurements can be dif f usion-wei ghted as well. As shown in FIG. 7, this results in 
two parallel data sets of complex time-domain data. Data sequence 310 generally 
corresponds to data obtained from the long recovery time T.sub.R2, while data 
sequence 32 0 corresponds to data obtained from the short recovery time T.sub.Ri. 
Between about 150 and 300 data points are used in each sequence. Preferably, the 
recovery times used are about 2 sec for T.sub.Ri and about 8 sec for T.sub.R2. Pairs 
of echo trains are formed by matching overlapping short and long TR intervals 
thereby minimizing the systematic variations introduced when formation bed 
boundaries are crossed. 



D<=>1-ai1ed Description Text (73) : 

Following the data acquisition process, in block 325 the two data sets are corrected 
using calibration data from blocks 315. Specifically, such calibration data is 
obtained from samples at room temperature, with 100% water saturation, for example 
in a wafpr-fi 11 ed tank. FIG. 8A shows two such calibrated data sequences as 
functions of time. As shown in FIG. 8A, the magnitude values for the measurement 
signals are conveniently calibrated in porosity units (p.u.). Skipping for a moment 
block 330, next the complex difference between the signals in each data pair is 
obtained in subtractor 335 to eliminate the brine contribution; the sum signal is 
obtained in adder 340 to estimate the input signal phase and phase-correct the 
difference signal in block 335 accordingly. Specifically, it has been observed that 
while the absolute phase of the nmr signal is subject to slow variations due to hole 
and tool conditions, it shows excellent short-term stability. Therefore, the phases 
of the sum and the difference signals- are approximately equal. In accordance with - 
the present invention this feature is used to correct the phase of the difference 
signal on the basis of a phase estimate for the sum signal which was found to be 
comparatively more accurate. 

Detailed Descr iption Tp.xt (74) : 

In particular, a depth^axerag^d signal phase is computed in block 350 from the 
complex sum signal. If proper frequency and phase cycling has been employed during 
data acquisition, all sum and difference echoes have the same avpragp. phase. The 
phase estimated in block 350 is used to rotate, in phase correction block 360, the 
phase of all data points in the complex difference signal into the real axis. Such 
rotation yields the true absorption mode (real-valued) signal component. The 
dispersion signal component (imaginary valued) can be discarded. 

Detailed Description Text (75): 

Based on the parametric representations for relaxation times and diffusion 
characteristics of the non-wetting hydrocarbon phases which are computed, as 
indicated in the section "Signal Modeling and Corrections", and stored in block 365 
for the oil and block 370 for the gas components. Matched filters representing the 
liquid and the gaseous phases are computed next in blocks 375, 380 in the echo- time 
domain, using the expressions: 

D^l-ailed Description Text (78) : 

In general, the filter functions in Eqs. (7a-7b) are not orthogonal and cannot be 
directly applied to the data. (See FIG. 9A) Instead, the ampl iturift responses to 
these filters are extracted from the phase -corrected difference signal d(t) by 
solving, in block 390, the overdetermined equation system, 



D^t-ail^d Description Text (79) : 
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in a least-squares sensed The solution is found by solving the following 2. times. 2 
equation system for the ampl -i t-ude vector x: 



n^t-ailpfl nPRrripfinn Tpyf (80) : 

If the difference signal d(t) was properly scaled in p.u., the first element of the 
solution vector x is_oil-f illed porosity .PHI.. sub. o, obtained in block 392, and the 
second is gas-filled porosity .PHI.. sub. g (block 394). Calculating backwards from 
these numbers, the properly scaled oil and gas signatures can be reconstructed in 
block 395, and subtracted from the complex sum signal in block 342. The remainder is 
the signal originating from brine, which, as wetting phase, is sensitive to the 
surface-to-volume ratios in the remaining pore space available to water. 

De-tailed Description Text (81) : 

In block 4 00, a t . aiib _ o. inversion mapping is constructed, as discussed, for example 
in the Prammer et al. paper above. The results are used, in block 410, to estimate 
the T.nuh.2 spprhnim of the signal and in block 420 to estimate the water -bound 
porosity. 

Dpt-ail^d npsrHptinn Text (82) : 

These ratios are indicative of the non- producible water volume held in place by 
capillary forces (BVI) , which is computed in block 440. On the other hand, the total 
area under t-h<=» t gnh ? curve is interpreted as watej:- filled porosity .PHI.. sub. w 
which is computed in block 420. The total nmp porosity can be computed in block 430 
using the expression: 

pRt-.ailed nPsrHptinn Text: (83) : 

The f r<qp> -f luid index as seen from the wat.pr phase is augmented by ail and gas 
porosity: 

Detailed Description Text (84) : 

From . PHI . . sub.NMR, BVI and FFI, a permeability estimate can be calculated in block 
450, which depends only on HMR-derived quantities. 

Detailed Description Text (85) : 
Processing of Dif fusionj^weighted Data 

Detailed Description Text (86) : 

Turning back to block 330 in FIG. 7, for a dif f usion-wei ghted measurement a 
complication arises from the different sampling grids employed in acquiring the data 
sets which make up a data pair. In this case, data from the shorter echo spacing is 
mapped onto the wider sampling grid by linear interpolation between complex echoes. 
Pi f fiiRio n-wpighting is taken into account to give the following matched-filters 
expressions : 

ppt-ailgd DRKcripl-ion Text (87) : 

The dif f usion^iHeighted data is next processed following the flow graph in FIG. 7. 
The combined t sub _ i -wpi ght-pri and dif f usionjjtfeighted measurement is advantageously 
used in cases where the gas filled porosity and HI are relatively low, and 
correspondingly the SNR of the measurement is relatively low. 

Dgt-ailgri Dpficri pfinn Tpyf (88) : 

Error estimates for t _ Hub . i -weighted and dif fusionzJateigliLed data acquisitions in 
accordance with a preferred embodiment of the present invention can be obtained 
using the following considerations. The input data consists of two data sets, 
weighed by different recovery times and possibly sampled with different echo 
spacings. Each set is individually calibrated for HI=1.0. To compute the uncertainty 
of the parameter estimates, it is assumed that the noise in each data set is random 
and has Gaussian distribution with standard deviation .sigma.=2. As indicated in 
blocks 350 and 360 above, a depth -averaged signal phase is computed from the complex 
second echo in the sum. The estimated phase is used to rotate all complex 
differences into the real axis. Allowing for a small error in phase estimation, the 
noise component in the real -valued difference signal (dt) is approximately 1.5 p.u. 
Gas and oil porosities are given as least-squares solutions in block 390 in FIG. 7. 
Formally, the least-squares solution can be written as: 

npt-ail ed Descr iption Tpxt (89) : 

The sensitivity of the solution to random errors in the input is given by the 
condition number of the square, positive-definite matrix A.sup.T A. The 
orthogonal i zed oil-sensitive filter function f(t)', is the first row of the 
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expression (A.sup.T A)^?up.-1 A.sup.T. The second row o^^his expression contains 
the orthogonal filter f(t)'.sub.g, which is sensitive to the gas component. The 
orthogonal filter responses are shown in FIG. 10. Filter response functions are 
computed as follows: ##EQU3## 

npl-ailp d Description Text (90) : 

The average output uncertainty was determined by Monte Carlo simulation. Using 
100,000 samples, and assuming parameters, HI. sub. g =0.5, T.sub.l,g =5 s, the 
standard output deviation in the answer for gas filled porosity is .apprxeq.2.5 p.u. 
The uncertainty in oil- filled porosity is substantially reduced and is approximately 
equal to 1 p.u., dependent on T . sub . 1 , o . 

Dp-l-ailpd Description Text (91) : 

As noted, in another aspect the invention uses low. field gradient tools resulting in 
reduced static field gradients on the order of 0.3-3 gauss/cm. The modified 
processing sequence illustrated in FIGS. 12A-C is used in such case in a preferred 
embodiment . 

Derailed Description Text (93) : 

Nomenclature A=design matrix for least-squares problem A.sup.T ^transpose matrix of 
A BVI=bulk volume irreducible water, p.u. D=restricted diffusivity, cm. sup. 2 /S 
D.sub.o =unrestricted diffusivity, cm. sup. 2 /S d( t) =dif ference function f (t)=filter 
function FFI=free fluid index, p.u. G=magnetic field gradient, G/cm HI=hydrogen 
indpy, relative to wafpr MPHI=apparent nmr porosity, p.u. T=absolute temperature, 
.degree. K. T.sub.l =longitudinal rplavafinn time, s T. sub. 1 . sup. .dagger . =pseudo 
hransvprgp relaxation time, S T snh 0 =frangyprfiP rel aval- ion time, S 
T. sub. 2. sup. .dagger. =apparent t ran sversft re 1 a xa f. i on time, s TE=CPMG echo-to-echo 
delay (TE=2.tau.), s TR=recovery time, s x=solution vector to least-squares problem 
.PHI.=porosity, p.u. . PHI . . sub.NMR =corrected_HMR porosity, p.u. .PHI.. sub. n =CPMG 
phase, n=l or 2 . gamma. =gyromagnetic ratio, rad.sup.-l G.sup.-l {character 
pullout}=density, g/cm. sup. 3 . sigma. = standard deviation .tau.=CPMG pulse-to-echo 
delay ( . tau. =TE/2 ) , s . tau . . sub . ef f =dif fusion-effective CPMG delay, s 

Derailed Description Text (94) : 
Subscripts g=gas o=oil 

Derailed Description Text (97) : 

FIG. 9A shows the response of the gas (block 375) and the oil (block 380) matched 
filters in accordance with the present invention; FIG. 9B illustrates the magnitude 
of the Fourier transform of the matched filters responses shown in FIG. 9A. 

Detailed Description Text. (98) : 

FIG. 10 shows the orthogonalized filter functions of the gas and the oil matched 
filters in accordance with the present invention. 

Detailed Description Text (100) : 

Logging tests were performed in single -frequency and dual frequency operating modes 
to assess the amount of vertical and lateral motion, which could affect the accuracy 
of the t ctnh i -wei ghted measurements conducted in accordance with the present 
invention. A similar operating procedure was already in place to set the optimum 
delay time between measurements in the pay zone. Data acquired under a variety of 
hole conditions and T.sub.l values, were examined; logging speeds in this particular 
mode were typically 300 ft/hr. In all cases, increasing the recovery time interval 
resulted either in a monotonic increase in nmr amplitude or in no increase. Sudden 
increases in ampl i rude at short recovery intervals, indicative of uncontrolled tool 
motion, were not observed. 

Derailed Description Text (101) : 

The effect of invasion was studied by comparing results from wells drilled with 
water -base muds (WBM) ancL_Q±l-muds (OBM) . The WBM-drilled formations generally 
suffer high invasion and residual <">i T and gas saturations are low. Nevertheless, in 
many cases gas quantities above the detectability threshold are present, possibly 
due to backsaturation of gas into the invaded zone. As indicated above, another 
factor aiding the MRIL.RTM. is the 4" blind zone into the formation. OBM filtrate 
generally invades less and is therefore better suited for near-borehole saturation 
measurements, oil filtrate mixes with the connate oil and replaces it to a certain 
extent. Because the filtrate has low viscosity, OBM aids the described hydrocarbon 
detection method by supplying a slowly relaxing component with known T.sub.l. We 
recommend performing T.sub.l anrf T.mih.2 measurements of filtrates in the laboratory 
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at 1 MHz to assess the^ffect of OBM invasion. 



net-ailed npsrnpUnn Text (103) : 

Log results are shown in FIG. 11. For the example illustrated in FIG. 11 nmr data 
for recovery times of TR.sub.l =3 s and TR.sub.2 =8 s were acquired in separate 
passes with a pulse-echo spacing . tau.=0.6 ms. Both apparent nmr porosities are too 
1 nw in the gas zone (shown in track 4) . Complex echo sums and differences from these 
echo sets were computed. The sum had a constant phase of -2.1 rad, which was used to 
convert the difference signal to real -type values . Matched filters for the oil 
component : 

Tie.t.a i led Description Text. (105) : 

were computed and applied to the data in a least-squares sense as described above. 
The oil-filter response was essentially zero (not shown) , the gas-filter response is 
plotted in track 1. 

DPit.flilPid npsdription Text. (106) : 

Gas -corrected HMR porosity indicates that mud filtrate did not invade the sampling 
diameter (15" at a probe temperature of 100 .degree . C. ) , or that gas did backsaturate 
into the invaded zone. In either case, the results clearly indicate the val hp of 
hydrocarbon saturation measurements near the borehole wall. 

Detailed nesr-Hpt-inn Tpyf (107) : 

Simulation Results for Embodiments Using t.ow Field Gradient 
Detailed Descri ption Text (108) : 

The NMR response can be computed by brute force from the Bloch Equations, or, more 
efficiently, by a Riemannian spherical affine projection onto the complex plane. The 
latter procedure ignores the effects of relaxation, which is appropriate here. Some 
results are shown in FIG. 13 and FIG. 14. 

Detailed DftRcripr.ion Text (109) : 

FIG. 13 treats the case of little or no diffusion. Switching TD to a 3falue larger 
than TE results in loss of signal in echo 2 and onward. This effect can be 
understood by considering that echoes 2, 3, etc. contain both direct spin echoes, as 
well as indirect ("stimulated") echoes (See, for example, Goelman and Prammer: "The 
CPMG Pulse Sequence in Strong Magnetic Field Gradients With Applications to Oil-Weil 
Logging," J. Magn. Reson. A 113 (March 1995), p. 11.) The significance of this 
observation is that in general calibrations must be performed separately for the 
case when TD=TE and TD>TE. Also, in order to recover signal-to-noise, preferably 
twice as many experiments per TD>TE scalue should be -performed compared with the case 
when TD=TE. 

Detaile d npsrHption Tpxt (110) : 

The simulation results illustrated in the figure show computed echo amplitudes for 
the case of no diffusion. Only the first 50 echoes are shown, more typically 
100-1000 echoes would be used in accordance with a preferred embodiment. The left 
diagram in FIG. 13 illustrates echo amplitudes for the case TD=TE=1 ms. The right 
diagram illustrates the case when TD=32 ms, TE=1 ms. Note that the first echo has 
the same ampl i 1-ntie as in the left figure, but all other echoes are attenuated. This 
is a result of the loss of stimulated echoes that are not carried over to the second 
echo and does not reflect attenuation due to diffusion. The amount of attenuation 
due to stimulated echo effects is predictable and can be calibrated in a preferred 
embodiment on a test sample. 

Detailed npffrnp^inn Text (111) : 

FIG. 14 shows the strong attenuation effect if fast diffusion is present. A (weak) 
gradient of 1 gauss/cm and a diffusion coefficient of 10. sup. -3 cm. sup. 2 /sec were 
assumed for these simulations. The simulation results show computed echo ampl i tndes 
for the case of rapid diffusion (D=10.sup.-3 cm. sup. 2 /s) , typical of vapor phase 
under reservoir conditions) . The left diagram illustrates echo ampl i tndfis for the 
case TD=TE=1 ms. The right diagram illustrates the case when TD=32 ms, TE=1 ms. 
After accounting for the loss of stimulated echoes, the gas signal has been 
attenuated to 14% of its original value.. The gas signal can be recovered in 
accordance with a preferred embodiment by individually fitting each echo train and 
comparing total amplitudes, or by subtracting the ampl i tndfts of corresponding 
echoes. Thus, by subtracting echo-by-echo the right diagram from the left diagram 
(ignoring all echoes occurring at times less than TD) , the gas phase signal can be 
isolated. 
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Derailed nesrr-jpti nn Text (112) : 

Appendix A provides listings of several NME Simulation Programs in the Matlab 
programming language, which can be used in accordance with the present invention to 
perform the above -referred computations. 

Detail pd ngRr-r-iption Tpxt (114) : 

Job Planning. In accordance with the present invention job planning consists of 
estimating dif fusivity_}calii£s for gas, ail and free water and selecting appropriate 
TD and TE values based on forward -mode ling of signal components to achieve maximum 
contrast between f lu-id phases. The simple calculation shown above is an example. 
Following is a brief description of a model-based linear inversion algorithm that 
can be used in accordance with a preferred embodiment for the forward modeling of 
signal components. 

Detail fid Description Text (116) : 

Let a,b,c, . . . denote the indices for the free fluid components, e.g., gas, water 
and oil . Let also S.sub.i (i=a,b,c . . . ) denote a signal component and let V.sub.i 
(i=a,b,c . . . ) be the volume occupied by a fluid component. S and V are related as 
follows : 

nptailpd nftsrripfinn Tf^vt (117) : 

where K is a calibration constant containing the system gain, the gyromagnetic 
ratio, the operating frequency, etc.; T is the absolute temperature; HI. sub. i is the 
hydrogen index (i.e., hydrogen density relative to watex) of the i-th signal 
component; P.sub.i is a component -specific polarization term, typically of the form 
(1-exp ( -Tw/Tl . sub. i) ) , where Tw is the pre-measurement wait time (polarization time) 
and Tl.sub.i is the Tl time constant of the i-th component. 

Detailed Description Text (118) : 

Measurements are made with at least three different val hrr for TD. In a specific 
example: TD=1, 32 and 70 ms . 

Detailed Description Text. (119) : 

The echo spacing qTE is held constant (example: TE=1 ms) . All echoes occurring before 
70 ms are ignored to eliminate the influence of fast -decaying components from water 
or oil phases. A relatively short acquisition window is selected (example: 10-100 
ms) . Within this time window, the inherent 12. decay is small for all free-fluid 
components and all echoes within the time window can be summed and averaged. This 
step improves the SNR by a factor of 3-10. The net result are three numbers, 
representing the avfiragpd observed signal* amplitude after 70 ms under the conditions • 
of TD=1, 32 and 70 ms. Let A.sub.j, (j=l,2,3, . . . ) be these numbers. The signals 
S.sub.i and the measured amp! -i tndea A.sub.j are related to each other as follows: 

Detailed Descr i pt- i on Tpyt (121) : 

The matrix elements M.sub.ji encode the response of each fluid component to a 
specific value, of TD. There are given by 

Detailed npRrHpl-inn Text (125) : 

And, finally, the f lni d volumes themselves: 

Detailed TV>firHpfinTi Tpvt (127) : 

Postprocessing. In accordance with a preferred embodiment, in postprocessing, 
calibrations are applied to eliminate stimulated echo effects that distort echo 
amplitudes . Next, the echo trains for various TD values are processed either 
simultaneously or in pairs. If processed in pairs, only the common time windows 
(from the longest TD interval onward) are considered. The subtraction of 
corresponding echoes enhances a particular fluid phase: the gas phase between the 
two shortest TD values; the water phase between the next TD pair, and the ail phase 
from the longest TD value, only without subtraction. Lastly, correction values are 
applied that account for (a) incomplete polarization, (b) incomplete suppression by 
diffusion, and (c) hydrogen density. 

Detailed Description Text (128) : 

The law- field gradient processing method in accordance with the present invention 
can be further combined with the frequency hopping method using variable wait times. 
FIG. 15 is an illustration of a multi- frequency TD/TE pulse sequence. Each box in 
FIG. 15 represents a TD/TE sequence. In a preferred embodiment, the TD parameter is 
systematically cycled through all of its values . The wait time Tw is kept constant 
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and is typically in the range of 8 to 12 seconds, although other values may be used 
in specific embodiments. FIG. 15 A shows an implementation of the method of the 
present invention using a modified TD/TE pulse sequence for two frequencies (i.e., 
two sensitive volumes) . FIG. 15B illustrates the use of the modified pulse sequence 
in the case of four operating frequencies (sensitive volumes) . Signal components 
with long T.sub.l rplaya^ion values can be selectively attenuated by choosing wait 
times between CPMG and modified CPMG trains that are short compared to T.sub.l. 

rtetfli led Description Paragraph Equation (22) : 

M.sub.ji =exp (- 1/12 ( .gamma. .sup. 2 G.sup.2 TD. sub . j . sup . 3 D.sub.i). .gamma, is the 
gyromagnetic ratio (2 .pi. 4258 sec. sup. -1 /gauss), G is the field gradient, and 
D.sub.i is the diffusivity of the i-th component. Continuing with the example: G=l 
gauss/cm, and D.sub.i =10. sup. -3 cm. sup. 2 /s for the vapor (gas) phase, D.sub.2 
=10. sup. -4 cm. sup. 2 /s for the free^Hatex phase, and D.sub.3 =10. sup. -5 cm. sup. 2 /s 
for the oil phase. 

Derailed nPsrHptinn Paragraph Table (2) : 

TABLE 2 Spin Lattice Relaxation Self Diffusion coeff . Eluid (sec) (cm. sup. 2 /sec) 

parameters/units Methane T.sub.l,g = 2 5. 10. sup. 3 . rho. /T. sup . 1 . 17 D.sub.Og = 

8 . 5 . 10 . sup. -7 T.sup.0.9 /.rho. Gas density .rho. in gramm/cc; T-abs temp in Kelvin 

Qi 1 T.sub.l,o = 1.2 (T/298) / .eta. D.sub.0,o = 1 . 3 (T/298 ) / . eta . .eta. is oil viscosity 

in cp; T-abs temp in Kelvin Water T.sub.l,w = 3 (T/2 98) / .eta. D.sub.0,w = 

1.2 (T/298) / .eta. Same as above 

ripl-ailpd Description Paragraph Table (3) : 

D = 10. sup. -3 cm. sup. 2 /s for the vapor (gas) phase, D = 10. sup. -4 cm. sup. 2 /s for 
the frp p-wat-pr phase, and D = 10. sup. -5 cm. sup. 2 /s for the oil phase. TE = 0.001 
sec, G = 1 gauss/cm, g = 2p4258 sec. sup. -1 /gauss (proton gyromagnetic ratio). 

netailed npflrr-ipfinn Paragraph Table (4) : 

TD value results in signals from. . . 1 ms (=TE) all fluid phases, wei ghted only 
hydrogen density and incomplete polarization (in case of short wait times) , 32 ms 
free water, nil and some bound water (attenuated by intrinsic T . auh . 2) , 70 ms oil 
only (attenuated by intrinsic T_ anh 2 ) . 

net-ailed Description Paragraph Table (5) : 

TABLE 3 parameters for example data, gas temperature: 100. degree. C. gas pressure: 
9000 psi gas type: CH.sub.4 CH.sub.4 density: .rho. = 0.26 g/cm.sup.3 CH.sub.4 
hydrogen index: HI. sub. g - 0.6 longitudinal relaxation T.sub.l = 6s time: 
unrestricted diffusivity: D.sub.O = 0.7 .times. 10. sup. -3 cm. sup. 2 /s diffusion 
restriction (est.): D/D.sub.O = 0.8 magnetic field gradient G = 18 G/cm effective 
pulse-echo . tau. . sub.ef f = 0.65 ms spacing: apparent transverse T.sub.2. sup. .dagger. 
= 60 ms relaxation: T.sub.l /T . sub . 2 . sup .. dagger . contrast: .about. 100 

Derailed nesrri prion Paragraph Table (6) : 

APPENDIX A % CPMG.M - - - CPMG sequence % % NMR Simulation Program. % % . COPYRGT . 
Copyright, M. G. Prammer, NUMAR Corp., 1991, 1992, 1993. % Modification History: % 
Oct. 1991, MGP, Written. % Jul. 1992, MGP, modified pulsim. % 

%%%%%%%%%%%%%%%%%%%%%%% experimental parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%% echo 
on; 12 = le6; % 12. sec Assumed true 12. TE = le-3; % TE sec Echo time NE = 50; % NE 
Number of echoes t90 = 100e-6; % t90 sec 90 -degree pulse width tl80 = 160e-6; % tl80 
sec 180 -degree pulse width AT = TE-tl80; % AT sec Acquisition time dt = 20e-6; % dt 
sec Dwell time A = l/{2*tl80); % A Hz Pulse amplitude fmin = -5000; % fmin Hz 
Leftmost freq sample fmax = 5000; % fmax Hz Rightmost freq sample df =40; % df Hz 
Frequency resolution f jitter = 10; % f jitter Hz Random frequency jitter figure (1); 
echo off; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% t=0:dt:AT; % create acquisition time axis 
f=fmin:df :fmax; % create frequency axis [dum,NF] =size (f ) ; % NF . . . no of frequency 
points sx=[]; % real channel signal sy=[]; % imaginary channel signal % establish 
equilibrium mx=zeros (size (f) ) ,- my=zeros {size (f ) ) ; mz=ones (size (f ) ) ; % 90 pulse; 
phase = 90 [mx, my , mz] =pulsim ( j *A, t90 , f , mx,my , mz) ; 

[mx,my,mz]=evolve(T2,t90+(TE-t90-tl80)/2,f ,mx,my,mz) ; %!!! for nr=l:NE, 
%fprintf ( "Echo %.0f . backslash . n" , nr*l . 0 ) ; % 180 pulse; phase = 0 
[mx,my,mz] =pulsim(A, tl80, f ,mx,my,mz) ; [tmpx, tmpy] =observe (dt , AT, f ,mx,my,mz) ; % 
compute response [mx,my,mz] =evolve(T2,TE-tl80, f ,mx,my,mz) ; % evolve 
[mxl,myl,mzl] =evolve (I2^(TE-tl80) /2, f ,mx,my,mz) ; % echo profile sx=[sx ; tmpx]; 
sy=[sy ; tmpy]; s=sqrt(sx. 2+sy. 2); plot(t,sx); grid; xlabel ("Time in sec after pi 
pulse."); drawnow; % introduce frequency jitter to mask numerical discretization f = 
f +f jitter*randn(size (f ) ) ; end % create entire timeline tall=[],- for nr=l:NE, 
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r^fct :AT/2) ] ; end % plot entire trf^ 



tall=[tall nr*TE+ (-AT/^dt :AT/2) ] ; end % plot entire traTri tmp=sx' ; sxall=tmp ( : ) ; 
figure{2); plot (tall , sxall) ; axis([0 0.05 0 200]); xlabel { "Elapsed Time (sec)"); % 
TDCPMG.M - - - Simulate DElayed CPMG sequence % % NMR Simulation Program. % % 
Modification History: % May 2001 % %%%%%%%%%%%%%%%%%%%%%%% experimental parameters 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% echo on; 12. = le6; % X2. sec Assumed true 12 TE = le-3; 
% TE sec Echo time TD = 32e-3; % TD sec Delay time from pi/2 to 1st echo NE = 50; % 
NE Number of echoes t90 = 100e-6; % t90 sec 90-degree pulse width tl80 = 160e-6; % 
tl80 sec 180 -degree pulse width AT = TE-tl80; % AT sec Acquisition time dt = 20e-6; 
% dt sec Dwell time A = l/(2*tl80); % A Hz Pulse amp! i tudft fmin = -5000; % fmin Hz 
Leftmost freq sample fmax = 5000; % fmax Hz Rightmost freq sample df =40; % df Hz 
Frequency resolution f jitter =10; % f jitter Hz Random frequency jitter figured); 
echo off; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% t=0:dt:AT; % create acquisition time axis 
f =fmin:df :fmax; % create frequency axis [dum, NF] =size ( f ) ; % NF . . . no of frequency 
points sx=[]; % real channel signal sy=[]; % imaginary channel signal % establish 
equilibrium mx= zeros (size (f) ) ; my=zeros (size (f ) ) ; mz=ones (size (f ) ) ; % 90 pulse- 
phase = 90 [mx,my,mz] =pulsim ( j *A, t90 , f ,mx,my,mz) ; % TD/2 - pi pulse (phase=????) 
[mx,my,mz] =evolve (T2, t90+ (TD-t90-tl80 ) /2, f ,mx,my,mz) ; 
[mx,my,mz] =pulsim( -A, tl80, f ,mx,my,mz) ; % TD/2 - first echo 
[mrnx, mmy , mmz] =evolve (I2^.(TD-tl80-AT) /2, f ,mx / my,mz) ; 

[tmpx, tmpy] =observe (dt, At, f ,mmx,mmy, mmz) ; % compute response sx= [sx ; tmpx] ; sy= [sy 
,- tmpy]; s=sqrt(sx. 2+sy. 2); plot(t,sx); grid; xlabel ( "Time in sec after pi 
pulse."); drawnow; %backtrack to right after 1st pulse and evolve 
[mx,my,mz]=evolve(T2,TD./2+TE./2-tl80,f ,mx f my,mz) ; for nr=2:NE, %f printf ( "Echo 
%. Of .backslash. n",nr*1.0) ; % pi pulse; phase = 0 

[mx,my, mz] =pulsim(A, tl80, f , mx,my,mz) ; [tmpx, tmpy] =observe (dt, AT, f ,mx,my, mz) ; % 
compute response [mx, my,mz] =evolve (T2 , TE-tl80 , f , mx,my, mz) ; % evolve 
[mxl,myl,mzl] =evolve (T2^.(TE-tl80) /2, f ,mx / my,mz) ; % echo profile sx= [sx ; tmpx]; 
sy=[sy ; tmpy]; s=sqrt(sx. 2+sy. 2); plot(t,sx); grid; xlabel("Time in sec after pi 
pulse."); drawnow; % introduce frequency jitter to mask numerical discretization f = 
f +f jitter*randn(size (f ) ) ; end % create entire timeline tall=[]; for nr=l:NE, 
tall=[tall TD+{nr-l) *TE+( -AT/2 :dt: AT/2) ] ; end % plot entire train tmp=sx' ; 
sxall=tmp( : ) ; figure (2); plot (tall , sxall) ; axis([0 0.05 0 200]); xlabel ( "Elapsed 
Time ( sec ) ") ; ******************************************************************* 
function [mx,my,mz] = pulsim( amp, tau, f, mxO, myO, mzO ) % % .COPYRGT. Copyright, 
M. G. Prammer, NUMAR Corp., 1991. % % PULSIM - - - calculate hard pulse response. % 
[mx,my,mz] = pulsim( amp, tau, f, mxO, myO, mzO) % amp . . . real or complex pulse 
ampl -i hidft in Hz. % tau . . . pulse duration in sec. % f . . . frequency sampling 
points in Hz. % mx,my,mz . . . cartesian magnetization components. % Method: 
Riemannian spherical affine projection onto % complex plane, followed by complex 
transform. Yao, JMR. % Fails for Mz<= ( -l+le-16) . Written by MGP, Oct. 1991. % 
Frequency sign corrected for right hand reference frame. MGP, 12/17/91. M=sqrt(mx0. 
2+myO. 2+mzO. 2); % magnetization vector length mx0=mx0./M; my0=myO./M; mz0=mz0./M; 
% normalize Om = sqrt (abs (amp) . 2+f. 2); % effective rotation freq in Hz ct = 
-f./Om; % cos of rotation angle theta +++MGP (-) . st = sqrt(l-ct. 2); % sin of 
rotation angle theta ep = amp. /abs (amp) ; % phasor for argument of amp W = 
(mxO+j *my0) . / (l+mz0) . *conj (ep) ; % project onto complex plane cotom = 
j . /tan(Om*tau*pi) ; % needed for rotation W = ep.*( st- (ct-cotom) . *W) . / ( 
ct+cotom+st . *W) ; % perform rotation wf = real ( W.*conj(W)); % backtransf orm mx - 
M.*real( (W+conj (W) ) ./ (1+wf)); my = M. *real ( (W-conj (W) ) . / ( j * (1+wf ) ) ) ; mz = 
M. * ( 1-wf ) ./ (1+wf ) ; 

******************************************************************* f unc t i on 

[mx,my,mz] = evolve ( T2, tau, f, mxO, myO, mzO) % % .COPYRGT. Copyright, M. G. 
Prammer, NUMAR Corp., 1991. % % EVOLVE - - - evolve spin system in time. % 

[mx,my,mz] = evolve ( T2, tau, f, mxO, myO, mzO) W = (mxO+j *my0) . *exp ( -tau/T2 + 
j*2*pi*tau. *f ) ; mx = real (W) ; my = imag (W) ; mz = mzO; % if any (isnan(mx) ) , % 
fprintf (" -backslash. nError in EVOLVE ---");% keyboard; % end 

****************************************************************** f unct ion [ sx , sy ] 
= observe { dt, AT, f, mxO, myO, mzO) % % .COPYRGT. Copyright, M. G. Prammer, NUMAR 
Corp., 1991. % % OBSERVE - - - calculate observable signal. % [sx, sy] = observe { dt, 
AT, f, mxO, myO, mzO) % dt . . . dwell time, % AT . . . (total) acquisition time, % 
f . . . frequency axis, % mx, my, mz . . . magnetization before acquisition (NOT 
updated). % %t=0:dt:AT; %W= (mxO+j *my0 ) . *exp ( j *2*pi* t ' *f ) ; %sx=[] ; sy=[]; 
% [dum,nt] =size(t) ; %for k=l:nt, % sx= [sx sum (real (W (k, :)))] ; % sy= [sy 
sum (imag (W(k, :)))] ; %end sx=[]; sy=[]; W0=mx0+j *my0 ; J2PI=j*2*pi; for t=0:dt:AT, 
W=W0.*exp(J2PI*t.*f ) ; mx=real(W); my=imag(W); sx=[sx sum(mx)]; sy= [sy sum(my)]; end 
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CLAIMS: 

1. A method for making nuclear magnpfir rp.qnnanrp (MM7? ) measurements of a geologic 
formation using a NMR logging tool, comprising the steps of: (a) providing a static 
magnetic field in a volume of said formation; (b) applying a modified CPMG pulse 
sequence according to the expression: 

W.sub.i - .pi./2-TD.sub. j /2 - . pi . -TD . sub . j /2-echo- [TE/2- .pi . -TE/2 -echo] . sub . k where 
W.sub.i is a variable delay with i.gtoreq.l; TD.sub.j is variable echo spacing, 
j.gtoreq.l; TE/2 is the Carr-Purcell echo spacing, TD.sub.j >TE ; k=l , . . . , N; 
and (c) measuring the induced NMR echo signals. 

2 . The method of claim 1 further comprising the step of processing the induced NMR 
echo signals to derive petrophysical properties of the formation. 

3. The method of claim 1 further comprising the step of selecting the variable echo 
spacing TD.sub.j such that the application of the modified CPMG pulse sequence 
causes the loss of pre -determined fast diffusion components of the measured NMR echo 
signals . 

5. The method of claim 3, wherein TD.sub.j is selected to cause the loss of signals 
associated with the gas phase and the wafpr phase. 

10 . The method of claim 1 further comprising the step of pairwise subtracting nmr 
echo signals associated with different TD.sub.j 's. 
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11. The method of cla£flF7 further comprising the step dWseparating the gas, oil and 
wat-pr signal components. 



12. A method for making nuclear magnpt-ir rpsnnanrp (nmr) measurements of a geologic 
formation, comprising the steps of: (a) providing a static magnetic field in a 
volume of said formation; <b) applying oscillating magnetic fields according to at 
least one modified CPMG pulse sequence characterized by having at least one first 
echo spacing TD and a second echo spacing TE; wherein the at least one first echo 
spacing TD is selected to correspond to diffusion characteristics of fluids in the 
formation and cause corresponding amp! i tude loss in induced NMR signals, and TE is 
relatively short, such that diffusion in the corresponding induced nmr echo signals 
is substantially negligible; (c) measuring the induced NMR echo signals; (d) 
determining the amount of —amplitude loss resulting from at least one TD interval; 
(e) computing diffusion properties of fluids in the formation based on the 
Hpt-prmi rtPd ampl i bidft loss . 

14. The method of claim 12 further comprising the step of processing the induced NMR 
echo signals to derive petrophysical properties of the formation. 

16. The method of claim 12, wherein TD is selected to cause the loss of signals 
associated with the gas phase and the wat-e-r phase. 

20. The method of claim 12 further comprising the step of pairwise subtracting NMR 
echo signals associated with different TD's to separate different f In-iri phases 
present in the formation. 

21. The method of claim 20 further comprising the step of separating gas, oi 1 and 
wafpr signal components present in the formation. 

22. The method of claim 12 further comprising the step of estimating diffusivity 
valups for gas, wat-pr and oil phases in the formation, and selecting values for the 
TD and TE echo spacings based on the estimated diffusivity valug.g. 

23. The method of claim 22, wherein the selection of val nen for the TD and TE echo 
spacings is based on forward -mode ling of signal components to achieve maximum 
contrast between fin-id phases. 

24. The method of claim 12 further comprising the step of applying calibrations to 
eliminate stimulated echo effects that distort echo amplitudes. . 

25. The method of claim 12, wherein in step (e) echo trains for different TD values - 
are processed either simultaneously or in pairs. 

26. The method of claim 25, wherein echo trains for different TD values are 
processed in pairs, and only common time windows from the longest TD interval onward 
are considered in the step of processing. 

27. The method of claim 26 further comprising the step of subtracting of 
corresponding echoes to enhance a particular fln-i phase. 

28. The method of claim 27, wherein the gas phase is enhanced by subtracting the two 
shortest TD szalnes; the wai-er phase between the next TD pair, and the oil phase from 
the longest TD valnf* only without subtraction. 

29. The method of claim 12, further comprising applying correct ion_icalii£S that 
account for (a) incomplete polarization, (b) incomplete suppression by diffusion, 
and {c) hydrogen density. 

30. A system for making nuclear magnpti r rpannanrp CNMR) measurements of a geologic 

formation, comprising: (a) means for providing a static magnetic field in a volume 
of said formation; (b) means for applying oscillating magnetic fields according to 
at least one modified CPMG pulse sequence characterized by having at least one first 
echo spacing TD and a second echo spacing TE; wherein the at least one first echo 
spacing TD is selected to correspond to diffusion characteristics of fluids in the 
formation and cause corresponding ampl i tudft loss in induced NMR signals, and TE is 
relatively short, such that diffusion in the corresponding induced NMR echo signals 
is substantially negligible; (c) means for measuring the induced NMR echo signals; 
(d) means for determining the amount of amplitude loss resulting from at least one 
TD interval; and (e) means for computing diffusion properties of fluids, in the 
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formation based on the^retermined ampl i t.uAet loss. 
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DOCUMENT- IDENTIFIER: US 614 0817 A 

TITLE: Magnetic: resnnance well logging method and apparatus 
Abstract Text (l) : 

Technique and apparatus are disclosed for identifying characteristics of 
hydrocarbons, such as distinguishing light and hpavy formation o±l^_ in formations 
surrounding an earth borehole. The technique and apparatus can also be used for 
determining residual oil saturation, permeability, and bound water volume. An 
embodiment of the method of the invention includes the following steps: providing a 
nuclear magnet- ir rpsnnanrp logging device that is moveable through the borehole; 
providing, from the logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnet j r T-F>xnT\ar\r.f* spin echoes from the 
formations; processing the received spin echoes to produce a t. anh. ? distribution; 
computing a magnet- ir rpRonanrp porosity from the t Knh 7. distribution; deriving a 
measure of total porosity of the formations; and comparing the magneti c rp.qnnanrp 
porosity with the total porosity. 

Brief Summary Text (2) : 

This invention relates to apparatus and techniques for determining characteristics 
of earth formations surrounding a borehole and, more particularly, to a method and 
apparatus for nnHpa r magnetic; reao nanre borehole logging to determine 
characteristics including one or more of the following: the nature of formation 
hydrocarbons, residual r>i T saturation, permeability, and bound water volume. The 
invention also relates to MR logging with improved logging speed. 

Brief Summary Text {4) : 

General background of nuclear magnetic resonance well logging [also called magnet i c. 
rPRnnanrp (MR) logging] is set forth in copending U.S. patent application Ser. No. 
08/873,582, assigned to the assignee hereof, and in U.S. Pat. No. 5,023,551. 
Briefly, in nuclear magnetic- rpsnnanrp operation, the spins of nuclei align 
themselves along an externally applied static magnetic field. This equilibrium 
situation can be disturbed by a pulse of an oscillating magnetic field (e.g. an RF 
pulse) , which tips the spins away from the static field direction. After tipping, 
two things occur simultaneously. First, the spins precess around the static field at 
the Larmor frequency, given by . omega sub . 0 = .gamma. B. sub. 0, where B.sub.O is the 
strength of the static field and .gamma, is the gyromagnetic ratio. Second, the 
spins return to the equilibrium direction according to a decay time Tl, the spin 
lattice rel ayation time. For hydrogen nuclei, .gamma. /2 .pi .=4258 Hz/Gauss, so, for 
example, for a static field of 235 Gauss, the frequency of precession would be 1 
MHz. Also associated with the spin of molecular nuclei is a second relaxation, T2 , 
called the gpin-gpin relaxation time. At the end of a ninety degree tipping pulse, 
all the spins are pointed in a common direction perpendicular to the static field, 
and they all precess near the Larmor frequency. However, because of molecular 
interactions, each nuclear spin precesses at a slightly different rate. 12 is a time 
constant of this "dephasing" . 

Rrief Summary Text (5) : 

A widely used technique for acquiring_UMR data both in the laboratory and in well 
logging, uses an RF pulse sequence known as the CPMG (Carr-Purcell-Meiboom-Gill) 
sequence. [See Meiboom, S., Gill, D., 1958, "Modified Spin Echo Method for Measuring 
Nuclear Pelaxat: ion Times", Review of Scientific Instruments, 29, 688-691]. As is 
well known, after a wait time that precedes each pulse sequence, a ninety degree 
pulse causes the spins to start precessing. Then a one hundred eighty degree pulse 
is applied to keep the spins in the measurement plane, but to cause the spins which 
are dephasing in the transvprsp plane to refocus. By repeatedly reversing the spins 
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using one hundred eighty degree pulses, a series of "spHF echoes" appear, and the 
train of echoes is measured and processed, for example to obtain a T . nub . 2 
distribution of fhn d components of the formations. 

RHpf Summary TPvf (6) : 

Magnetir ypRnnanrp logging has added a new dimension to formation evaluation. The MR 
measurement is sensitive to total f inid content, to the intrinsic properties of the 
fluids, and to the environments in which the fluids reside in the pore space of 
porous rock. There are many applications of magnetic reso nanrp. tools, and among the 
most common and important are the determination of capillary and clay bound water 
volumes and the estimation of permeability. Bound water analysis is of central 
interest to the log interpreter because it helps predict the potential water cut of 
a formation volume, and is a major factor in the perforation decision. A continuous 
permeability log can provide an estimation of potential production rate and also 
indicates flow barriers and thief zones. These log outputs are difficult to obtain 
reliably with other logging tools, and hence may be termed "MR-unique". 

Brief Summary Text (7) : 

One limitation of MR tools is their slow logging speed. Magnetic reso nanrp logging 
uses a cyclic measurement consisting of a wait time followed by an echo acquisition 
period. Conventional MR data are presently acquired with 

Brief Summary Text (8) : 

wait times sufficiently long to substantially polarize all fluid protons in the 
formation. Protons in gas, light oil, oil base mud filtrates, and vug water polarize 
very slowly. Therefore, very long wait times have been used when logging formations 
in which those fluids are present. The long wait times necessarily contribute to 
slow logging speeds. Thus, while MR tools provide important additional information, 
it can be at the cost of wellsite efficiency. 

Brief Summary Tftyf. (10) : 

Heretofore it has been considered difficult or impossible to estimate residual oil 
saturation in wells drilled with nil base mud. This is particularly true when the 
native nil ha« a low downhole viscosity, which renders it difficult or impossible to 
distinguish native fluids from invaded fluids on the basis of existing magnetic? 
resonance T.mih.2 measurements. 

Brief Summary Text. (11) : 

It is among the objects of the present invention to overcome limitations of prior 
art MR techniques by determining residual ail saturation and other formation 
characteristics in wells drilled with ni 1 base mud. It is also among the objects- of 
the invention to improve logging speed (e.g. for a given precision and/or vertical 
resolution) of MR logging. It is also among the objects of the invention to provide 
a technique for determining when formation permeability may be overestimated in 
certain types of formations. 

Brief Summary Text (14) : 

Magnetic resonanpp tools are extreme examples of cyclic tools: the measurement cycle 
is longer and more complicated than for any other logging tool. In a typical MR 
measurement cycle, and as referenced in the background portion hereof, permanent 
magnets within the tool polarize the hydrogen nuclei of formation fluids during a 
"wait time" . Then, an antenna broadcasts a series of pulses of hxgh frequency 
magnetic field. Return spin echoes from the formation are detected after the second 
and each subsequent pulse in a CPMG pulse sequence. 

RHpf Summary Text (15) : 

In a known type of MR logging device to be referenced hereinbelow, the 
wait -acquisition cycles are collected in phase alternated pairs (PAPs) . A phase 
shift of the RF carrier of one pulse in the sequence reverses the sense of the CPMG 
acquisition [See Sezginer, A., Kleinberg, R. L. , Fukuhara, M. , Latour, L. L., 1991, 
"Very Rapid Simultaneous Measurement of Nuclear Magnetic Resonance. Spin-Lattice 
Relaxation Time and Spin-sp-in Relaxation Time", Journal of Magnetic Resonance 92, 
504-527; U.S. Pat. No. 5,023,551.] Taking the difference between the positive and 
negative CPMG's echo-by-echo doubles the size of the signal and eliminates baseline 
and several spin-dynamics errors. 

Brief Su mmary Text (19) : 

The sandstone logging mode is generally satisfactory when the formation hydrogen can 
be polarized during the 1.3 second wait time. Then the logging tool can provide a 
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complete suite of MR measurements, including porosity, ^Pre size distribution, and 
other useful log outputs. However, in some circumstances, the 1.3 second wait time 
is inadequate to fully polarize the formation. This is the case when formations have 
long longitudinal rpl avatinn times, T.sub.l. As noted in the background hereof, 
examples are reservoirs with light ni 1 or gas, those drilled with oil base muds, or 
those with vugs. In such cases, T.sub.l is usually equal to its bulk fluid value, 
which can be several seconds. [See Kleinberg, R. L. Vinegar, H. J., 1996, "NME 
Properties of Reservoir Eliiids" , Log Analyst November -Dec ember 1996. pg. 20-32.] The 
wait time must be several times longer than the longest formation T.sub.l to obtain 
complete information. This increases the cycle time, T. To maintain precision and 
vertical resolution, the logging speed must be reduced. See, for example, Morriss, 
C. E . , Deutch, P., Freedman, R. , McKeon, D. , Kleinberg, R. L., 1996, "Operating 
Guide for the Combinable Magnetic: Resonance. Tool", Log Analyst, November -Dec ember 
1996, pg. 53-60, which describes an MR tool operating in an oil base mud mode with a 
logging speed of 200 ft/hr. 

Brief Summary Text (20) : 

Magnetic rpflnnanrp logging tools are characterized by considerable flexibility in 
measurement modes. Measurement can be tailored to the objectives of the logging 
program. In an illustrated embodiment and examples hereof, wells drilled with oil 
base mud are considered. Neutron, density, resistivity and gamma ray tools can 
efficiently provide porosity, lithology, and deep saturation information. 
Permeability and precise_jcalii£S of bound wai-pr volume used for the determinations 
hereof are obtained h y magnetic rpsnnanrp logging. Rig time expense dictates the 
most rapid acquisition possible consistent with requirements of enhanced precision 
and good vertical resolution, and rapid acquisition is an advantageous feature of 
the invention. 

Rripf S ummary Text (21) : 

As is well known, the volume of bound wat.pr at a given air-brine capillary pressure 
is associated with mirlpa r magnerir? rpRn nanrp signals that decay with a T . .qnb . 2 less 
than a format ion -dependent T. sub. 2cutof f . In sandstones, for an air-brine capillary 
pressure of 100 psi (700 kPa) , the T. sub. 2cutof f is often about 33 msec, but can be 
higher or/ lower [see Straley, C. , Rossini, D, Vinegar, H., Tutunjian, P., Morriss, 
C, 1997, "Core Analysis by T,ow v\ r1 d nmt?" r Log Analyst, March-April 1997, pg. 
84-94] . The T.sub.l associated with this part of the porosity is often about 1.5 
times longer than t ruT^ 7. [see Kleinberg, R. L. , Straley, C, Kenyon, W. E., Akkurt, 
R., Farooqui, S. A., 1993, "Nuclear Magnetic Resonance of Rocks: T.sub.l vs. T . sub . 2 
", Society of Petroleum Engineers Papers 26470; Straley, C, Morriss, C. E., Kenyon, 
W. E., Howard, J. J., 1995, " NMR in Partially Saturated Rocks", Log Analyst, 
- January- February 1995, pg. 40-56]. A reduced wait time, of duration a few times 
longer than longest T.sub.l of the bound fluid, will completely polarize the bound 
wafpr in the formation. Note also that fewer echoes are needed to characterize 
r-el avafinn time components shorter than T. sub. 2cutof f than are required to 
characterize the entire t guri ? distribution. 

Brief Summary Text. (22) : 

When the wait time is short, those formation fluids having long T.sub.l ' s are not 
completely polarized, i.e. they are not fully represented in the subsequent CPMG 
echo train. As is known in the art, for a fluid component characterized by T.sub.li, 
the apparent signal ampl i hirip is given by ##EQU2## where W is the wait time, A.sub.i 
is the signal strength (in porosity units) associated with a formation f 1 ui d i 
having polarization time T.sub.li, and A.sub.Oi is the signal strength for that 
fluid in the limit W. f wdarw. . inf in. . if W is not several times longer than T.sub.li, 
the signal strength A.sub.i is not representative of the volume of that fluid. 

Rrief Summary Text (23) : 

Laboratory measurements on water saturated rocks have shown that T.sub.l can be 
estimated by T.sub.l = .xi . T . sub. 2 , where l<.xi.<3. Then, equation (3) can be written 
as ##EQU3## Conventional processing with the type of MR tool described herein 
automatically accounts for this factor. The result is functionally equivalent to 
correcting the T . sub . 2 distribution; that is ##EQU4## where P (T . sub. 2i) . sub.uncorr 
is the uncorrected t snh . 2 distribution and P(T.sub.2i) is the displayed 
distribution. The default value, .xi.=1.5, is appropriate for many water saturated 
rocks . 

Brief Summary Text (24) : 

In the presence of light oil or oil base mud, when the T - suh . 2 distribution is 
distorted by the effects of reduced cycle time, the standard polarization correction 
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is inappropriate. It ^rpref erable to use equation (3) W recover the true amp! i t-nrip 
of fluid components with long T.sub.l 's. The polarization correction would then be 
##EQU5## Note that a fixed vai»p of T.sub.l --that is, the longer rpiavatinn time of 
the ni 1 phase- -is used in this equation. This correction should only be applied to 
that part of the t. mih - 7 distribution which represents the oi 1 phase . The bound 
wat-er fin water wet rock) has a very short rpl ayatinn time, and clearly the T.sub.l 
of the o-i 1 phase should not be applied to it. T.sub.l is not measured directly by MR 
tools, but in many cases it can be estimated from the bulk properties of the fluids 
at reservoir conditions [see Kleinberg, R. L. Vinegar, H. J., 1996, "NMR Properties 
of Reservoir Elnida" , Log Analyst November -Dec ember 1996. pg. 20-32]. 

R-Hpf Summary Text (25) : 

Cycle time can also be reduced by reducing the number of echoes acquired in each 
CPMG. A long acquisition time, N.times.TE, is required if good definition of long 
t . snh 7 components is required. To obtain accurate values of bound water volume, it 
is only necessary to resolve components having t anh.9 less than T . sub . 2 cutof f . For 
TE=0.32 msec, as few as 100 echoes are adequate to define the T.snh 7 distribution 
below 30 msec. 

Brief Summary Text (27) : 

In accordance with an embodiment of the invention, there is disclosed a method for 
identifying characteristics of hydrocarbons (e.g. distinguishing light and heavy 
formation oil) in formations surrounding an earth borehole, comprising the following 
steps: providing a nuclear magnpt- ir rpannanrp. logging device that is moveable 
through the borehole; providing, from the logging device, a polarizing magnetic 
field and cycles of a magnetic pulse sequence, and receiving magnetic: resonance spin 
echoes from the formations; processing the received spin echoes to produce a T ■ nub . 2 
distribution; computing a magnpfir resonance porosity from the T. sub. 2 distribution; 
deriving a measure of total porosity of the formations; and comparing the magnet i c 
rpffnnanrp porosity with the total porosity. Residual_oil saturation can also be 
obtained as a function of the comparison. 

RT"if>f Summary Tftvt (28) : 

In a preferred embodiment of the invention, the polarization wait time, W, of said 
cycles of magnetic pulse sequence is short enough to incompletely polarize the 
protons in an investigation region of the formations. In this embodiment, the wait 
time, W, is in the range 0.1 second to 1 second. Also in this embodiment, the step 
of comparing the magnet. -i f rpsnnanrp porosity with the total porosity comprises 
comparing the respective porosities in a jarater- bearing zone of the formations and 
comparing the respective porosities in an oil-bearing zone of the formations. In a 
form of the preferred embodiment, the following steps are also performed: developing 
a free/bound fluid mit-.nff of thp T.mih.a distribution and applying a polarization 
correction to the free fluid components of the t anh 7 distribution before computing 
the magnpfir rpsnnanrp porosity from the T , suh ■ 2 distribution. In this form of the 
preferred embodiment, the step of comparing the magnetic: resonance porosity with the 
total porosity comprises comparing the respective porosities in a watex- bearing zone 
of the formations and comparing the respective porosities in an nil -bearing zone of 
the formations . 

Drawing Description Text (3) : 

FIG. 2 is a cross -sectional diagram of a portion of the nuclear magnetic resonance 
logging device of FIG. 1. 

Drawing Description Text (4) : 

FIG. 3 is a perspective view of the RF antenna of the FIG. 1 nuclear magnetic 
rpgnnanrp logging device. 

Drawing nPRrription Text (7) : 

FIG. 6 shows t g^h 7. distributions [graph (a)] and cumulatives of t gnh 7 
distributions [graph (b) ] for laboratory NMR measurements on a mixture of hexane and 
nickel chloride for different wait times. 

Drawing Description Text (8) : 

FIG. 7 shows nmr fluid volume (normalized by true volume) versus T.sub.l for three 
values: of wait time W. 

Detaile d Dpnrripti nn Tpxt (2) : 

Referring to FIG. 1, there is shown an apparatus for investigating subsurface 
formations 31 traversed by a borehole 32, which can be used in practicing 
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nv^ntion. The borehole 32 is typicJ^^ 



embodiments of the invention. The borehole 32 is typicifWy filled with a drilling 
f 1 ni d or mud which contains finely divided solids in suspension, and mudcake 16 is 
shown on the walls of the borehole. The invention has particularly advantageous 
application when the mud is oi 1 base mud. 

DRtailpd Dpficr ipt-i nn Tpyt (3) : 

A magnRt-ir; rpsnnancp investigating apparatus or logging device 30 is suspended in 
the borehole 32 on an armored cable 33, the length of which substantially determines 
the relative depth of the device 30. Also shown in FIG. 1 as part of the logging 
string are further logging devices represented by reference numerals 60, 70, and 80. 
In the present embodiment the device 60 represents a density and/or neutron logging 
device (although it may alternatively be a sonic device or other logging device for 
obtaining total porosity) , the device 70 is a resistivity logging device, and the 
device 80 is a gamma ray logging device. It will be understood that measurements 
from these and/or other logging devices can, if desired, be obtained from separate 
logging runs with these and/or other logging devices or from other sources. The 
length of cable 33 is controlled by suitable means at the surface such as a drum and 
winch mechanism (not shown) . Surface equipment, represented at 7, can be of 
conventional type, and can include a processor subsystem and communicates with the 
all the downhole equipment. It will be understood that processing can be performed 
downhole and/or uphole, and that some of the processing may be performed at a remote 
location. Also, while a wireline is illustrated, alternative forms of physical 
support and communicating link can be used, for example in a measurement while 
drilling system. As described for example in the U.S. Pat. No. 5,055,787, the 
magnptir rftaonannp logging device 30 has a face 14 shaped to intimately contact the 
borehole wall, with minimal gaps or standoff, and a retractable arm 15 which can be 
activated to press the body of the tool 13 against the borehole wall during a 
logging run, with the face 14 pressed against the wall's surface. Although the tool 

13 is shown as a single body, the tool 

Dpta-il^d Description Text (5) : 

The magnetic rpsnnanre logging device 30 includes a magnet array 17 and an RF 
antenna 18 positioned between the array 17 and the wall engaging face 14. Magnet 
array 17 produces a static magnetic field B.sub.O in regions surrounding the tool 
13. The antenna 18 produces, at selected times, an oscillating magnetic field 
B.sub.l which is focussed into formation 12, and is superposed on the static field 
B.sub.O within those parts of formation opposite the face 14. The "volume of 
investigation" of the tool, shown in dotted lines in FIG. 3, is a vertically 
elongated region directly in front of tool face 14. As described in the referenced 
Patent, the tool 13 can make measurements by magnetically tipping the nuclear spins 
of particles in formation 12 with a pulse of oscillating field B.sub.l, and then 
detecting the precession of the tipped particles in the static field B.sub.O within 
the volume of investigation over a period of time. Reference can also be made to 
Morriss, C. E . , Deutch, P., Freedman, R. , McKeon, D., Kleinberg, R. L. , 1996, 
"Operating Guide for the Combinable Magnphir Rpsnnanrp Tool", Log Analyst, 
November -December 1996, pg. 53-60. 

npfailpH nRKrrription Text (6) : 

FIG. 2 shows a magnet array 17 of the type disclosed in an embodiment of the 
referenced '787 Patent. The magnet array includes two permanent magnets 24 and 26, 
which are mounted generally parallel to each other within a metal alloy body 27. The 
body 27 should be of a material having "1 r>w magnetic permeability, so as to not 
interfere with the static magnetic field. Magnets 24 and 26 are slabs which are 
elongated in the longitudinal direction of the borehole. The magnetic poles of each 
magnet are not on the smallest faces of the slab, commonly viewed as the ends of a 
bar magnet. Instead, the poles appear on the two opposing edges of the slab magnet 
and point to the left and right, respectively, in the Figure. Therefore, within the 
formation 12, the. magnetic field B.sub.O surrounding the magnets remains fairly 
constant along the longitudinal direction of the borehole axis. In the illustration 
of FIG. 2, magnets 24, 26 are symmetrically mounted in the two sides of the body 27 
with the north poles facing the same direction, that is, the direction of the face 

14 of the tool. One or more further permanent magnets can be used. 

Detailed Deficr ipl-ion Tpyt (8) : 

In the referenced '787 Patent, the antenna 18 is used both as an RF transmitter to 
produce an oscillating magnetic field in the formations, and as a receiving antenna 
to detect coherent magnetic signals emanating from precessing protons (spins) after 
the oscillating field is terminated. The antenna, which has a body 29 and an 
elongated center probe 42, across which signals are applied and detected, serves 
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effectively as a currenx loop which produces an oscillaWhg magnetic field B.sub.l 
(see FIG. 4) within the volume of investigation that is perpendicular to the static 
magnetic field, B.sub.O (which is radial in the volume of investigation). The body 
29 is trough-shaped and has end plates 40, 41 with the center conductor or probe 42 
extending from one end plate 40 to the other end plate 41, parallel to and centered 
in the semi -cylindrical trough 29. The U.S. Pat. No. 5,153,514 discloses that the 
trough antenna, which can be filled with a ferrite, can have an inner conductive 
shell that is separated from a steel body by a rubber layer, which suppresses 
magnetoacoustic ringing. It will be understood that various other types of magnetic 
rpsnnanrp logging equipment can be used in practicing the invention. 

Detailed Description Text (10) : 

To illustrate the effect of cycle time reduction, roughly equal amounts of hexane 
(T.sub.l = T. sub. 2 =1 . 80 sec) and water doped with NiCl (T.sub.l =T. sub. 2 =0.022 sec) 
were placed together in a sample tube and nuclear magnetic resonance measurements 
performed. The results are shown in FIG. 6. Run #1 (solid line) used wait=4 sec and 
1800 echoes, Run #2 (dotted line) used wait=0.25 sec and 300 echoes. Graph (a) shows 
t snh 7 distributions and graph (b) shows cumulatives of the distributions. No 
polarization correction was applied. The t. snh. s> distribution found from the long 
cycle time measurement shows two distinct peaks corresponding to water and hexane. 
The short cycle time measurement underestimates the quantity of long-re 1 axati on-time 
fluid while capturing all the information associated with short rplavatinn time 
fluid. The t snh. o. distribution for the hexane (like an oil signal) is considerably 
distorted when only 300 echoes are collected, as this number of echoes is not 
sufficient to define the details at the long-T. suh.^ end of distribution. The 
Rt-andard polarization correction of Equation (5) can yield erroneous results unless 
artificial valuer of .xi. are used. Instead, the fixed T.sub.l correction of 
Equation (6) should be applied to the 1 nng-T snh 2 part of the rel axati on time 
distribution . 

TV»1-*-ilRri nf>gm-ipf.ion Text (11) : 

On a local basis, permeability sometimes correlates with porosity. In general, 
however, permeability depends on rock fabric, which can be independent of porosity. 
MagnRtir rpsnnannp. is an accepted method of probing rock microgeometry, and a number 
of permeability relations based on magnpf irr -resonance have been found useful as 
follows [see Morriss, C. E. , Deutch, P., Freedman, R. , McKeon, D. Kleinberg, R. L. , 
1996, "Operating Guide for the Combinable Magnetic: Resonance Tool", Log Analyst, 
November -December 1996, pg. 53-60] ; 

Deta i 1 ed De.qrrr-i pt i on Text (12) : 

Permeability can also be estimated by* Equation (8) from short cycle MR measurements, 
using FFI= .phi . -BFV, where the bound fluid volume is measured by magnetic resonance, 
and porosity is measured by other tool(s) . ##EQU7## There is some ambiguity in 
exactly which porosity is to be used; this point having been discussed by Singer et 
al. [see Singer, J. M. , Johnston, L. , Kleinberg, R. L. , Flaum, C, 1997, "Fast NMR 
Logging for Bound Fin-id and Permeability", SPWLA 38th Annual Logging Symposium, 
Paper YY] . In the present embodiment, the porosity used is density porosity 
.phi.. sub. D, or, in the presence of gas, the neutron-density crossplot porosity 
.phi. .sub.T. 

ngtiailpd Descri ption Text: (13) : 

Recent advances in MR tools, measurement techniques, and signal processing [see e.g. 
Prammer, M. G. , Drack, E. D., Bouton, J. C. , Gardner, J. S., 1996, "Measurements of 
Clay Bound Hater and Total Porosity by Magnshi r Rpsnnanrp Logging", Long Analyst, 
November-December 1996, pg. 61-69; Freedman, R. , Boyd, A., Gubelin, G. , McKeon, D. , 
Morriss, C. E., Flaum C, 1997, "Measurement of Total NMR Porosity Adds New Value to 
nmr Logging", SPWLA 38th Annual Logging Symposium, Paper 00] have dramatically 
improved the measurement of bound fluid volumes. Now, fluids having T.suh.2 
rplavafinn times of 0.5 msec or lower are included in the bound fluid volume (BFV) . 
With these improvements, MR porosity agrees with density porosity in most shale 
intervals. Since all shale porosity is bound fluid, it is usually found that 
.phi.. sub. D =BFV in shales. Thus, Equation (9) estimates zero permeability in 
typical shales. 

n^l-ailed neptrri pti on Text (14) : 

Despite recent improvements in measuring porosity components with very short 
t gnh o. f magnpUr rpan-nanrp* measurements can still underestimate bound fluid volume. 
For example, it has been observed that some mature, dewatered shales have rpi ayai-i on 
time components below 0.3 msec. Bitumen volumes are also considerably underestimated 
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by MR measurements. Then .phi.. sub. D >BFV, even when ttierB is no free fluid. Under 
these circumstances, Equation (9) will erroneously indicate the presence of 
substantial permeability. 

Detailed npHcHpfinn Text (15) : 

In order to automatically detect formations in which permeability is likely to be 
overestimated, a t,ow BFV Flag can be computed as: ##EQU8## Here FFI.sup. (-) is the 
underestimatecL_mliie of free fluid volume. FFI is underestimated in reduced cycle 
time logging because short wait times result in partial polarization of fluid 
components having T.sub.l >W. .epsilon. is a small constant {typically 0.25 p.u.) 
which prevents the denominator from going to zero in tight zones and shales. 

Detailed Description Text (16) : 

The flag is zero in tight zones and in typical shales. In these zones, no correction 
to the computed permeability is required. In most other formations; .phi.. sub. D 
-BFV. about .FFI. sup. (-) , so the Flag is near unity, which is denoted "low"; no 
correction to the computed permeability is required. In dewatered shales and 
bitumen, .phi.-BFV>0 and FFI.sup. (=) 0, so the Flag is high. In these zones the 
permeability can be automatically set equal to an arbitrarily low value. A flow 
diagram for this processing is described in conjunction with FIG. 9 below. 

Dpfailpd Descr iption Text (17) : 

Oi 1 a of different viscosities have different nuclear magnetic rpsn nanrp rplavaf inn 

times [see Kleinberg, R. L. Vinegar, H. J., 1996, "NMR Properties of Reservoir 
Fluids n , Log Analyst November -Dec ember 1996. pg. 20-32]. T. mih. 2 distributions from 
the type of MR logging device described herein are routinely used to determine the 
viscosity, .eta., of oil in place, and to distinguish productive oil zones from 
heavy and immobile hydrocarbons, such as in tar mats. Moreover, oi 1 -water contacts 
in wells drilled with watex base mud can frequently be unambiguously identified 
using t _ snh . 7. measurements [see White, J., 1997, "Applications of Downhole Magnet i c 
Bpsnnanrp Imaging in the North Sea", Society of Petroleum Engineers Offshore Europe 
Conference Paper 38551] . 

Detail pri Description Text (18) : 

Finding the oi 1 -water contact (OWC) poses a challenge to magnet i c resonance well 
logging when the well is drilled with ail base mud and the formation oi 1 is light. 
In the watex leg,_oil base mud (OBM) filtrate displaces movable water in the volume 
of investigation of the MR tool, while in the oil leg, the filtrate mixes with 
formation oil. Because the nuclear magnetic resonance T.suh.2 is relatively 
insensitive to viscosity when .eta. . Itoreq. 1 cP, the T . snh. 2 distribution often 
shows no change at the OWC . - 

Detailed Description Text (19) : 

However, T.sub.l remains sensitive to viscosity over the entire range of crudes and 
mud filtrates, as shown in FIG. 11 of Kleinberg, R. L. Vinegar, H. J., 1996, "NMR 
Properties of Reservoir Fluids " f Log Analyst November-December 1996. pg. 20-32. 
Thus, T.sub.l can be used to discriminate OBM filtrate from formation oil. The 
apparent wmr porosity . phi . . sub . app is the sum of contributions from the fluids in 
the formation. For true porosity .phi. and saturations S.sub.j ##EQU9## where W is 
the wait time and T.sub.lj is the relaxation time of fluid j. If the wait time is 
substantially longer than the T.sub.l ' s of all f in-i ds, . phi . . sub . app =.phi., so the 
magnetic resonance signal ampl i tude will be unchanged as the tool crosses the OWC. 
On the other hand, reduced wait times can result in substantial signal ampl i tude. 
changes. When combined with other porosity logs, use of a single MR pass will 
differentiate OBM-invaded wai-er and oil zones. 

Detailed Description Text (20) : 

FIG. 7 illustrates a basis of the technique. The factor in square brackets in 
Equation (11) is plotted versus T.sub.l for three wait times: 8 sec, 1.3 sec, and 
0.45 sec. Even 8 sec is not adequate to fully polarize fluids with T.sub.l>2 sec, 
which includes many OBM filtrates and light oils . Using any of the three wait times 
will produce ampl itnde contrast between fluids of differing T.sub.l. While W=0.45 
sec is shorter than optimal for the OWC application, this aggressively reduced wait 
time permits faster logging, improved vertical resolution and/or increased data 
stacking over the entire logged section while yielding good ampl i t.urie contrast at 
the OWC. 

Detailed Descri ption Text (21) : 

Referring to FIG. 8, there is shown a flow diagram of a routine for controlling a 
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processor in accordar^^ with an embodiment of the inve^^on. The processor can be a 
downhole processor, the uphole processor in equipment 7 (FIG. 1) , or a remote 
processor which operates subsequently on stored and/or transmitted 

log- representative signals, as well as combinations of downhole, uphole, and remote 
processors. In the flow diagram of FIG. 8, the block 805 represents selecting of the 
wait time W. The wait time will be relatively short for incomplete polarization. The 
preferred range of W is 0.1 to 1 second. Next, as represented by the block 815, the 
echo spacing, TE, and the number of echoes are selected for the pulse sequence. The 
pulse sequence is then applied and the resultant echo train is recorded. It will be 
understood that any suitable technique can be used for applying the pulse sequence, 
a CPMG pulse sequence being used for the present embodiment. A T. sub. 2 distribution 
is then generated from the echo train data, as represented by the block 840, and a 
ahanrfarrf free/bound fin-id cutoff is developed, as represented by the block 850. 
Reference can be made, for example, to Straley, C, Rossini, D. , Vinegar, H., 
Tutunjian, P., Morriss, C., 1997, "Core Analysis by Low Field NMR" , Log Analyst, 
March-April 1997, pg. 84-94. The block 860 represents application of correction to 
the free fin-id components of the T . nub . 2 distribution, which are affected by the 
incomplete polarization, as previously described. In the present embodiment, the 
correction is in accordance with Equation (6) . The density (or neutron-density) 
porosity, such as measured by device 60 (FIG. 1) is indicated as being available 
(block 1055) , and can be utilized in a wafpr zone to select the correction such that 
corrected MR porosity will match the total porosity .phi.. sub. T in the wafpr zone 
flushed with oil base mud filtrate. The decision block 870 represents the comparison 
of the magnPt-.-ir rpsnnanrp porosity with the total porosity in an oil bearing zone. 
In a zone of light hydrocarbons, the MR porosity will be less than the 
neutron-density derived total porosity (since the applied correction will be too 
small) , whereas in a zone of h^avy n-i 1 the MR porosity will be greater than the 
neutron- density derived total porosity (since the applied correction will be too 
large) . The block 875 represents the determination of residual oil saturation, which 
can be implemented e.g. by quantification of the comparison. Reference can be made 
to equation (A-10) below and the derivation thereof. The block 880 represents 
determination of permeability, e.g. by using Equation (9) . 

Detailed Description Text. (22) : 

FIG. 9 is a flow diagram of a routine for programming a processor in accordance with 
an embodiment of the invention for producing and using the Flag for determining 
potentially overestimated permeability in certain formations wherein at least some 
of the bound water is not "seen" by the MR logging technique. Again, the processor 
can be a downhole, uphole, or remote processor. The blocks 8 05 through 850 (only 
block 850 being shown in FIG. 9) can correspond to their counterparts in the flow 
diagram of FIG. 8, with a short wait time again being used. After the free/bound 

Detail pd DgRcripl-inn Text (23) : 

fluid cutoff is developed (block 850), the volumes of bound fluid (BFV) and free 
fluid (FFI .sup. ( -) , which, as noted above, is the uncorrected and possibly 
underestimated volume of free fluid) are computed from the respective components of 
the t gun •? distribution (block 910) . Using the computed valnfis, and the porosity 
from block 855, the Flag can be computed (block 920) . Determination is then made 
(decision block 930) as to whether the Flag is above a predetermined t-hrgghold. If 
so, the zone can be marked as one of potentially overestimated permeability (block 
940) . 

DRtaile d nPRnription Tpyf (25) : 

In examples hereof, reduced cycle time logging is employed. The reservoirs of the 
examples were drilled with oi 1 base mud and logged using an MR logging tool of the 
type described herein (and designated as the "CMR" , which is a mark of Schlumberger) 
and other logging devices as in FIG. 1. [Regarding the CMR tool and processing, 
reference can be made to U.S. Pat. Nos. 5,055,787 and 5,2 91,137 and to Freedman, R. , 
Morriss, C. E., 1995, "Processing of Data from an NMR Logging Tool", Society of 
Petroleum Engineers Paper 30560".] Log Example #1 was a well in the Andrew field in 
the North Sea. A logging objective was to acquire data at high speed and with 
enhanced vertical resolution. Neutron, density, resistivity and gamma ray logs were 
also run. Using the CMR tool to find yet another porosity measure is inefficient 
because formation nil and oil base mud filtrate both have long T.sub.l, so very slow 
logging would have been required. Even so, detailed log analysis shows that 
comparing reduced- cycle -time CMR apparent porosity to neutron-density porosity 
allows precise determination of the residual oi 1 contact. The CMR logging device and 
pulse sequence parameters were (with SI being sampling interval) : 
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net-ailed Description TExt (26) : 

The porosity analysis is as follows: The porosity logs are shown in FIGS. 10A and 
10B, track 3. Density and neutron tools were used to determine porosity assuming a 
sandstone matrix and a f Inid density of 0.9 g/cc. Log porosity is 

Detail e ri Descri ption Text (27) : 

The second term comes into play only in the presence of gas, where .phi.. sub. D 
< .phi . . sub.N. This porosity is in excellent agreement with core measurements, shown 
as dots, water volume, .phi.S.sub.W is derived from deep resistivity (90 inch 
array). R.sub.W =0.07 .OMEGA. -m was selected so that resistivity-derived porosity 
matched crossplot porosity in the water leg below x760 feet: ##EQU10## 

Detail pd Description Text (28) : 

Porosity derived from the CMR tool is also shown, as is the t. Rnh. 9 distribution 
from 0.3 msec to 3000 msec. The wait time of 0.45 sec was insufficient to fully 
polarize the native oil or oil base mud filtrate. Therefore the standard 
polarization correction was turned off and a fixed value of T.sub.l =2.5 sec was 
used in a custom polarization correction using Equation (6) . This correction was 
applied only to the free fluid part of the signal. The value of T.sub.l was picked 
so that the corrected CMR porosity matched the density porosity in the water leg 
(flushed with oi 1 base mud filtrate) below x760. The excellent correlation can be 
noted . 



Detailed Description Text (29) : 

The bound fluid volume, BFV, is shown in Track 3. The formation is a very clean 
sand, in which there is little bound water. The T. sub. 2cutof f =100 msec was selected 
so that BFV matched .phi.S.sub.W high in the hydrocarbon zone. The one parameter fit 
of CMR data yields excellent curve -character agreement with irreducible water volume 
determined by resistivity above the transition zone, FIG. 10, Track 3. A 
conventional sharp nifnff was used. The use of a tapered cutoff [see Coates, G. R. , 
Marschall, D. , Mardon, D., Galford, J., 1997, "A New Characterization of Bulk Volume 
Irreducible Using Magner i r Bpgnnanrp " f SPWLA 38th Annual Logging Symposium, Paper 
QQ; Kleinberg, R. L. , Boyd, A., 1997, "Tapered Cutoffs for Magnetic Resonance Bound 
Water Volume", Society of Petroleum Engineers Paper 38737] would not be appropriate 
here, because most of the CMR signal is from the non-wetting phase (oil and oil base 
mud filtrate) and therefore is not reflective of the pore size distribution of the 
rock. 



Detailed Description Text (30) : 

In the hydrocarbon analysis, the CMR log, combined with other logs, reveals five 
distinct zones (FIG. 10) . (1) Below x760 the sand is water saturated; oil- base mud 
filtrate has displace d water in the CMR tool volume of investigation. (2) The 
original ( "paleo" )_oil^s£atex contact is marked by heavy oil bands at x760 and x750. 
(3) A layer of dead residual oil f remaining after partial depletion of the 
reservoir, is found between x740 and x705- (4) Live crude (mixed with OBM filtrate 
in the flushed zone) is found between x705 and xS30. (5) Gas mixed with OBM filtrate 
is found above x530. 

Detail pd Description Text (31) : 

In the wahpr sand, the free fluid is nil base mud filtrate which has displaced 
movable cnnnar e water . OBM filtrates are often alkanes of moderate chain length 
(e.g. hexadecane, C16) , and lack dissolved gas. As explained above, porosity from 
the CMR tool matches nuclear porosity in this interval when T.sub.l =2.5 sec is used 
in the polarization correction. 

Detailed Description Text (32) : 

When formation crude has a viscosity greater than OBM filtrate, the polarization 
correction used to match CMR porosity to nuclear porosity below the oi 1 -water 
contact will be too large. The overcorrection will make CMR porosity larger than 
nuclear porosity. This is observed in streaks at x760 and x750, where polarization 
corrected CMR porosity (dashed line) is higher than nuclear porosity (heavy solid 
line) . The presence of heavy ni 1 is sometimes associated with the original oi 1 -water 
contact in North Sea reservoirs, and is confirmed by the operator in this example. 

Detailed Description Text (33) : 

Above the paleo OWC, from x740 to x705, is a depleted zone in which formation oil 
has been stripped of gas and light ends. The dead ail here, mixed with mud filtrate, 
appears to have magnetic resonance properties similar to those of the unmixed 
filtrate below x760. The polarization correction using T.sub.l =2.5 sec again 
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perfectly matches CMR^ffrosity to nuclear and core porosities. 



Dei-ailed DpRrr ipf -inn Tpxh (34) : 

CMR porosity departs from the other porosity measurements above x705 feet. This is 
interpreted as the depth at which there is the first appearance of live oil in the 
CMR volume of investigation. This is the native formation n-i 1 , composed of short 
chain alkanes, and having a significant amount of dissolved gas. Its viscosity at 
reservoir conditions is 0.29 cP. The low viscosity of live formation crude means its 
rel axatinn time is significantly longer than that of dead crude and OBM filtrate. 
When live native ni 1 is present, T.sub.l =2.5 sec is no longer large enough to 
provide an adequate polarization correction, and CMR porosity is too lxiw. 

Detail fid Dfianript.ion Text (35) : 

Resistivity shows a very long transition zone spanning more than 80 feet, a 
consequence of the rock fabric: small pore throats trap formation water a t h igh 
capillary pressure. However, because there is no capillary pressure difference 
between OBM filtrate, dead_oJLL^ and live oi 1 , the current f "hii d contact detected by 
CMR at x705 is sharp. Volumetric analysis is shown in Track l. Total hydrocarbon 
saturation is determined by deep resistivity. The live oil volume is found from CMR 
data; the method by which it is determined being set forth in the Appendix below. 

Detailed Description Text (36) : 

Neutron -density crossover indicates gas above x530. Native gas has mixed with the 
dead OBM filtrate, producing a mixture with a reduced viscosity and therefore long 
T.sub.l. This effect can be helpful in finding gas-oi 1 contacts, but hydrogen index 
effects and variability in the degree of mixing make it difficult to quantitatively 
estimate volumes. 

npfailprf Description Text (37) : 

Permeability was estimated by ##EQU11## where .phi. is the Equation (12) porosity 
and BFV is the CMR bouncL_£luid volume. The first term is the Equation (9) 
permeability. The second term is a permeability function appropriate for clean sands 
[see Herron, M. M. , 1987, "Estimating the Intrinsic Permeability of Clastic 
Sediments from Geochemical Data", SPWLA 28th Annual Logging Symposium, Paper HH] . 
The permeability was set to 0.01 mD whenever the permeability flag of Equation (10) 
exceeded 4 with . epsilon. =0 . 25 p.u. Agreement with core permeability is good, within 
the limitations of depth matching and resolution (FIG. 10, Track 2) . 

Detailed npffpn'ptinTi Text (39) : 

Reduced wait time MR logging was used in a second field on the UK continental shelf 
to enhance logging speed and measurement precision. Oil base mud was used to drill 
into a shaly sand light ail reservoir. Conventional porosity tools were used for 
porosity and lithology determination, as above. Porosity was computed as for Log 
Example #1, using a f ln-id density of 0.8 g/cc, appropriate for the oil base mud 
filtrate in this well. CMR logging obtained measurements of bound water and 
permeability. The formation of interest was logged with two passes of the CMR. 
Logging and pulse sequence parameters were (with S=PAP's stacked per 6 inch depth 
frame) 

Detail p:d Description Text. (40) : 

The first pass was designed to obtain a full suite of MR measurements, including 
porosity and free f TiHd porosity. The second pass was designed to determine the 
quality of CMR logs acquired at higher speed, where measurement precision was an 
important consideration. Porosity and t . suh - 2 distribution were derived from the CMR 
tool. 

Detailed TteflrHphinn Text (41) : 

For the slow pass, the standard polarization correction was applied using T.sub.l 
/T . mih . 7. =3, which is typical for formations invaded with oil base mud filtrates. 
For the fast pass, less than 0.1 sec of data are collected in each CPMG acquisition. 
Therefore, the longer t . sub . 2 components in the T. sub. 2 distribution cannot be 
accurately estimated. This leads to errors when the conventional polarization 
correction method, based on the T.sub.l /T. nuh. 2 ratio, is applied. Therefore 
T.sub.l =0.5 sec was assumed for the mud filtrate, and a fixed T.sub.l correction 
applied to the free fluid part of the measured porosity. T.sub.l of the bound fiiiid 
is much shorter, and no polarization correction is required. 

Detail ed DPRrript-inn Text (42) : 

The hydrogen i ndex of the OBM filtrate was determined by obtaining a sample of the 
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base ni 1 used to make np the whole mud. In the shop, th^^MR was used to compare the 
ampi -i i-nrte of the signal from the base oil to an identical volume of the aqueous 
master calibration solution. The hydrogen index was found to be 0.83, a surprisingly 
i nw valiip . This correction was applied to the free fluid component; the hound fin in 
is formation wat.f>r r with HI=1. It was assumed that the neutron tool, which has a 
significantly deeper depth of investigation, was substantially unaffected by the low 
HI of the filtrate. The measured mud filtrate density, used for bulk density 
computation of porosity, was 0.8 g/cc. The logs are shown in FIG. 11. 

pet-ailed D e ^c ri pt i on Text (43) : 

In the slow pass analysis, a wait time of 4 sec was employed to substantially 
polarize the oi l base mud filtrate which partially displaced the connate free fluid 
in the flushed zone. 1800 echoes were collected to define the long T . sub . 2 
components associated with OBM filtrates. One phase alternated pair (PAP) was 
collected in each 6 inch depth frame. 

Detailed Descr iption Text: (44) : 

Total CMR porosity matches density porosity. The T. mih. 2 distribution in the shales 
(Track 3) is generally well resolved and is centered around 3 msec. Note several 
permeable streaks, for example at x734 meters, indicated by high amplitude at long 

T.mih.2 . 

Ttet.a i 1 prf DpRrri pt i on Text (45) : 

In the sand, x678~x705 meters, CMR shows that reservoir quality varies with depth. 
Bound fluid volume (BFV) , the amplitude, below t. mih. a-nnt-.of £ =33 msec, correlates 
with gamma ray, which has been plotted on a scale of 30 to 110 API units to 
emphasize the similarity of the curves (Track 1) . The correlation coefficient (from 
the MatLab function "corrcoef") is 0.80 over the entire section shown. 

Detailed Description Text. (46) : 

In the fast pass analysis the goals were to achieve both better precision in bound 
f Iniri volume and enhanced logging speed. A 0.25 sec wait time and acquisition of 
only 300 echoes permitted a logging speed of 600 ft/hr, with four phase alternated 
pairs stacked per 6 inch depth frame. Thus, the effective signal to noise ratio was 
doubled while the logging speed was improved by a factor of three, compared to the 
conventional pass. 

DptailRr i DPRrription Text (47) : 

In the shales, the fast pass results correlate well with those of the slow pass,- 
free fluid near zero, CMR total porosity in excellent agreement with density 
porosity, and good correlation between bound fluid -volume and gamma ray. The T . fluh . ?. 
distributions for the two runs are similar below 33 msec; the fast pass results 
(Track 4) are more stable due to superior signal to noise ratio. 

npt-a-ilp d Defirription Text (48) : 

In the sands, the fast pass suppressed long rel axation time components and moved 
them to somewhat short er_X*suli J ^2_^ but the use of a polarization correction (T.sub.l 
=0.5 sec) restored free and total CMR porosity (Tracks 2 and 4). The correction is 
large and distorts the t anhi *> distribution, but excellent agreement between fast 
and slow pass porosity demonstrates that the polarization correction is accurate, 
and that no long -T. sub. 1 components are unaccounted for. 

Detail pri Dpsnript-.inn Text (49) : 

As expected, the fast pass improved the definition of the bound fluid in the T.Finh,2 
distributions. The improved signal to noise ratio of the fast pass manifests itself 
in visibly improved correlation between BFV and gamma ray. This conclusion is 
supported by statistical examination. The correlation coefficient is 0.91. 

DPtailf^H Dppirriptinn Text (53) : 

Residual oi 1 saturation in Well #1 was determined by comparing the partially 
polarization-corrected CMR porosity to total porosity determined from density and 
neutron tools. The wait time W=0.45 sec. At the CMR depth of investigation, the 
wat.pr leg is saturated with irreducible water, oil base mud filtrate, and dead 
residual oil. The irreducible water has very short T.sub.l, and no polarization 
correction was applied to it. The T.sub.l 's of the OBM filtrate and dead crude were 
assumed equal and were determined by applying the polarization correction, Equation 
(6) , to the free fluid to match CMR porosity to nuclear porosity in the wat.pr leg; a 
val hp of T.sub.ldead =2.5 sec was found. 
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Dprailpd Descr ipti nn r ^gt (54) : 

The n-i 1 leg contains irreducible nil base mud filtrate, and live native 

hydrocarbon. This mixture, which is rich in gas and has a viscosity at reservoir 
conditions of .eta. =0.2 9 cP, Longitudinal rplayatinn time is estimated by ##EQU13## 
giving approximately T. sub. Hive =6 sec at reservoir temperature. 

Detailed Description Text (55) : 

The uncorrected CMR signal is ##EQU14## For wait time W=0.45 sec, and T.sub.l for 
mud filtrate and <">i 1 as given in the table, 

Dptailpd npfif!ript-.ir>n Text (56) : 

The polarization corrected CMR signal is ##EQU15## where FFI.sup.(-) is the apparent 
free f 1 in d volume, which is less than its true valnp due to insufficient wait time. 
Assuming the bound fluid volume and the water volume are the same, the free f 1 uid is 
corrected by 

D-Ptailpd ne.qrri pti nn Text (60) : 

Solving (A-8) for the volume of live oil 

npi-ailprl Dparrri ptinn Paragraph Equation (1) : 

k=a.phi. .sup. 4 T. sup. 2. sub. 2, log (7) ##EQU6## The constants a and a' are formation 
dependent, with default valuta a=4 mD/msec . sup . 2 and a' =10. sup. 4 mD. The use of 
Equation (7) calls for a porosity measurement and computation of the logarithmic 
mean of the entire t. «nh. 2 distribution, while the use of Equation (8) calls for 
measurements of porosity .phi., bound fin-id volume BFV, and free fluid volume FFI. 
Thus, both relationships ostensibly require the measurement of the full nuclear 
magnpfir rpannanrp t snh 9 distribution, necessitating a long measurement cycle 
time . 

Detailed Tiff arr-ipt inn Paragraph Table (3) : 

Pi iH d Volume T.sub.l (sec) 

Irreducible watex V.sub.W <<0.45 Filtrate + 
Dead Residual V. sub. dead T.sub.ldead =2.5 Live Oil V. sub. live T. sub. Hive = 6 
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Sezginer, A., Kleinberg, R.L., Fukuhara, M. , Latour, L.L., 1991, "Very Rapid 
Simultaneuous Measurement of Nuclear Magnetic Resonance, Spin-Lattice Rel axation Time 
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Nov. -Dec. 1996, p. 61-69. 

Other Reference Publication (10) : 
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White, J., 1997, "Applications of Downhole Magnet -i n Rpannanrp Imaging in the North 
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Symposium, Paper QQ. 
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CLAIMS : 

1. A method for determining residual <">i 1 saturation in formations surrounding an 
earth borehole, comprising the steps of: 

providing a nuclear magnetic resonance logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnetic resonance, spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
being short enough to incompletely polarize the protons in an investigation region 
of the formations ,- 

processing the received spin echoes to obtain a magnetic resonance porosity; 
deriving a measure of total porosity of said formations; and 

determining residual oi 1 saturation from said magneti c resonance porosity and said 
total porosity. 

4. The method as defined by claim 1, further comprising applying a polarization 
correction to said magnetic resonance porosity from said T . sub . 2 distribution. 

5. The method as defined by claim 1, wherein said step of processing the received 
spin echoes to obtain a magneti c resonance porosity includes processing the received 
spin echoes to produce a t snb o distribution, and computing said magnetic resonance 
porosity from said t snb 7. distribution. 

6. The method as defined by claim 5, further comprising applying a polarization 
correction to the free fluid components of said T . nnh ■ 2 distribution before 
computing said ma gnetic resonance porosity from said T.fiuh.2 distribution. 
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7. The method as defined by claim 6, wherein said polarization correction is of the 
form ##EQU17## where P(T.sub.2i) is the corrected T. sub. 2 distribution, 
P{T.sub.2i) .sub.uncorr is the uncorrected t. snh. 2 distribution, and T.sub.l is a 
selected longitudinal rplavatinn time of the free fluid. 

8. The method as defined by claim 5, further comprising the steps of developing a 
free/bound fluid mtnff of Raid T.Hiih.2 distribution and applying a polarization 
correction to the free fluid components of said t. sub 7. distribution before 
computing said magnpfir rpsnnanrp porosity from said T.suh.2 distribution, and 
wherein said step of determining residual oil saturation from said magnet i c 
rpsnnanrp porosity and said total porosity comprises comparing said respective 
porosities in a water -bearing zone of said formations and comparing said respective 
porosities in an ail-bearing zone of said formations. 

9. A method for determining an indication that a formation surrounding an earth 
borehole has a potentially overestimated permeability, comprising the steps of: 

providing a nuclear magnet if rpRnnanpp. logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnpt-i r rpsnnanrp. spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
being short enough to incompletely polarize the protons in an investigation region 
of the formations ; 

processing the received spin echoes and deriving a bound fluid volume, BFV, and a 
free fluid volume FFI.sup.(-) from said received spin echoes; 

deriving a measure of total porosity .phi. of said formations; and 

deriving a flag indicative of potentially overestimated permeability from said BRV, 
FFI . sup . { - ) , and . phi . . 

13 . Apparatus for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole, comprising: a nuclear magnetic resonance logging 
device that is moveable through the borehole ; 

means, in said logging device, for providing a polarizing magnetic field and cycles 
of a magnetic pulse sequence, and for receiving magnetic resonance spin echoes from - 
the formations; 

means for processing the received spin echoes to obtain a magnetic resonance, 
porosity; 

means for deriving a measure of total porosity of said formations; and 

means for comparing said magnpi-ir rpRnnanrp porosity with said total porosity, said 
comparing means comprising means for comparing said respective porosities in a 
watex- bearing zone of said formations and for comparing said respective porosities 
in an oil-bearing zone of said formations. 

15. Apparatus as defined by 14, further comprising means for determining the 
residual ni 1 saturation of said formations from the output of said comparing means. 

16 . Apparatus for determining residual ni 1 saturation in formations surrounding an 
earth borehole, comprising: 

a nuclear magnetic rpsonanrp logging device that is moveable through the borehole; 

means, in from said logging device, for providing a polarizing magnetic field and 
cycles of a magnetic pulse sequence, and for receiving magnetic resonance spin 
echoes from the formations, the polarization wait time, W, of said cycles of 
magnetic pulse sequence being short enough to incompletely polarize the protons in 
an investigation region of the formations; 

means for processing the received spin echoes to obtain a magnetic resonance 
porosity ; 
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means for deriving a measure of total porosity of said formations; and 

means for determining residual ni 1 saturation from said magnpfir rpsn nanrp porosity 
and said total porosity. 

17. A method for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole, comprising the steps of: 

providing a nuclear magnetic resonance logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnetic resonance spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
being short enough to incompletely polarize the protons in an investigation region 
of the formations ,- 

processing the received spin echoes to obtain a magnetic rfisonanne porosity; 

deriving a measure of total porosity of said formations from a density and/or 
neutron porosity measurement in said investigation region of the formations; and 

comparing said magnetic resonance porosity with said total porosity. 

19. A method for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole, comprising the steps of: 

providing a nuclear magnet i n rpsnnanrp logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnpt-i r rpsnnanrp spin echoes from the 
formations ; 

processing the received spin echoes to obtain a magnetic resonance porosity; 
deriving a measure of total porosity of said formations; and 

comparing said magnetic rpsnnanrp. porosity with said total porosity, said comparing 
step comprising comparing said -respective porosities in a_water- bearing zone of said 
formations and comparing said respective porosities in an oil-bearing zone of said 
formations . 

22. A method for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole, comprising the steps of: 

providing a nuclear magnpti r TPsnnanrp logging device that is moveable through the 
borehole; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnetic resonance spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
being short enough to incompletely polarize the protons in an investigation region 
of the formations; 

processing the received spin echoes to obtain a magnetic rpsn nanrp porosity by 

processing the received spin echoes to produce a t snh distribution, applying a > 

polarization correction to the free flni d components of said t mih ? distribution, 

said polarization correction being of the form ##EQU19## where P(T.sub.2i) is the 

corrected t_ mih. 2 distribution, P(T. sub. 2i) . sub.uncorr is the uncorrected t. sub s> 

distribution, and T.sub.l is a selected longitudinal r^T ayat-i on time of the free 

fluid, and computing said magnetic resonance porosity from said polarization 

corrected t. snh. ? distribution; 

deriving a measure of total porosity of said formations; and 
comparing said magnetic resonance porosity with said total porosity. 



15 of 17 



04/14/2003 2:48 PI* 



23. A method for iden^ffying characteristics of hydrocaWons in formations 
surrounding an earth borehole, comprising the steps of: 



providing a nuclear magnptir rpsnnanrp logging device that is moveable through the 
borehole; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnetic rpsonanrp spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
being short enough to incompletely polarize the protons in an investigation region 
of the formations; 

processing the received spin echoes to obtain a magnetic resonance porosity by 
processing the received spin echoes to produce a T. auh. 2 distribution, developing a 
free/bound fluid ctitnff of said t distribution, applying a polarization 

correction to the free fluid components of said T . sub . 2 distribution, and computing 
said magnetic: rp.qnnance porosity from said T . Riib. 2 distribution; 

deriving a measure of total porosity of said formations; and 

comparing said magnet ir rpsnnanrp porosity with said total porosity, said comparing 
step comprising comparing said respective porosities in a water-bearing zone of said 
formations and comparing said respective porosities in an oil-bearing zone of said 
formations . 

24. The method as defined by claim 23, wherein said earth borehole was drilled with 
ni 1 base mud, and wherein said step of comparing said magnetic resonance porosity 
with said total porosity comprises comparing said respective porosities in a 
watex-bearing zone of said formations that has been flushed with oil base mud 
filtrate and comparing said respective porosities in an_oil- bearing zone of said 
formations . 

25. The method as defined by claim 24, further comprising the step of identifying 
the presence of light ni 1 in said formations when the mrrflnhfi d magnetic r pRnnanrp 
porosity is less than the total porosity. 

26. The method as defined by claim 24, further comprising the step of identifying 
the presence of heavy oil in said formations when the corrected magnetic resonance 
porosity is greater than the total porosity. 

27. The method as defined by claim 23, -further comprising the step of identifying 
the presence of light ni 1 in said formations when the rnrrprfpd magnptir resnnannp 
porosity is less than the total porosity. 

28. The method as defined by claim 23, further comprising the step of identifying 

the presence Of heavy ni 1 in said formations when the rnrrprhpd magnpfir rpflnnanrp 

porosity is greater than the total porosity. 

29. A method for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole that has been drilled with oil base mud, comprising 
the steps of: 

providing a nuclear magneti c rpsnnanpp logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnet \r. rpsonanrp spin echoes from the 
formations; 

processing the received spin echoes to obtain a magnetic resonance 
porosity; 

deriving a measure of total porosity of said formations; and 

comparing said magnetic; rpRnnanpp porosity with said total porosity. 

31. The method as defined by claim 29, further comprising the step of computing the 
residual r>i ~\ saturation of said formations as a function of said comparison. 



32. A method for identifying characteristics of hydrocarbons in formations 
surrounding an earth borehole, comprising the steps of: 

providing a nuclear m^gn^Mr- rpsnnanrp logging device that is moveable through the 
borehole ; 

providing, from said logging device, a polarizing magnetic field and cycles of a 
magnetic pulse sequence, and receiving magnptir rpsnnanrp spin echoes from the 
formations, the polarization wait time, W, of said cycles of magnetic pulse sequence 
is short enough to incompletely polarize the protons in an investigation region of 
the formations; 

processing the received spin echoes to obtain a magnetic resonance porosity; 
deriving a measure of total porosity of said formations; and 

comparing said magnpt-i c rpsnnanrp porosity with said total porosity and computing 
the residual ni 1 saturation of said formations as a function of said comparison. 
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Ahstract Text (1) : 

The present invention is a method to rapidly determine the fluid-flow permeability 

of porous media with nnrlpa r magnpti r. rpsnnanrp (HMR) . The method can be applied to 

measurements of permeability in fluid- saturated earth formations using NMR logging 
tools. 

Tnvpntnr City (1) : 
High Bridge 

Tnvpnfnr flronp (1) : 

Jerosch-Herold; Michael High Bridge NJ 

Rrifif fin mmary T^yt {2) : 

The present invention relates to a method for determining the fluid flow 
permeability of porous media. In particular, the present invention relates to 

measuring the permeability by nuclear magneti c rpsnnanrp (NMR) using the T , sub ■ 2 

trI ayat-inn time for the decay of the transverse magnetization of f 1 uiris saturating 
the medium under investigation. 

Brief Summary Text (3) : 

The properties of flni d in porous media are of great relevance in many fields of 
science and engineering. There are numerous measurements which bear some importance 
on characterizing fin-id properties in confined geometries like sandstone rocks. We 
list here a few: porosity, fluid flow permeability (both dc and ac) , electrical 
conductivity, wettability, etc. Quantities like porosity and fluid flow permeability 
in porous rocks are of great relevance for determining, the producibility of 
petroleum reservoirs. 

Brief Summary Text (4) : 

For porous media it has become customary to speak of the solid material which forms 
the "backbone" as the matrix and its complement as the pore space. Porosity is 
defined as the ratio of pore space volume inside the porous material to the total 
volume of the porous medium. Permeability is a measure for the ability of porous 
materials like e.g. porous rocks to permit fluid flow through the pore space. It 
generally increases with porosity, but also depends on other parameters of the rocks 
as e.g. the specific surface area of the pore space, the pore size distribution and 
the pore shape. The fluid flow permeability can vary by about 8 orders of magnitude 
in loose sediments and sedimentary rocks. It has the dimension of area and is 
defined by Darcy's law which relates the rate of fluid flow to the pressure 
differential between two parallel planes for inflow and outflow. The f Inid flow 
permeability is measured in the laboratory by fitting sleeves to core samples which 
are often cylindrically shaped. The top and bottom of the core samples are connected 
to fluid inlets and outlets and a known pressure difference is applied across the 
sample. The f Tin d flow rate is measured for a set of different pressure gradient. 
Liquids or gases can be used as flowing medium, although the measurement using a 
liquid is generally easier as in most cases the liquid can be considered 
incompressible. The laboratory procedure therefore requires first to drill core 
plugs from core samples, which have to be cleaned with various solvents. In contrast 
the method of the present invention can be carried out with a nuclear magnet i r. 
rPRnnanrp logging tool to measure in situ the tranflwrHP relaxation time of the 
f lnidg saturating an earth formation to accurately predict the fluid flow 
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permeability of the earth formation. 



R-r-jpf Summary Tpyt (5) : 

Nuclear magnet- i < - rpsnnanrp (tmR) has been employed for some time to study f lindn 
permeating the pore space of porous media [see J. R. Banavar and L. M. Schwartz, 
"Molecular Dynamics in Restricted Geometries", chapter 10, edited by J. Klafter and 
J. M. Drake, J. Wiley (1989)] . The fl nid supplies the probe particles which diffuse 
in the pore space. Since the classic paper by Brownstein and Tarr (BT) [see K. R. 
Brownstein and C. E. Tarr, Physical Review A, 19, 2446(1979)] it has been realized 
that nuclear spin rplayati nn can provide information about the pore space geometry. 
BT discussed the case of T.sub.l and T.fluh.2 relaxation in an isolated pore where 
the nuclear spins are rel aypd by collisions with the pore walls. The interpretation 
of T.sub.l measurements with this model for fluids in porous media can present 
several problems. In the limit where the nuclear spins diffuse at a fast rate to the 
pore surface and the surface relaxation is in comparison relatively slow, the 
avpragpd relaxation curve can be related to the pore size probability distribution. 
In this so called fast diffusion limit where the lowest order relaxation mode 
dominates one still has to assume that the surface relaxation strength is uniform 
and the pores are isolated to relate the distribution of relaxation times uniquely 
to the pore size distribution. It is conceivable to have porous samples with the 
same pore size geometry but different levels of paramagnetic impurities which 
influence the mirfarp relaxation velocity while the f lui d flow permeability would 
remain unchanged. To obtain a reliable estimate of the fluid flow permeability with 
NMR one therefore has to perform an experiment which directly probes fluid transport 
in the porous medium like for example the diffusion of fln-i d molecules in the pore 
space. For T.sub.l measurements the nuclear spin rp-T arafinn depends on the rate at 
which magnetization is carried to the surface but also on the surface relaxation 
velocity .rho. . As the surface re»1 axation strength .rho. has no bearing on 
permeability one can therefore hope to correlate T.sub.l and the fluid flow 
permeability only for classes of materials with similar «nrfarp> rpl axati on 
properties. 

Brief Summary Text (6) : 

There is an increasing interest in applying NMR in well-bore environments to 
determine the properties of fluid carrying earth formations [see P. N. Sen, C. 
Straley, W. E. Kenyon and M. S. Whittingham, Geophysics, 55, 61-69(1990)]. This 
interest has been spurred by the introduction of a new generation of NMR logging 
tools by NUMAR [see M. N. Miller, A. Paltiel, M. E. Gillen, J. Granot and J. C. 
Brouton, Society of Petroleum Engineers, SPE 20561, 321(1990)], which are already 
being used in the field. The new wmr logging tools are very well fitted to carry out 
the physical* measurements required for our method of invention. 

Brief Summary Text (7) : 

In the present invention, a measurement of the f. ran averse re laxa t i on time T . .sub . 2 
for fluids in porous media is used to determine the permeability of the medium by 
taking advantage of magnetic field inhomogeneities across pores. For strong magnetic 
fields and in the fast diffusion limit the relaxation is determined to first order 
by the transport of magnetization through the pore space and not the surface 
rplavaUnn velocity. It will be shown that it is possible to correlate T . sub . 2 to a 
length characteristic of the pore space geometry which can also be determined 
independently from mercury injection experiments and thereby relate T.mih.2 tn the 
fluid flow permeability. It is also feasible to study the degree to which the 
diffusion of fluid molecules is restricted by the pore space geometry. T.mih.2 for 
f ini da in porous media is in general orders of magnitude shorter than T.sub.l in 
marked contrast to the situation for bul k fluids . The main mechanism for T FHih.2 
relaxation of t-he fluid spins in strong magnetic fields is due to the internal 
random magnetic field gradients generated by the difference in magnetic 
susceptibility for the fluid filling the pore space and the material making up the 
matrix of the porous medium. At„_low fields surface relaxation can not be be 
neglected but the .tau. dependence of t . mih . 7. (.tau.) is still primarily due to 
diffusion in the internal magnetic field gradients. Surface rel axation will under 
standard experimental conditions not lead to a .tau. dependence of T. nnh. 2 in a CPMG 
experiment. This is confirmed by recent experimental results at low field with an 
NMR logging tool [see M. N. Miller, A. Paltiel, M. E. Gillen, J. Granot and J. C. 
Brouton, Society of Petroleum Engineers, SPE 20561, 321(1990)]. The spatial 
dependence of the internal gradients is determined by the pore space geometry and 
pore size distribution. The internal gradients in turn determine the rate at which 
the spins diffusing through the pore space loose their phase memory. The loss of 
phase memory can be monitored with a multi spin-echo pulse sequence like the 
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Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence [see S.^eiboom and D. Gill, Rev. 
Sci. Instr., 29, 688 (1958)]. 

Brief Summary Text (8) : 

The phenomenon of spin echoes essential to the present invention was first 
discovered in MME by Erwin Hahn [see E. L. Hahn, Phys. Rev., 77, 297 (1950)]. In an 
inhomogeneous magnetic field nuclear spins will precess at a Larmor frequency, 
.nu..sub.L, determined by the local field. After an initial radiof requency pulse 
which tips the spins into a plane transvprsp to the direction of the applied static 
magnetic field the spins are all in phase and the sum of the total transverse 
magnetization is at the maximum possible scaliie.. Due to the spread in precession 
frequencies the spins will dephase and the macroscopic magnetization measured with 
the NMR instrument will decay. It is useful to remember here that the macroscopic 
magnetization is a vector sum of the magnetic moments of the spins which vanishes 
when the phases of the magnetic moments are random. One can reverse the dephasing 
process by applying a 180 degree pulse a time .tau./2 after the initial 
radio- frequency pulse which tipped the nuclear spins into the transverse plane. 
Immediately after this pulse a spin which precesses at a faster frequency than the 
avsrage lags behind by an angle which is exactly the same angle by which it was 
ahead of the avprage immediately before the 180 degree pulse. Similarly spins 
precessing at a frequency slower than the_ai£exage. are now ahead. A time .tau./2 
after the 180 degree pulse the spins will be again be in phase and one can observe a 
spin-echo. Spins diffusing will be subject to different local fields between the 
time the first pulse was applied and the detection of the spin echo. As their Larmor 
frequency is not constant the refocusing of magnetization will be incomplete and the 
echo will be attenuated. The degree of attenuation depends on the displacement and 
field inhomogeneity. This attenuation can be used to measure diffusion constants in 
■Fin-idfi and to probe the diffusion of fluid spins in the pore space of porous media. 

^Hpf Rurnmary Tpxt (10) : 

The present invention provides a method for determining the permeability of porous 
media saturated with a liquid using nuclear magnetic resonance (NMR) . The steps of 
the method include: (a) applying a radiof requency pulse sequence which after an 
initial pulse generates successive spin echoes with a train of radio frequency 
pulses spaced apart by a time interval of length .tau. wherein all pulses have a 
carrier frequency corresponding to the Larmor frequency of the f 1 in d spins filling 
the pore space of the medium for which the f luid flow permeability is to be 
determined; (b) measuring the decay of the transverse magnetization at each of the 
successive regularly spaced midpoints between the 18 0 degree pulses where the 
midpoints coincide with the peak of the spin echoes; (c) repeating steps a and b at 
least one more time wherein each repeat of step (a) uses said radio frequency pulse 
train with a different valnp of the pulse spacing .tau.; (d) determining the 
rranRvprfip rplayatinn time T . sub . 2 (.tau.), from the transverse magnetization decay 
for each valnp of .tau. and determining one of a pref actor .DELTA, or a restricted 
diffusion length l.sub.nmr from said T . suh. 2 (.tau.); (e) measuring the porosity of 
said porous media; (f ) determining the permeability of said media from the porosity 
and either a pref actor .DELTA, or a restricted diffusion length l.sub.nmr from said 
T. sub. 2 ( .tau. ) . 

Rrief Summary Text (11) : 

In a preferred embodiment the method is performed as a down-hole wellbore 
measurement to swiftly and accurately determine the fluid -producing potential of an 
earth formation using a magnetic; resonance logging system which employs static and 
radio -frequency magnetic fields to perform the spin- echo CPMG pulse experiment in a 
wellbore environment. 

Drawing DPRrripHrtn Teyr (3) : 

FIG. lb shows the numberically evaluated CPMG decay for a certain .tau. value and 
the pore size distribution shown in (la) ; 

Drawing Description Text (4) : 

FIG. lc shows the t. mih. 2 values calculated from computer simulated CPMG decays for 
the same pore size distribution. 

Drawing Description Text (5) : 

FIG. 2 shows a fit of t . aub . 2 (.tau.) for a restricted diffusion model. 
Drawing npsnH pti on Tpxt (7) : 

FIG. 4 shows the valuta o f t . F»ih . 2 ( . tau . . f wdarw . 0 ) obtained from numerical 
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calculations for a set of pore size distributions versu^the width .DELTA, of the 
pore size distribution. 

Drawing Description Text (8) : 

FIG. 5: Schematic diagram of the NME spin-echo pulse sequence used for the 
measurement of t. snh.a - 

Drawing Description Text (9) : 

FIG. 6 shows an example of stretched exponential fit to CPMG decay for a water 
imbibed sandstone sample . 

Drawing Pegf-ri pt i on Text (10) : 

FIG. 7 shows the dependence of T.snh.2 on CPMG pulse spacing .tau. for a set of 
sandstone samples. The t wih ? valnps can be extrapolated to T . .sub . 2 
{ . tau. . fwdarw. 0} with an exponential function. 

Drawing npRcnphinn Text (11) : 

FIG. 8 shows a graph of extrapolated value of T.mib.2 versus l.sub.c determined from 
mercury injection for sandstone rock samples. 

drawing Dfiscriphinn Text (12) : 

FIG. 9 shows a graph of T.sub.2 ( .tau. =0) .times. .phi. (.phi. is the porosity) versus 
permeability k for sandstone rock samples. 

Deta i 1 ed Desrr i pt i on Text (2) : 

The present invention is a method to determine the permeability of porous media 
using pulsed nnrlpar magnpfir rpsnnanrp . This method uses the relaxation decay of 
the t-ransvprsp. component of the magnetization (T.mih. a rplavat.inn ) measured with the 
Carr-Purcell- Meiboom-Gill (CPMG) pulse sequence. 

neta -iled DP-snr i pt i on Text (4) : 

where .beta.. sub. 2 is the stretch exponent. A nonlinear least squares fit of the 
measured values of magnetization m(t) to the stretched exponential function 
determines th e relaxation time t suh. 7. and the stretch exponent, .beta. . sub. 2 . If 
the pulse spacing, .tau., is changed in incremental steps, a set of stretched 
exponential curves is produced yielding different values for T.fnih.2 (.tau.). For 
the case of free diffusion, the pore size distribution introduces a spent nam of 
t ctnh o decay values because the internal magnetic field gradient is related to the 
pore size. The analysis of the transverse magnetization decay in this case yields an 
effective T.Knb.2 f whirh is a weighted average over the range of T.suh.2 values.. In 
this case a characteristic length can be determined from the pref actor of the - 
functional form that describes the .tau. dependence of T.suh.2 (.tau.). When T . sub . 2 
is only measured for a limited range of .tau. values this pref actor is dominated by 
the pore size distribution and is not very sensitive to the exact functional form of 
t - <=nih , i (.tau.). A characteristic length, .DELTA. , is derived from this pref actor. 

Detailed Description Text (5) : 

where F(.tau.) is a fitting function which reproduces the .tau. dependence of 
t. gun 2 ■ Specific forms for F(.tau.) are discussed in the section on Theoretical 
Background . 

Detailed Description Text (6) : 

For a sufficiently wide range of .tau. values, the exact functional form of T. suh . 2 
(.tau.) is determined by the cross-over from free to restricted diffusion. For a 
given range of .tau. vainer, this cross-over is also dependent on the pore size 
distribution. Analysis of the functional form of T . suh . 2 (.tau.) in the region of 
this cross-over yields a characteristic length, l.sub.nmr, according to: ##EQU1## In 
either case, the characteristic lengths, .DELTA, or l.sub.nmr, derived from the 
decay of t . suh . 2 (.tau.) with .tau., are a measure of the f tnid flow permeability. 
With l.sub.nmr the flni d flow permeability is determined with 

Pet-aile d Descri ption Text (7) : 

If a prefactor .DELTA, has been determined from t. gnh. ? (.tau.) the fluid flow 
permeability is determined with: 

Detailed Description Text (8) : 

If T . sub . 7. (.tau.) is a simple exponential then the prefactor .DELTA, and the 
extrapolated value of T. sub .2 (.tau.) are related by T . sub . 2 (. tau. . fwdarw. o) 
.varies. .DELTA. . The value of T.anb.2 extrapolated to .tau.=0 combined with the 



4 of 10 



04/14/2003 2:48 P^ 



porosity, .phi., of the material gives the permeabilit}^^:, of the porous medium 
according to (see FIG. 9) : 



Del-aileri Deacripti nn Tpxt ( 9) : 

The permeability relations 4, 5 and 6 follow from the relations: k .varies. 

1 . sub . c . sup . 2 . phi . . sup . 2 and . DELTA . . varies . 1 . sub . nmr . varies . 1 . sub . c and T. anh.? 
{ . tau. . fwdarw. 0) .varies . .DELTA. . The proportionality constants in these 
permeability relationships 4, 5 and 6 can be determined from a calibration 
experiment with a porous material of known porosity and permeability. The exponents 
n and m may vary depending on which of the above relations 4, 5 and 6 is used to 
determine the permeability. The vahips of n and m are not necessarily the same for 
each of the relationships. The porosity is determined by first running an nmr 
experiment with a wat.pr jacket which surrounds the radiof requency probe {e.g. of the 
NMR logging tool) and therefore calibrates the NMR signal for 100% porosity. To 
prevent the interference from local radio signals one can enclose the probe and 
wafpr jacket in a Faraday cage. As the NMR signal is directly proportional to the 
number of hydrogen nuclei in the sensitive volume of the probe it will scale 
linearly to 1 nwpr porosities. Only the contribution from liquid phases or 
non-adsorbed fluid will be recorded if the dead time of the instrument is set 
accordingly. The temperature of the fluid in the sensitive volume is of secondary 
influence. In a bore hole environment this varies with depth and can be compensated 
for by measuring the bottom temperature and temperature gradient. The NMR signal 
amplitude changes as a function of temperature as 1/T. 

Detailed Description Text. (11) : 

The magnetization of the fluid spins diffusing in the pore space satisfies the 
following modified Bloch equation: ##EQU2## T.sub.b is the bulk rel axati nn rate of 
the f 1 tii d in the pore space and D is the diffusion constant which is on the order of 

2 . times . 10 . sup. 9 [m.sup.2 /sec] for water at room temperature. The boundary 
condition is: 

Detailed Defirri pti on Text (12) : 

where .rho. is the mirfarft relaxation velocity, which has units of length over time 
and can be thought of as the -relaxation rate at the surface multiplied by the 
thickness of the layer of fln-i d spins rplaxi ng near the surface. The bulk relaxation 
rate can always be factored out of the solution. 

Detailed Description Text (13) : 

The Bloch equation becomes: ##EQU3## Brownstein and Tarr [K. R. Brownstein and C. E. 
Tarr, Physical Review A, 19, 2446-2453 (1979)] expressed the general solution as a 
sum of normal modes: ##EQU4## Two limiting cases can be considered for the 
magnetization decay which are characterized by decay constants .tau. .sub.c : 
##EQU5## In the fast diffusion case the magnetization is approximately uniform 
across an isolated pore and only the 1 nwpsf mode in the above eigenmode expansion 
contributes significantly to the magnetization decay. In this limit, . rho.r/D<<l, 
and when the pores are approximately isolated (narrow throat limit) one obtains for 
m(t) : ##EQU6## S.sub.p and V.sub.p are the surface area and the volume of a pore and 
V.sub.p /S.sub.p .apprxeq.r. For a pore size distribution the -relaxation is 
described by: ##EQU7## 

Detailed Description Text (14) : 

r. sub.c is a lower- r*ut-off on the pore size distribution. Thompson et al [see A. H. 
Thompson, S. W. Sinton, S. L. Huff, A. J. Katz, R. A. Raschke and G. A. Gist, 
Journal of Applied Physics, 65, 3259 (1989)] found that the pore size distribution 
of many porous rocks is well represented by the following class of functions: 

Detailed Description Text (15) : 

and where .DELTA, is a measure of the width of the pore size distribution and .beta, 
lies between 0 and 1. A .beta, value, of 2/3 will yield a Guassian pore size 
distribution. This type of pore size distribution will lead to a stretched 
exponential decay in the fast diffusion limit: 

Detailed Description Text (17) : 

Equation 16 is of a form observed experimentally for T.sub.l and T. nub. 2 
magnetization decays for f 1 in ds in porous media for a broad set of experimental 
conditions, although it was derived for the case of surface induced rel aval- inn we 
wish to arrive at an expression which explicity accounts for the contribution to the 
t siih 7 decay from diffusion of fluid spins in porous media in the presence of 
internal magnetic field gradients. The magnetic field gradients lead to dephasing of 
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the nuclear spins which can only be partially compensator with a spin- echo sequence 
as they diffuse. As a starting point we use an expression derived by Robertson and 
Neuman and experimentally confirmed by Wayne and Cotts for the decay of the 
hranHVPrgp magnetization in the presence of a uniform gradient G [see B. Robertson 
Physical Review, 151(1), 273 (1966), R. C. Wayne and R. M. Cotts, Physical Review, 
151(1), 263 (1966), C. H. Neuman, The Journal of Chemical Physics, 60 (11), 4508 
(1974)] . Neuman produced expressions for the case where the spins are assumed to 
diffuse in a bounded medium of spherical geometry and the Carr-Purcell spin-echo 
sequence with spacing of the .pi. pulses given by .tau. is being used: ##EQU8## 
. theta.=2.tau.D/r.sup.2 and the .alpha, are determined from tan . alpha .. sub . i 
=2 . alpha. . sub. i / (2- .alpha. .sub. i. sup. 2) whose solutions asymptotically approach: 
.alpha. .sub. i =i .multidot . .pi . . For small .pi. the above expression agrees with the 
well-known expression for transverse rel ayat-.inn of spins due to unrestricted 
diffusion in a uniform gradient: 

netaiTed npscript-ion Text (19) : 

. DELTA .. sub .. chi . is the susceptibility difference, .mu..sub.o the magnetic 
permeability of vacuum and H.sub.o the magnetic field strength. This means that 
under this model the gradient is uniform over individual pores but varies from pore 
to pore. The contribution of the spins in each pore has to be weighted by the volume 
of the pore. For a pore size distribution we must average, the magnetization decay, 
m(t) , over the pore size distribution in a fashion similar to the example in 
equation 14. We therefore arrive at the following expression for the contribution to 
the t-ranavfirsg magnetization of the spins diffusing in random internal gradient 
fields: ##EQ09## where the last exponential explicitly accounts for the possibility 
of surface relaxation. This expression is the basis for numerical calculations 
performed for a series of pore size distributions. This expression was also 
previously, independently derived by Kleinberg and Horsfield [see R. L. Kleinberg 
and M. A. Horsfield, Journal nf Magnptir Re sonance, 88, 9-19 (1990)]. We numerically 
calculated the CPMG echo decay curve with the sum over the roots .varies. .sub. i 
being approximated up to the 20th term. The pore size distribution is of the form 
shown in Equation 15 . The simulations were carried out in the same manner as the 
experiment. The CPMG decays were calculated for a set of valuta Q f .tau.. We observe 
that the calculated magnetization decays are well described by stretched exponential 
functions. This suggests that equation 16 applies to a broader class nf relaxation 
decays than suggested by the derivation which only considered surface induced 
relaxation. Indeed the numerical results show that the decay of the magnetization 
due to diffusion in the presence of random magnetic field gradients and averaged 
over a pore size distribution also leads to stretched exponential decays. 

Detailed Description Text (20) : 

FIG. la shows the pore size distribution calculated from equation 15 for typical 
values of .DELTA, and .beta. . FIG. lb shows the magnetization decay calculated for 
the pore size distribution of la from equation 20. This calculation was repeated for 
several_3£alues of .tau. and the nonlinear least squares fits to the stretched 
exponential function yield values fr>r t Bnh.2 and .beta . . sub. 2 . FIG. lc shows the 
t giih> • ? values thus obtained for a set of .tau. values 



De^-iled Desrrription Text (21) : 

Restricted diffusion can be expected at longer valnp.q of .tau.. The diffusion 
dynamics were therefore studied by systematically varying the parmameter .tau. in 
the CPMG sequence. The crossover from free to restricted diffusion will then be 
evident in the dependence of both t. wih. ? and .beta. .sub. 2 on the diffusion time set 
by .tau. . CPMG experiments were repeated for a series of .tau. values. Values, for 
t hiiH •? and .beta., sub. 2 were obtained by non-linear least squares fits to a 
stretched exponential function for each coherence decay. For free diffusion is a 
linear gradient G the CPMG T . sub . 2 . sup . - 1 should be . varies . ( . gamma . G) . sup . 2 
D. tau. .sup. 2, where D is the diffusion constant. For restricted diffusion D becomes 
effectively time dependent. To explain the .tau. dependence of t mih o. we use a 
model which interpolates between the limits of free diffusion considered above and 
restricted diffusion at long .tau. . The cross-over from unrestricted to restricted 
diffusion is well characterized by: 

Dp tailed Descri ption Text (22) : 

a result which was derived first in a pioneering paper by Stejskal on NUR 
measurements of restricted diffusion [E. O. Stejskal, Journal of Chemical Physics, 
43, 3597-3603 (1965)]. l.sub.nmr should be on the order of V.sub.p /S.sub.p, the 
volume -to -surface ratio of a pore. This gives a .tau. dependence for T . sun . 2 of the 
form: ##EQU10## This approximation for T . sub . 2 (.tau.) leads to the correct limiting 
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behavior for . tau . . f wdarw . 0 . The experimental data poi^^ f or T. nub. 2 ( . tau. ) . sup. -1 
fit well to the above equation. We observe that for typical parameters. 4D.sub.o 
. tau. /I . sub.nmr . sup. 2 <<1. A power series expansion of equation 22 is therefore 
adequate as an approximation for t Rnh 7 {.tau.): ##EQU11## 

Detaile d Description Text (23) : 

The fit of equation 23 or 22 to the experimental t. <=mb. s> (.tau.) data for several 
sandstone rocks are shown in FIG. 2. The fits to the experimental data yield 
coefficients for the two terms in the expression for T . nub . 2 (.tau.) .sup.-l of 
opposite sign as suggested by equation 23. The resulting values, for 1. sub.nmr are 
seen in FIG. 3 to correlate with a quantity l.sub.c derived from mercury injection, 
l.sub.c provides an estimate of the fluid flow permeability [see A. H. Thompson, S. 
W. Sinton, S. L. Huff, A. J. Katz, R. A. Raschke and G. A. Gist, J. Appl . Phys., 
65(8), 3259 (1989)]. Therefore since l.sub.c is proportional to 1. sub.nmr, 1. sub.nmr 
provides an estimate of the flu-id flow permeability. Mercury injection data indicate 
that the exponents n and m in equation 4 should be on the order of 2 . 0 . 

Detailed npsrripfinn Tevl- (24) : 

Equation 22 is one approximation to the crossover from free to restricted diffusion 
and other expressions have been proposed in the literature. All have in common that 
they provide a characteristic length 1. sub.nmr which is of relevance to determine 
the transport properties of the porous medium. The restricted diffusion analysis 
presented above is therefore not confined to the use of equation 22 but can be used 
for a variety of models describing the crossover from free to restricted diffusion. 
An example of an alternative model is equation 20 (used for numerical calculations) , 
where the sum in the exponent can be used to define an effective diffusion 
coefficient. In each case the crossover from free to restricted diffusion is probed 
by measuring the .tau. dependence of T . sub . 2 (.tau.) as outlined in this invention. 

Detail fid Description Te.xt (26) : 

In actual practice the ability to determine the exact functional form which 
describes the .tau. dependence of t. mih.2 (.tau.) is limited by the time available 
for signal collection. This limits both the range of .tau. values which can be 
measured and the signal-to-noise ratio of the data. However in each case there is a 
characteristic length which determines the dominant contribution to T . nub. 2 (.tau.). 
This length is related to the mean life time <t> for the magnetization decay. The 
mean life time is defined by [see D. J. Wilkinson, D. L. Johnson and L. M. 
Schwartz", Phys. Rev. B, 44, 4960-4973 (1991)]: ##EQU12## For stretched exponential 
decays the mean life time and t anh. 7. from a nonlinear least squares fit to a 
stretched exponential are simply related: <t->=r «uh. 7 .GAMMA. (1/ .beta. . sub. 2) and 
for the range of . beta sub . 2_scalii£S typical for magnetization decays- in rocks the 
dependence of the mean life time on .beta.. sub. 2 is weak. Using the approximation of 
free diffusion and with the assumption that the magnetic field gradient is inversely 
proportional to the pore radius the mean life time for the t-ranRvprsp magnetization 
decay is given by: ##EQU13## where .delta. .omega. = .gamma. Gr 

.apprxeq. .gamma. .DELTA. .chi.H. sub. o. The free diffusion expression is a good 
approximation when t. is measured only for a limited set of short .tau. valuer. 

Using equation 15 for the pore size distribution and a change of variable 
( (r/. DELTA. ) .sup. .beta./ (1-. beta. ) =x) , the mean lifetime is given by: ##EQU14## 
Significantly the mean life time and therefore also T . mih . 2 are directly 
proportional to .DELTA, and only weakly dependent on .beta, and r.sub.c /.DELTA, for 
typical valnpfi of these two parameters. This relation between T. sub. 2 and .DELTA, is 
a consequence of the pore size distribution. This result remains unchanged if we 
generalize the expression for the mean life time to include a different .tau. 
dependence for t ptnh. 7 (.tau.). This is evident by substituting a general functional 
dependence t. snh.a .varies. F(.tau.) in place of the t snh 7. .varies, .tau. -.sup. 2 
in the above integral. ##EQU15## 

Detailed DeRrrriptinn Text (27) : 

Empirically we have observed that at a magnetic field of 7T and for the range of 
experimental .tau. values used, the exponential function is a good approximation for 
F(.tau.). At other magnetic field strengths the simple exponential may not 
necessarily be a good approximation for F(.tau.. For the numerical simulations, we 
repeated the calculation of the T . rsuh . 2 (.tau.) data set for different pore size 
distributions where we arbitrarily changed .DELTA, and .beta, to change the shape of 
the pore size distribution. An exponential dependence of T.mih.2 (.tau.) is seen to 
be also an excellent approximation for the numerically calculated results for 
T . Hub . 2 (.tau.) as is evident from the example of FIG. lc. For this case of a simple 
exponential dependence of T . sub . 2 (.tau.), the pref actor of the exponential is 



7 of 10 04/14/2003 2:48 P^ 



f 




equivalent to the valnp o f t snh 2 extrapolated to . tau^u. Thus the exponential 
extrapolation of t mih. 2 to .tau.=0 is a measure of the prefactor of the 
exponential. We therefore expect that the extrapolated value o f Tsnb.2 to .tau.= 0 
is proportional to .DELTA, and this is confirmed by numerical simulations as shown 
in FIG. 4. 

Detailed Dpsrriptinn Text (28) : 

The procedure of estimating . DELTA, from the .tau. dependence of T.anh. 2 either by 
prefactor analysis or extrapolation provides a simple estimate of the fl n-i d flow 
permeability: 

Detailed Description Text (29) : 

The proportionality r>f t snh 2 ( . tau. . fwdarw. 0) to .DELTA, (and therefore l.sub.c) 
allows one to use T.sub.2 (. tau .. fwdarw. 0) to predict the fluid flow permeability 
from NMR measurements: 

Detailed Description Text (33) : 

The method of this invention was tested in the laboratory using a nuclear magnet i c. 
rpsnnanpp spectrometer whose functionality can replicate the capabilities of an NMR 
logging tool instrument. The samples used were sandstone core plugs from various 
geological formations in North America. For all porous samples we used the following 
procedure to prepare the samples: The samples were placed in a sealed container and 
imbibed with Mater. After imbibing a core plug with water for several hours it was 
taken out, sealed with Teflon tape and transferred to an wmp glass tube for 
measurements. In between measurements the samples were kept under_wat£r. in a sealed 
container. 

neta-iTed Descr -ipt ion Text (34) : 

Several samples of sandstone rock were used in developing the present invention. The 
permeability and porosity of these samples was determined by standard methods. In 
addition on samples obtained from the same batch of rock cores mercury injection 
experiments were performed to determine a characteristic pore (throat) size l.sub.c. 
Mercury is a non-wetting fluid and under an applied pressure the mercury will first 
penetrate the largest pores of the medium. For a certain threshold pressure the 
first continuous path of mercury will be formed between the two ends of the sample 
and this will be detectable as a jump in the electrical conductivity across the 
sample. The length l.sub.c can be calculated with the Washburn equation from the 
threshold pressure at which this first conducting mercury path is established across 
the sample. The valnes for permeability, porosity and l.sub.c for the set of 
sandstone samples used for the nmr experiments are listed in Table 1. These values 
of permeability and l.sub.c obtained independently from the NMR measurements will 
have to be compared with the predictions using T ■ sub ■ 2 (. tau. . fwdarw. 0) . 

Detailed Desoript inn Text (35) : 

All laboratory NME experiments were performed at a field strength of 7.05 Tesla 
corresponding to a Larmor frequency of the hydrogen nucleus of .nu..sub.L =300.13 
MHz . The T.sub.l measurements were made with a standard inversion recovery pulse 
sequence. The T.sub.l measurements were only made for completeness and are not 
important for the development of this invention. To acquire t. «nh. 2 data we use the 
Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence which is comprised of an initial 90 
radio frequency pulse which tips the longitudinal magnetization into the transverse 
plane. It is followed a time . tau./2 later by a train of 180 pulses with constant 
pulse spacing .tau.. At the midpoint between the 180 pulses the nuclear spin 
magnetization of stationary spins is refocused and a Hahn echo is formed. The echo 
will be attenuated for spins moving along the orientation of the magnetic field 
gradients. The magnetization is sampled at the center of the Hahn echoes. This means 
that for a train of n 180. degree, pulses we acquire n data points. We repeat this 
experiment for a set of .tau. values - This procedure corresponds exactly to one 
outlined above in the section titled "Method of Invention". FIG. 5 shows a schematic 
diagram of the pulse sequence. 

Detailed Description Text (36) : 

When measuring the t anh 7. decay with the CPMG sequence as a function of pulse 
spacing, t. mih. 2 decreases rapidly for porous samples imbibed with water or oil . By 
using the CPMG pulse sequence we can ascribe this decay of T . sub . 2 to the effects of 
internal magnetic field gradients in the porous sample, which are mainly due to the 
magnetic susceptibility difference for the fin-id imbibing the sample and the matrix. 
For most samples we recorded the decay of the transverse magnetization with the CPMG 
sequence for a set of 8-16 values of .tau. . 
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Detailed Df».qcriptinn Text (38) : 

The T. snh. o. magnetization decays were fit to a stretched exponential function with a 
nonlinear least squares fitting procedure. FIG. 6 shows an example of the stretched 
exponential fit to a t wnh ? CPMG decay curve. The stretched exponential function 
approximates the experimental data well over 2-3 orders of magnitude as can be seen 
from FIG. 6b. No data points were discarded before fitting the data to the stretched 
exponential function. 

Detailed Inscription Text (39) : 

We acquired for each sample a set of CPMG decay curves by varying .tau.. The 
stretched exponential fit yields a valnp nf T.snh.p fnr parh valup of .tau.. The 
composite set of data point for t anh o {.tau.) is shown in FIG. 7 for three 
representative sandstone samples. We fit the resulting set of T . sub . 2 (.tau.) data 
to a single exponential function. FIG. 7 shows the nonlinear least squares fits as 
dotted curves. This allows one to extrapolate T . .sub ■ 2 (.tau.) to .tau.=0. T. sub. 2 
( . tau. .fwdarw. 0) should therefore give a value of the, trans/ye. roe, magnetization 
rplavahinn time t snh 7. independent of the strength of the internal magnetic field 
gradients and be proportional to some characteristic length of the pores. FIG. 8 
shows a plot of t .qnh 2 (.tau.=0) versus the characteristic pore diameter l.sub.c 
determined from the threghnl rf pressure for mercury injection via the Washburn 
equation. Such a relationship is of great value for a rapid order of magnitude 
determination of permeability based on the above mentioned relationship proposed by 
Katz and Thompson : k=l/126 1 . sub . c . sup . 2 . phi . . sup . 2 . And indeed as FIG . 9 shows 
there exists a correlation between the product of porosity .phi. and t wnh s> 
{. tau. . fwdarw. 0) and the independently determined permeability of the porous rock 
samples as predicted by the relation k .varies. [T.sub.2 (. tau. . fwdarw. 0) ]. sup. 2 
.phi.. sup. 2. In this equation the exponent nf T.mih.2 (. tau. . fwdarw. 0) is predicted 
to be 2. However the range of the data in FIG. 9 would indicate that the exponent 
falls into a range of about 1.8 to 2.2. T.wnh. s values obtained with the CPMG pulse 
sequence (which should be the method of choice) are strongly dependent on the pulse 
spacing .tau.. Only when the gradient and diffusion induced relaxation is taken into 
account, do the resulting T . suh . 2 (.tau.) values extrapolated to .tau.=0 show the 
discussed correlations with permeability and l.sub.c. 

Detailed Description Text (40) : 

For the restricted diffusion analysis of T. sub. 2 (.tau.) we used a nonlinear least 
squares fitting algorithm to fit the experimental T , sub , 2 (.tau.) data to the 
expression of equation 22. The pref actor of .tau. in the exponent of the exponential 
is one of the adjustable parameters, l.sub.nmr was directly determined from this 
adjustable fitting parameter obtained from the nonlinear least squares -fit. 

Detaile d DeKoription Text (41) : 

The present invention takes advantage of the correlations between T . nub . 2 determined 
with NMR and a characteristic pore or throat size which determines the jfluid flow 
permeability. This characteristic length l.sub.c can also be determined by mercury 
injection but mercury injection experiments cannot be carried out in a down-hole 
environment. For in situ nmr rfilayahinn time measurements in rock formations with an 
nmr logging tool it is feasible to carry out a t snh ? measurement described here 
with the present generation of logging tools and predict the f Inid flow permeability 
using NMR in a manner which is much less time-consuming than other methods known to 
date. 

Detailed Descr iption Paragraph Equation (2) : 
T . sub . 2 ( .tau. ) .varies. .DELTA. /F (. tau. ) (2) 

Detailed Description Paragraph Equation (14) : 

k vanpfl. T.Hiih.2 ( . tau . . f wdarw. 0 ) . sup . m .phi., sup. n (29) 

Detailed Description Paragraph Table ( 1) : 

TABLE 1 Sandstone samples and their porosity, 

permeability and l.sub.c deter- mined independently of the NMR measurements, 
sandstone type l.sub.c [.mu.m] porosity [%] permeability [md] 

Berea 14.6 20.5 273 Marsing No. 2 118.0 29.5 

54,000 Red Navajo 23.5 23.6 1138 Nugget 10.77 10.9 4.16 Silver No. 1 23.0 12.2 14.1 
Layered Navajo 28.5 25.1 883 Marsing No. 1 78.6 23.9 1276 Silver No. 2 79.2 30.2 
21,000 Table No. 2 91.6 24.1 3000 



CLAIMS : 
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1. A method for determining the permeability of porous media saturated with a liquid 
using nuclear magnsl:ir rpsnnanre OMR) and comprising: 

(a) applying a radiof requency pulse sequence which after an initial pulse generates 
successive spin echoes with a train of radio frequency pulses spaced apart by a time 
interval of length .tau. wherein all pulses have a carrier frequency corresponding 
to the Larmor frequency of the fluid spins filling the pore space of the medium for 
which the fin-id flow permeability is to be determined; 

(b) measuring the decay of the transverse, magnetization at each of the successive 
regularly spaced midpoints between the 180 degree pulses where the midpoints 
coincide with the peak of the spin echoes; 

(c) repeating steps a and b at least one more time wherein each repeat of step (a) 
uses said radio frequency pulse train with a different 3£alii£ of the pulse spacing 

. tau . ; 

(d) determining the transverse relaxation time T . sub . 2 (.tau.), from the transverse 
magnetization decay for each value of .tau., and determining one of a pref actor 
.DELTA, and a restricted diffusion length, l.sub.nmr, from said T.sub.2 {.tau.): 

(e) measuring the porosity of said porous media; 

(f ) determining the permeability of said media from the porosity and one of the 
prefactor .DELTA, and a restricted diffusion length, l.sub.nmr, from said T. sub . 2 
(.tau.) : . 

2. The method of claim 1 wherein said t ? (.tau.) is of a form which includes an 
exponential function of .tau.. 

3. The method of claim 2 wherein said t. snh o. (.tau.) is obtained from an expansion 
of the exponential functions. 

4. The method of claim 2 wherein said T . sub . 2 (.tau.) is a simple exponential 
function of .tau. and said prefactor .DELTA, is obtained by extrapolation of T . sub. 2 
( . tau . ) to . tau . . f wdarw. 0 . 

5. The method of claim 1 wherein said porous media is an earth formation and the 
measurement is performed in a well-bore environment with an NME logging tool. 

7. The method of claim 1 wherein the rplaxatinn time, T . sub . 2 is obtained from the 
franRvprsp magnetization decay in step (d) by relating the fransvprRP magnetization 
to time by a stretched exponential function. 

13. The method of claim 11 wherein prefactor .DELTA, is determined as prefactor of 
the observed t suh 2 ( . tau . ) dependence . 

14. The method of claim 4 wherein said permeability of step (f) is obtained from the 
product of the porosity and the extrapolated transversa relaxation time, t <=mh> 2 at 
.tau.=0 . 

15. The method of claim 14 wherein said porosity and the extrapolated transverse 
relaxation time, t snh 7 at .tau.=0 are each raised to a power. 
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DOCUMENT- IDENTIFIER: US 5370901 A 

TITLE: Compositions for increasing the image contrast in diagnostic investigations 
of the digestive tract of patients 



RHpf Summary Text (2) : 

The present invention concerns innocuous ingestible or enterally administrable 
compositions which, depending on the contrast agent incorporated thereto, can be 
used as contrast enhancer media for imaging, on the first hand in ultrasonic 
echography, and on the second hand, in nuclear magnetic resonance imaging (NMRI) , 
both of the gastro- intestinal tract of animal and human patients. 

Rrief Summary Text (5) : 

It is also well known that NMRI techniques comprise subjecting a patent to a main 
static magnetic field combined with a linear gradient magnetic field, both being 
directed to some parts of the body to be investigated. The magnetic fields act on 
the nuclei of atoms with fractional spin quantum numbers and encode them into 
various degrees of statistical alignment with different resonant frequencies in a 
selected direction of orientation; the nuclei of concern here are mainly that of 
hydrogen atoms, i.e. protons, these being predominantly that of molecules present in 
relatively high concentration in or around the organs to be investigated, viz, the 
protons of water and lipids. For doing the measurements, one will apply to the parts 
of the body under investigation pulses of radio -frequency that matches with the 
resonance energy of the protons involved in the tissues or f "hiidg of said parts of 
the body. When the protons under consideration are excited by a pulse of resonant 
energy, they are raised to a higher energy state which causes them to flip f,rom the 
avpragp orientation direction controlled by the magnetic field. Thereafter, the 
protons will return to their original state by relaxation in an exponential time 
dependent fashion, the corresponding energy then reemitted (spin- echo) forming a 
response signal typical of the protons under consideration, i.e. depending on their 
immediate environment. 

Brief S ummary Tpyf (6) : 

NMRI techniques are actually based on the detecting, acquiring and electronically 
processing of this signal {according to Fourier transforms) and thereafter 
displaying it spatially on a screen, thus forming an image whose various patterns 
correspond to areas having protons in different environments, i.e. to protons 
belonging to organ tissues or body fluids being subjected to investigation. 

Rr-ie>f Summary Tpyf (7) : 

Among the critical factors pertaining to MRT , one usually distinguishes two mutually 
perpendicular components of the proton-distinctive relaxation time parameter, namely 
the spin-lattice component along the axis of magnetization (called T.sub.l), which 
corresponds to the release of energy to the nuclear environment, and The 
perpendicular or transverse f ,qpi n-spi n) relaxation component (called T . nub . ?.) , that 
corresponds to the returning of the nucleus to the initial statistical energy level. 
Either T.sub.l or t anh ? can contribute to the definition of the NMR images 
depending on the kind of organ selected and the measurement conditions. 

Rrigf Summary Text (8) : 

It should be noted that when the measurements are carried out in the absence of 
agents added for increasing image contrasts, the differences in reTaxati nn time 
constants between protons in various parts of the organs are small and the image is 
of poor to bad quality. The contrast effect can however be enhanced by the presence, 
in the environment of the hydrated molecules under excitation, of a variety of 
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' magnetic species, e .g .^paramagnetic (which mainly af f ec^f . sub. 1) and ferromagnetic 
or superparamagnetic (which mainly affect the T . sub . 2 response) . The paramagnetic 
substances include some metals in the ionic or organo-metallic state (e.g. 
Fe.sup.+3, Mn.sup.+2, Gd.sup.+3 and the like, particularly in the form of chelates 
to decrease the intrinsic toxicity of the free metal ions) . Ferromagnetic contrast 
substances preferably include magnetic aggregate particles of micronic or 
submicronic size, i.e. not smaller than about 100-2 00 nm, for instance particles of 
magnetite (Fe.sub.3 O.sub.4), .gamma. -Fe . sub. 2 O.sub.3, ferrites and other magnetic 
mineral compounds of transition elements. Superparamagnetic materials are usually 
very small magnetic particles (below about 100-150 nm) which, because their size is 
under a critical value, do not behave any longer as small autonomous magnets, i.e 
they will align in a preferential direction only when subjected to an external 
magnetic field. The advantage of the superparamagnetic materials (also defined 
sometimes as superparamagnetic £lnida) over the ferromagnetic particles is mainly of 
efficiency density, i.e. being smaller, the number of available magnetic particles 
for a given wpi ght of metal is greater in the case of superparamagnetic particles 
than with ferromagnetic particles and the magnetic efficiency on the neighboring 
protons is further enhanced. 

Ftripf Summary Text (13) : 

It has been indicated above that paramagnetic contrast agents in which the metals 
are in the ionic state or in the form of metal -organic compounds are often 
metabolizable and toxic and, although this toxicity can be controlled to some extent 
by using very strong chelatants and non-metabolizable polymer carriers, it is 
desirable to further minimize possible hazards by using less toxic materials, e.g- 
non-metabolizable magnetic particles of sufficient size not to diffuse through the 
intestinal membrane; the micronic ferromagnetic and nanometric superparamagnetic 
aggregate particles typically fulfill such requirements. For instance, in U.S. Pat. 
No. 4,770,183 (ADVANCED MAGNETICS), there is recommended to use biodegradable 
sub-micron sized superparamagnetic metal oxide particles (1-50 nm) which may be used 
uncoated or coated with a polysaccharide (like dextran) or serum albumin. Coating is 
effected by precipitating the particles with alkali, starting with water solutions 
of metal salts in the presence of the polymer. These products are suitable for 
intravenous applications as well as for gastrointestinal applications, in which case 
they are administrable by intubation or enema, presumably because otherwise 
biodegradation by the stomach flnidR would be too fast and toxicity might become a 
problem. 

PT-if*f Summary Text (14) : 

In W085/04330 (NYCOMED) , there is disclosed the use of ferromagnetic particles as 
contrast agents for NMRI. As mentioned before, ferromagnetic particles are bigger 
than superparamagnetic particles and behave as small permanent magnets which also 
achieve a significant reduction of t wih ? For direct administration into the 
digestive tract, the ferromagnetic particles are preferably embedded in a cellulose 
matrix or coated with this matrix. Cellulose derivatives can also be added as 
viscosants but the reference indicates that contrast enhancement is not readily 
achieved beyond the stomach, presumably because the embedding cellulose matrix does 
not protect sufficiently the particles from attack by the stomach fluids . 
Non -biodegradable embedding or coating matrices are therefore recommended to 
minimize absorption of toxic materials by the body. 

Brief Summary Text (15) : 

EP-A-186.616 (SCHERING) discloses the use of complexes of particles of magnetite 
(Fe.sub.3 O.sub.4), .gamma. -iron oxide (Fe.sub.2 O.sub.3) and metal ferrites as 
contrast agents for NMRI. The cited complexants include oligo- and polysaccharides, 
proteins, polycarboxylic acids, protective colloids and other compounds- Examples of 
such compounds comprise polyvinyl -alcohol (PVA) , polysilanes, polyethylene -imine, 
dextran, dextrin, oleic acid, gelatin, globulin, albumin, insulin, peptides and 
antibodies. The particles can also be encapsulated in liposomes. For enteral 
administration, the contrast agents are suspended in a water medium which may 
contain further ingredients such as salt or excipients like methylcellulose, 
viscosants, lactose, mannitol and surfactants like lecithin, Tween . RTM . , Myrj.RTM. 
and the like. 

Brief S ummary Tpxt (18) : 

After testing the compositions of the cited prior art consisting of coated or 
uncoated magnetic particles in admixtures with polymer carriers, the present 
inventors noted that the contrast effect in nmr imaging is generally unstable and 
rapidly vanishes, presumably because despite the presence of the carrier phase the 
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1 magnetic particles tenc^to coalesce or coacervate togetrW: under the influence of 
the external magnetic field which strongly reduces their controlling effect on the 
.qpi n-rplayatinn of the neighbouring protons. They however also found that such 
undesirable coalescence of the magnetic particles and vanishing of the T ■ mih . 2 
relaying effect can be prevented by either selecting as the carrier phase 
substantially_water-insoluble hydrophilic water- swelling substrates which tend to 
form gels with water or, when using water-soluble polymer carriers, raising the pH 
of the aqueous medium containing the magnetic particles to at least 13 when admixing 
with the polymer solution, and keeping the dry weight ratio of said polymers to 
magnetic particles not below 5:1 and, preferably, in the range of 100:1-10:1. 

Brief Summary Text (19) : 

Although the exact reason of these findings is not definitely explained, it can be 
postulated that using carrier matrices which form nearly insoluble gels upon 
admixing with water (thixotropic or pseudo-plastic solutions) will locally raise the 
viscosity at the particle/carrier interface to such extent that the particle 
mobility is impeded and agglomeration is prevented. 

Rrief Summary Text (22) : 

The foregoing findings form the basis of the present invention as recited below and 
in the annexed claims, i.e. the invention mainly concerns diagnostic compositions to 
be used as contrast enhancing agents, on a first hand in the imaging by echography, 
and on a second hand in the nmr imaging, of portions of the digestive tract of 
humans and animals, said compositions comprising, depending on the needs, i.e. on 
the kind of investigation, at least one member of a group consisting of two types of 
contrast enhancing particles, viz. magnetically and ultrasonically responsive 
particles, in admixture with, or chemically bonded to, at least one physiologically 
acceptable carrier phase The said carrier phase is substantially water insoluble, 
although swellable by hydration and gel forming, and, when hydrated, has 
differential affinity for the gastro- intestinal mucosa, i.e. it adheres preferably 
to some areas thereof, thus forming luminal linings or coatings having specific 
magnetic contrast response patterns which enable improved visualization of specific 
portions of the oeso-gastro-duodenal tract. 

RH pf Summary Teyf (23) : 

The carrier phase of the present composition can be selected from water-soluble 
polymers which can form homogeneous solutions in aqueous media, in which case the 
proportion of signal generating particles therein preferably not exceeds 10% by 
wpi ght of solution, more preferable 5%, for optimalized stability against 
coalescence when subjected to a steady magnetic field. The polymers which are 
convenient are for instance dextran, polyvinyl -pyrrol idone (PVP) , 

carboxymethyl- cellulose (CMC) and the like. It is quite surprising that, contrary to 
products of the prior art with a high particle/carrier weight ratio, homogeneous 
dispersion of magnetic particles in relatively low concentrations in the carrier 
polymer remain stable in a magnetic f i eld with virtually no or very little 
agglomeration of the magnetic particles, and hence no substantial weakening of the 
imaging response under the conditions of nmt? imaging. Although the compositions with 
water-soluble polymers have no recognized specific affinity for the internal mucosa 
membrane of the digestive tract, they are useful for controlling the transit 
therethrough of the magnetic particles, either alone or when preferably used in 
association with bioadhesive carrier phases. 

Brief Summary Text (25) : 

These carriers have the general property to get hydrated with wat.ftr and form 
therewith substantially insoluble gels which strongly adhere to the membrane mucosa. 
Depending on the selected carrier phase, the compositions of the invention can be 
made to have controllable transit time through the gastrointestinal tract and 
particular adhesion to selected portions thereof, thus enabling efficient NMRI of 
selected organs. Normally, the compositions using polymer carriers having no or 
little bioadhesion will permit faster transit than the compositions using carriers 
with strong bioadhesive character, hence controlled transit can be achieved with 
compositions in which non- bioadhesive and bioadhesive carriers are admixed in 
suitable proportions. It has also been surprisingly found that the particular 
selection of carriers used in the compositions of the invention, i.e. the carriers 
with enhanced affinity for mucosa membranes at low pH and/or soluble carriers with 
high dispersive capacity for magnetic particles, provide generally excellent 
t-ransvprsp and longitudinal luminal imaging homogeneity. 

Brief Summary Text (26) : 
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One preferred type of THToadhesive polymers to be used a^the carrier phase in this 
invention is a polymer or copolymer of acrylic acid containing a proportion of other 
monomers (generally di- or polyfunctional allyl ethers or acrylates) to impart a 
degree oil_watex- insolubility and swelling capacity to form gels. Suitable polymers 
of this type are available on the market under the names ail Carbopol . RTM . or 
Polycarbophil . RTM . , for instance from the Goodrich Company. Other suitable polymers 
of comparable type are disclosed in EP-A-309.404 (included for reference) and 
comprise copolymers of acrylic acid with allyl ethers of glycols or sugars. Upon 
addition of water these polymers will form viscous dispersions of microgels which 
have strong affinity for internal mucous membranes. Furthermore, the gelling and 
swelling properties of such polymers are pH dependent; hence the volume, bulk and 
adhesive properties of the carrier can be controlled by adjusting the pH to a 
desired value - 

Brief Summary Text (27) : 

The present ingestible NMRI contrast compositions are stable even at pH 1 and this 
is another advantage over similar products of the prior art as they are not attacked 
by stomach f liHd« which normally dissolve magnetite particles into Fe.sup.+2 and 
Fe.sup.+3 ions. This is obviously undesirable for both toxicity and imaging 
efficiency reasons. Bioadhesivity is particularly effective under acid or near 
neutral conditions; at high pH the bioadhesivity is decreased and may become 
negligible. Hence control of bioadhesivity by pH control is another asset of the 
present invention. 

RriPf Summary Tp.xt (28) : 

For manufacturing the contrast compositions according to the invention one usually 
admixes the signal generating particles (the echogenic or magnetically responsive 
particles, or both) with the carrier phase. In regard to the MRX aspect of the 
invention, one may use magnetic particles which can be ferromagnetic or 
superparamagnetic. The nature of both and the distinction between them has been 
stressed before in this specification. Hence, one prefers superparamagnetic 
particles for opmimalized density efficiency, i.e. magnetic particles of size not 
exceeding about 150 nm and preferably in the 1-100 nm range (10-1000 A) . The 
preparation of such magnetic particles is well known and disclosed in many 
references, e.g. the aforementioned references and further references, e.g. U.S. 
Pat. No. 4,554,088; EP-A-125.995 (Advanced Magnetics) ; U.S. Pat. No. 4,267,234; U.S. 
Pat. No. 4,157,323 (California Institute); WO-A-78/00005 (Mosbach) ; WO-A-83/01738 ; 
WO-A-83/03426; WO-A-84/00294 (Schroder) . They are also available commercially, e.g. 
from Ferrof luidics Corp., Burlington, Mass. They normally result from the 
alkalinization of an aqueous solution containing Fe.sup.+2 and Fe.sup.+3 salts in 
correct proportions. 

Brief Summary Text (30) : 

In one embodiment of a method for preparing the compositions according to the 
present invention, an aqueous suspension of the magnetic particles is treated with 
alkali to raise the pH to 13 or more, then the obtained alkaline suspension is 
admixed with an aqueous solution or dispersion of the carrier phase at the same pH 
and finally a_water- compatible organic solvent in which the carrier phase is 
insoluble is added, whereby a precipitate of the desired composition is formed. The 
composition can thereafter be separated, for instance by filtration, and washed free 
from alkali, then it can be stored dry. When used for internal NMRI investigation, 
the composition can be rehydrated with wafpr or an aqueous solution of innocuous 
excipients, this being for oral or enteral administration. 

Brief Summary Text (31) : 

In the foregoing embodiment, the magnetic particles are not covalently linked to the 
carrier phase, the latter being Mater -soluble, but still they are immobilized 
therein presumably due to the existence of electrostatic or Van der Waals forces 
between the carrier phase constituents and the magnetic parricides or, if the 
carrier is insoluble in wat-pr but in hydrated gel form, to immobilization of the 
hydrated magnetic particles within the carrier gelled structure. 

Rr-ipf Summary Tpvt-. (34) : 

For using the present compositions as an echogenic or NMRI contrast agent (or both) 
in the study of the digestive tract, the compositions (when stored dry) are admixed 
with an aqueous phase suitable for oral administration in order to provide a gel 
which will adhere to a predetermined extent to the mucous membrane of the digestive 
tract and therefore will carry the particles through said tract to the portion 
thereof to be visualized. The rate of transfer and the transit time can be adapted 
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- by properly selecting nature and degree of bioadhesHh of the carrier phase. For 
this, variable proportions of contrast compositions including water-soluble carrier 
polymer can be admixed with the bioadhesive carrier, these proportions being 
selected to reach up to 99% by wpi ght or more in the case where little restricted 
transit, unrestricted transit or accelerated transit is desired. For instance, for 
NMR investigations, using mainly dextran as the carrier phase in rats gave a transit 
rate similar to that observed with no carrier, while using mainly CMC as carrier 
gave accelerated transit. Using Carbopol .RTM. as the carrier phase gave strongly 
retarded transit, particularly in the duodenum portion of the intestinal tract. The 
present contrast compositions may also advantageously comprise isoosmolarity agents 
which minimize loss or gain of_water during intestinal transit, i.e. diffusion by 
osmosis; such agents may comprise carbohydrates such as sorbitol, mannitol, xylitol 
and the like. 

Detailed Descr iptinn Text (3) : 

A water solution (30 ml) of FeCl. sub. 3 .6H. sub. 2 O (185.3 mg; 0.685 mmol) and 
FeCl. sub. 2. 4H. sub. 2 O (80.5 mg; 0.405 mmol) labeled with .sup. 59 Fe was brought 
slowly to pH 13.4 by the dropwise addition of 3% aqueous NH.sub.4 solution; then it 
was heated to 75. degree. C. for about 10 min. This resulted in the formation of a 
dark brown suspension of magnetite in large and coarse particles. 

Detailed Description Text (4) : 

To this was added under agitation a solution of 5 g of Dextran in 200 ml of wahar 
and the pH was reajusted to 13.4 by adding some more NaOH solution. After 15 min 
more stirring, alcohol was added, whereby a very fine precipitate of 
Dextran-magnetite was formed. This precipitate was drained under suction, washed 
with alcohol and dried at 50. degree. C. in air. 

Detailed Description Text. (6) : 

The proton magnet ir rpsnnanrp frarifivPrsP rplayivi fy. expressed as R.Sub.2, i.e. 

R.sub.2 =i /t .giih. o I mmol Fe) was measured on aqueous dispersions of the 
aforementioned precipitate by means of a 60 MHz RMN RppntrnniRtp.r. The measurements 
were effected at mid-height of the absorption peaks. Values of R.sub.2 in the range 
of 500-550 (mM. s) . sup. -1 were obtained which did not substantially change with time, 
thus showing that reagglomeration of the particles did not occur. 

Dp1-ai1pd npscr ipt-ion Text (9) : 

To a solution of Dextrin (80 g) in 180 ml of H.sub.2 O at 70. degree. C. were added a 
solution of 19 g of FeCl . sub. 3 . 6H. sub. 2 O (70 mmol) and 7.01 g FeCl . sub . 2 . 4H. sub . 2 O 
(35 mmol) in 105 ml of water. A tracer quantity of .sup. 59 FeCl. sub. 3 was also 
added. The pH was raised to 2.4- by first adding 20% aqueous Na.sub.2 CO. sub .-3, then 
it was brought to pH 11 with 10N NaOH added dropwise, whereby the color of iron 
oxides developed. The mixture was refluxed for 30 min, then it was cooled to room 
temperature and the pH was brought to 6.2 with 6N HCl. Then 0.8 liter of ethanol was 
added which caused the formation of a precipitate. This was collected, resuspended 
in watpr and the suspension was dialyzed for 24 hrs in running water (membrane 
mf -of f-i n ' nnn) . Then the suspension was freeze-dried which gave a blackish powder. 

TlPta i 1 f*d Dpsrri pr. i on Text (10) : 

Measurements of R.sub.2 as described above gave none or useless responses. Analogous 
results were obtained by replacing Dextrin by a comparable amount of Dextran. 
Obviously, the ratio of metal to polymer in this formulation is much too high to 
prevent reagglomeration of the particles subjected to a magnetic field. 

Detailed Description Text (13) : 

To a solution of 13.56 g FeCl . sub . 3 . 6H . sub . 2 O (50 mmol) and 5.3 g 
FeCl. sub. 2. 4H. sub. 2 O (27 mmol) in 60 ml of water labeled with .sup. 59 Fe (about 
10. sup. 6 cpm) were added dropwise 60 ml of 5M aqueous NaOH until the pH was above 
13 . The solution was left to stand until a blackish sediment accumulated at the 
bottom of the container. This was washed carefully by decantation ten times with 
portions of 300 ml H.sub.2 0. After the last portions, the pH of the water had 
dropped to 9.3. 

Detailed Dencriprinn Text (15) : 

Sixty ml of glycerol were added and the mixture was heated under reduced pressure in 
a rotavapor until a black viscous residue has formed. This residue was diluted with 
240 ml of water r agitated for 20 min and allowed to rest whereby it separated into 
two phases. This was dialyzed for 24 hrs against a citrate buffer (10 mM, pH 8 . 3 
with ammonia) in order to remove low molecular material (glycerol, unreacted silane 
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and acrylate oligomer: 



Detail f* d Descri ption Text (16) : 

A suspension was made containing 20 ml of acrylic acid 10 ml H.sub.2 O and 10 ml of 
the silanized magnetite prepared as described above. This suspension was heated to 
50. degree. C. and a 10% aqueous ammonium persulfate solution was added dropwise. 
After the polymerization was complete, the polymer was ground in 500 ml of watf*T- and 
dialyzed against running water . Then it was freeze-dried to give 20.1 g of silvery 
powder. 

Detailed Description Text (17) : 

A suspension of 1 g of this powder in 100 ml of water was prepared and 1 ml of this 
suspension was mixed with 9 ml of a 1% aqueous solution of Carbopol R at DH 7-8. The 
rplaYflfinn time t gnh s> was measured as described in Example 1, the valnp found 
being R.sub.2 =160 (s . mM) . sup . -1 . 

Detailed Description Text (20) : 

The Table provides the names of the polymers, the yield of the preparation 
(calculated on the basis of the iron converted to magnetite) and the relaxi vi f.y 
R . sub . 2 in terms of s . sup . - 1 mmol 

Detailed Description Text (24) : 

Samples of the suspensions (0.1 ml) were diluted with 9.9 ml (dilution 1:100) of 1% 
by weight aqueous polymer solutions also containing 0.3 mol/1 of mannitol to 
preserve isotonicity. The polymers used were: 

npt-ailftd nftfirriptinn Tsxt (25) : 

a) Carboxymethylcellulose (CMC) Na salt (0.1 g/9.9 ml of water). 
HPl-ailpd rtefirriH pt i on Tpyt (29) : 

e) Control; this was made as above by diluting the magnetite preparation to 1/100 
but using only watp-r without polymer. 

Detailed Description Text (30) : 

Five ml of samples a) to e) were administered to Sprague-Dawly rats (about 180 g) 
which were kept fast except for water 24 hrs before administration. The samples were 
introduced intragastrically with a ball-point syringe. The tracer amount in each 
sample accounted for about 10. sup. 6 cpm. 

Detailed Description Text (36) : 

A 3% aqueous solution of NH.sub.4 OH -(about 25 ml) was slowly added dropwise under 
agitation to a 100 ml water solution of 619.2 mg of FeCl . sub. 3 . 6H. sub. 2 O (2.29 
mmol) and 262.1 mg of FeCl . sub. 2 . 4H. sub. 2 O (1.32 mmol) . The mixture was then heated 
to 75. degree. C. for 10 min. There were thereafter added 200 ml of an aqueous 2.5% 
by weight solution of sodium alginate and the pH was raised to 13.5 with 40% aqueous 
NaOH. After further agitating for 15 min at room temperature, alcohol was added to 
precipitate the magnetite -polymer composition; the powder was drained on a filter, 
washed with alcohol and dried; the yield was 4.96 g containing 0.216 mmol of Fe 
(theory=0.22 mmol) . The mmolar t-r-angy^-rgf* relativity R.sub.2 was measured at 60 MHz 
on a water- suspension of the powder with a perkin-Elmer/Hitachi nmr apparatus. The 
measured valnp was R.sub.2 =465 (s.mM) .sup.- 1. 

DPf-a-n^d nPRrripfinn Tpyf (39) : 

The magnetite suspension was allowed to come back to room temperature, the pH was 
lnwsrp.fi to 6.5 with IN HC1 and there were added two grams of 

aluminum- sucrose-octasulf ate (sucralfate) sold under the name of Keal.RTM. by 
Laboratories SINBIO, 75116 Paris (France) . The mixture was diluted with 50 ml 
H.sub.2 O and further agitated for 10 min after which it was left to stand, whereby 
a precipitate of magnetite- sucralfate settled at the bottom of the vessel. This was 
drained on a filter, washed wit h water and dried, yield 2.03 g. R.sub.2 (in 
water) =283 (s.mM) .sup. -1. 

Detailed Description Text (43) : 

Montmorillonite clay of the Beidel type ( BEDELIX . RTM . ) available from Laboratories 
Beaufour, 28100 Dreux (France) was tested for NMRI contrast response under the 
following conditions: three g of BEDELIX. RTM. were stirred in 100 ml of a 0.5% 
solution of Carbopol. RTM. . The l/T.sub? value of this preparation was measured as 
usual with a Perkin-Elmer/Hitachi NMR apparatus at 60 MHz and found to be 435.5 
s.sup.-l without the addition of magnetite. This value was stable (due to the 
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' * presence of the Carbopol . RTM . ) and the preparation couJ^^e used as such for the 
NMRI of the digestive tract of animals. 

Detailed DPRcript-.-inn Text (45) : 

A suspension of magnetite was prepared as disclosed in Example 1, first paragraph. A 
0.2 ml portion of this suspension was added to 99.8 ml of a 0.5% by weight solution 
of Carbopol .RTM. . The iron concentration in this preparation was 0.178 mmol/ml; the 
1 /T.finh 9 valnp was 144.5 s.sup.-l, i.e. R.sub.2 =812 ( S . mM) . sup . - 1 . 

Detailed Descri ption Text (46) : 

Aliquots of the magnetite preparation were taken and to each were added a different 
quantity of BEDELIX.RTM. . The R.sub.2 values were measured as usual and reported 
below in function to the % by weight, of added BEDELIX.RTM. . 

Dp tailed Description Text (49) : 

Male Sprague Dawley rats were used for the experiment after 24 hours of fasting. 
Eight ml of the Carbopol -magnetite solution were administrated orally to the 
conscious animals which were anesthetized 10 min later by 30 ml/kg of Pentotal . 
Transversal MR images were then taken on an Esaote ESATOM 5000 imager equipped with 
a special 8 cm i.d. RF receiver coil, using a slice of 2 mm, a F.O.V. of 
15. 9. times. 15.9 cm, and a matrix size of 128 . times . 256 . T.sub.2 weighted images (SE 
200/70/1), were acquired as well as T.sub.l weighted images (SE 500/16/2), and 
intermediate scan images {SE 350/50/2) . For the T2 weighted and intermediate scan 
images gradient moment hulling techniques were used in order to minimize respiratory 
artifacts. The T2 wfnghtpd and intermediate scan images show a clear delineation of 
the darkened and expanded bowel loops. Especially the 350/50/2 image showed very 
clearly the single loops of the small intestine. The wall of the loops could be 
clearly observed. A cross section of the colon and of a kidney were also seen as 
well as abdominal and dorsal muscles. The contrast media was distributed evenly over 
the whole GI tract. 

Detaile d Description Text (51) : 

A suspension of echogenic microballoons was prepared as described in Example 4 of 
reference EP-A-458 745 using poly-L- lactic acid (commercially available under the 
name of Resomer.RTM. R-207 from Boeringer Ingelheim, Germany) . There was obtained a 
suspension of microballoons in distilled water (concentration 10. sup. 9 /ml; average 
size 5.2 .rau.m) . 

Ppt-aileri Description Text (52) : 

On the other hand, an aqueous suspension of gelled bioadhesive polymer was prepared 
-by dispersing 1& b y weight- of Carbopol .RTM. resin (Goodrich Company) in watex. The 
pH of this suspension was brought to 13 with concentrated NaOH and to 100 ml of the 
alkalinized solution were added 10 ml of the aforementioned microballoon suspension. 
After homogenizing the mixture under agitation, a quantity of ethanol was added 
sufficient to precipitate the solids; the polymer was drained on a filter under 
suction, washed with alcohol and dried under vacuum. The dry powder was then 
resuspended in 100 ml of a 0.3 M mannitol aqueous solution and the pH adjusted to 



3-4 with HC1. 



DPt-ail^d npRnriphion Paragraph Tahl e (3) : 



R.sub.2 (s .multidot. mM).sup.-l 
918 1 251.3 1412 1.5 402.1 2259 



% Of BEDELIX .RTM. 1 /T.mih.2 (s.sup.-l) 

0.5 163.4 
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